-
Tuning the exchange bias on a single atom from 1 mT to 10 T
Authors:
Kai Yang,
William Paul,
Fabian D. Natterer,
Jose L. Lado,
Yujeong Bae,
Philip Willke,
Taeyoung Choi,
Alejandro Ferrón,
Joaquín Fernández-Rossier,
Andreas J. Heinrich,
Christopher P. Lutz
Abstract:
Shrinking spintronic devices to the nanoscale ultimately requires localized control of individual atomic magnetic moments. At these length scales, the exchange interaction plays important roles, such as in the stabilization of spin-quantization axes, the production of spin frustration, and creation of magnetic ordering. Here, we demonstrate the precise control of the exchange bias experienced by a…
▽ More
Shrinking spintronic devices to the nanoscale ultimately requires localized control of individual atomic magnetic moments. At these length scales, the exchange interaction plays important roles, such as in the stabilization of spin-quantization axes, the production of spin frustration, and creation of magnetic ordering. Here, we demonstrate the precise control of the exchange bias experienced by a single atom on a surface, covering an energy range of four orders of magnitude. The exchange interaction is continuously tunable from milli-eV to micro-eV by adjusting the separation between a spin-1/2 atom on a surface and the magnetic tip of a scanning tunneling microscope (STM). We seamlessly combine inelastic electron tunneling spectroscopy (IETS) and electron spin resonance (ESR) to map out the different energy scales. This control of exchange bias over a wide span of energies provides versatile control of spin states, with applications ranging from precise tuning of quantum state properties, to strong exchange bias for local spin doping. In addition we show that a time-varying exchange interaction generates a localized AC magnetic field that resonantly drives the surface spin. The static and dynamic control of the exchange interaction at the atomic-scale provides a new tool to tune the quantum states of coupled-spin systems.
△ Less
Submitted 7 June, 2019;
originally announced June 2019.
-
Engineering the eigenstates of coupled spin-1/2 atoms on a surface
Authors:
Kai Yang,
Yujeong Bae,
William Paul,
Fabian D. Natterer,
Philip Willke,
Jose L. Lado,
Alejandro Ferrón,
Taeyoung Choi,
Joaquín Fernández-Rossier,
Andreas J. Heinrich,
Christopher P. Lutz
Abstract:
Quantum spin networks having engineered geometries and interactions are eagerly pursued for quantum simulation and access to emergent quantum phenomena such as spin liquids. Spin-1/2 centers are particularly desirable because they readily manifest coherent quantum fluctuations. Here we introduce a controllable spin-1/2 architecture consisting of titanium atoms on a magnesium oxide surface. We tail…
▽ More
Quantum spin networks having engineered geometries and interactions are eagerly pursued for quantum simulation and access to emergent quantum phenomena such as spin liquids. Spin-1/2 centers are particularly desirable because they readily manifest coherent quantum fluctuations. Here we introduce a controllable spin-1/2 architecture consisting of titanium atoms on a magnesium oxide surface. We tailor the spin interactions by atomic-precision positioning using a scanning tunneling microscope (STM), and subsequently perform electron spin resonance (ESR) on individual atoms to drive transitions into and out of quantum eigenstates of the coupled-spin system. Interactions between the atoms are mapped over a range of distances extending from highly anisotropic dipole coupling, to strong exchange coupling. The local magnetic field of the magnetic STM tip serves to precisely tune the superposition states of a pair of spins. The precise control of the spin-spin interactions and ability to probe the states of the coupled-spin network by addressing individual spins will enable exploration of quantum many-body systems based on networks of spin-1/2 atoms on surfaces.
△ Less
Submitted 11 December, 2017;
originally announced December 2017.
-
Reading and Writing Single-Atom Magnets
Authors:
Fabian D. Natterer,
Kai Yang,
William Paul,
Philip Willke,
Taeyoung Choi,
Thomas Greber,
Andreas J. Heinrich,
Christopher P. Lutz
Abstract:
The highest-density magnetic storage media will code data in single-atom bits. To date, the smallest individually addressable bistable magnetic bits on surfaces consist of 5-12 atoms. Long magnetic relaxation times were demonstrated in molecular magnets containing one lanthanide atom, and recently in ensembles of single holmium (Ho) atoms supported on magnesium oxide (MgO). Those experiments indic…
▽ More
The highest-density magnetic storage media will code data in single-atom bits. To date, the smallest individually addressable bistable magnetic bits on surfaces consist of 5-12 atoms. Long magnetic relaxation times were demonstrated in molecular magnets containing one lanthanide atom, and recently in ensembles of single holmium (Ho) atoms supported on magnesium oxide (MgO). Those experiments indicated the possibility for data storage at the fundamental limit, but it remained unclear how to access the individual magnetic centers. Here we demonstrate the reading and writing of individual Ho atoms on MgO, and show that they independently retain their magnetic information over many hours. We read the Ho states by tunnel magnetoresistance and write with current pulses using a scanning tunneling microscope. The magnetic origin of the long-lived states is confirmed by single-atom electron paramagnetic resonance (EPR) on a nearby Fe sensor atom, which shows that Ho has a large out-of-plane moment of $(10.1 \pm 0.1)$ $μ_{\rm B}$ on this surface. In order to demonstrate independent reading and writing, we built an atomic scale structure with two Ho bits to which we write the four possible states and which we read out remotely by EPR. The high magnetic stability combined with electrical reading and writing shows that single-atom magnetic memory is possible.
△ Less
Submitted 13 July, 2016;
originally announced July 2016.
-
Interplay between Orbital Magnetic Moment and Crystal Field Symmetry: Fe atoms on MgO
Authors:
S. Baumann,
F. Donati,
S. Stepanow,
S. Rusponi,
W. Paul,
S. Gangopadhyay,
I. G. Rau,
G. E. Pacchioni,
L. Gragnaniello,
M. Pivetta,
J. Dreiser,
C. Piamonteze,
C. P. Lutz,
R. M. Macfarlane,
B. A. Jones,
P. Gambardella,
A. J. Heinrich,
H. Brune
Abstract:
We combine density functional theory, x-ray magnetic circular dichroism, multiplet calculations, and scanning tunneling spectroscopy to assess the magnetic properties of Fe atoms adsorbed on a thin layer of MgO(100) on Ag(100). Despite the strong axial field due to the O ligand, the weak cubic field induced by the four-fold coordination to Mg atoms entirely quenches the first order orbital moment.…
▽ More
We combine density functional theory, x-ray magnetic circular dichroism, multiplet calculations, and scanning tunneling spectroscopy to assess the magnetic properties of Fe atoms adsorbed on a thin layer of MgO(100) on Ag(100). Despite the strong axial field due to the O ligand, the weak cubic field induced by the four-fold coordination to Mg atoms entirely quenches the first order orbital moment. This is in marked contrast to Co, which has an out-of-plane orbital moment of $L_z = \pm 3$ that is protected from mixing in a cubic ligand field. The spin-orbit interaction restores a large fraction of the Fe orbital moment leading a zero-field splitting of $14.0 \pm 0.3$~meV, the largest value reported for surface adsorbed Fe atoms.
△ Less
Submitted 25 June, 2015;
originally announced June 2015.
-
Transient adhesion and conductance phenomena in initial nanoscale mechanical contacts between dissimilar metals
Authors:
William Paul,
David Oliver,
Yoichi Miyahara,
Peter Grütter
Abstract:
We report on transient adhesion and conductance phenomena associated with tip wetting in mechanical contacts produced by the indentation of a clean W(111) tip into a Au(111) surface. A combination of atomic force microscopy and scanning tunneling microscopy was used to carry out indentation and to image residual impressions in ultra-high vacuum. The ~7 nm radii tips used in these experiments were…
▽ More
We report on transient adhesion and conductance phenomena associated with tip wetting in mechanical contacts produced by the indentation of a clean W(111) tip into a Au(111) surface. A combination of atomic force microscopy and scanning tunneling microscopy was used to carry out indentation and to image residual impressions in ultra-high vacuum. The ~7 nm radii tips used in these experiments were prepared and characterized by field ion microscopy in the same instrument. The very first indentations of the tungsten tips show larger conductance and pull-off adhesive forces than subsequent indentations. After ~30 indentations to a depth of ~1.7 nm, the maximum conductance and adhesion forces reach steady-state values approximately 12x and 6x smaller than their initial value. Indentation of W(111) tips into Cu(100) was also performed to investigate the universality of tip wetting phenomena with a different substrate. We propose a model from contact mechanics considerations which quantitatively reproduces the observed decay rate of the conductance and adhesion drops with a 1/e decay constant of 9-14 indentation cycles. The results show that the surface composition of an indenting tip plays an important role in defining the mechanical and electrical properties of indentation contacts.
△ Less
Submitted 7 August, 2013; v1 submitted 5 August, 2013;
originally announced August 2013.
-
FIM tips in SPM: Apex orientation and temperature considerations on atom transfer and diffusion
Authors:
William Paul,
David Oliver,
Yoichi Miyahara,
Peter Grütter
Abstract:
Atoms transferred to W(111) and W(110) tip apices from the Au(111) surface during tunneling and approach to mechanical contact experiments in STM are characterized in FIM at room temperature and at 158 K. The different activation energies for diffusion on the (111) and (110) tip planes and the experiment temperature are shown to be important factors controlling the extent of changes to the atomic…
▽ More
Atoms transferred to W(111) and W(110) tip apices from the Au(111) surface during tunneling and approach to mechanical contact experiments in STM are characterized in FIM at room temperature and at 158 K. The different activation energies for diffusion on the (111) and (110) tip planes and the experiment temperature are shown to be important factors controlling the extent of changes to the atomic structure of the tip. W(111) tips are much better suited to scanning probe studies which require the characterization of an atomically defined tip and subsequent verification of its integrity in FIM. The statistics of the observed spikes in the tunneling current when the tips are approached to Au(111) are interpreted using a simple model of adatoms diffusing through the STM junction.
△ Less
Submitted 3 August, 2013;
originally announced August 2013.
-
Comment on 'Field ion microscopy characterized tips in noncontact atomic force microscopy: Quantification of long-range force interactions'
Authors:
William Paul,
Peter Grütter
Abstract:
A recent article by Falter et al. (Phys. Rev. B 87, 115412 (2013)) presents experimental results using field ion microscopy characterized tips in noncontact atomic force microscopy in order to characterize electrostatic and van der Waals long range forces. In the article, the tip radius was substantially underestimated at ~4.7 nm rather than ~8.1 nm due to subtleties in the application of the ring…
▽ More
A recent article by Falter et al. (Phys. Rev. B 87, 115412 (2013)) presents experimental results using field ion microscopy characterized tips in noncontact atomic force microscopy in order to characterize electrostatic and van der Waals long range forces. In the article, the tip radius was substantially underestimated at ~4.7 nm rather than ~8.1 nm due to subtleties in the application of the ring counting method. We point out where common errors in ring counting arise in order to benefit future experimental work in which the determination of tip radius by FIM is important.
△ Less
Submitted 1 August, 2013; v1 submitted 24 April, 2013;
originally announced April 2013.
-
Estimation of Persistence Lengths of Semiflexible Polymers: Insight from Simulations
Authors:
Hsiao-Ping Hsu,
Wolfgang Paul,
Kurt Binder
Abstract:
The persistence length of macromolecules is one of their basic characteristics, describing their intrinsic local stiffness. However, it is difficult to extract this length from physical properties of the polymers, different recipes may give answers that disagree with each other. Monte Carlo simulations are used to elucidate this problem, giving a comparative discussion of two lattice models, the s…
▽ More
The persistence length of macromolecules is one of their basic characteristics, describing their intrinsic local stiffness. However, it is difficult to extract this length from physical properties of the polymers, different recipes may give answers that disagree with each other. Monte Carlo simulations are used to elucidate this problem, giving a comparative discussion of two lattice models, the self-avoiding walk model extended by a bond bending energy, and bottle-brush polymers described by the bond fluctuation model. The conditions are discussed under which a description of such macromolecules by Kratky-Porod worm-like chains holds, and the question to what extent the persistence length depends on external conditions (such as solvent quality) is considered. The scattering function of semiflexible polymers is discussed in detail, a comparison to various analytic treatments is given, and an outlook to experimental work is presented.
△ Less
Submitted 11 March, 2013;
originally announced March 2013.
-
Scattering function of semiflexible polymer chains under good solvent conditions
Authors:
Hsiao-Ping Hsu,
Wolfgang Paul,
Kurt Binder
Abstract:
Using the pruned-enriched Rosenbluth Monte Carlo algorithm, the scattering functions of semiflexible macromolecules in dilute solution under good solvent conditions are estimated both in $d=2$ and $d=3$ dimensions, considering also the effect of stretching forces. Using self-avoiding walks of up to $N = 25600$ steps on the square and simple cubic lattices, variable chain stiffness is modeled by in…
▽ More
Using the pruned-enriched Rosenbluth Monte Carlo algorithm, the scattering functions of semiflexible macromolecules in dilute solution under good solvent conditions are estimated both in $d=2$ and $d=3$ dimensions, considering also the effect of stretching forces. Using self-avoiding walks of up to $N = 25600$ steps on the square and simple cubic lattices, variable chain stiffness is modeled by introducing an energy penalty $ε_b$ for chain bending; varying $q_b=\exp (- ε_b/k_BT)$ from $q_b=1$ (completely flexible chains) to $q_b = 0.005$, the persistence length can be varied over two orders of magnitude. For unstretched semiflexible chains we test the applicability of the Kratky-Porod worm-like chain model to describe the scattering function, and discuss methods for extracting persistence length estimates from scattering. While in $d=2$ the direct crossover from rod-like chains to self-avoiding walks invalidates the Kratky-Porod description, it holds in $d=3$ for stiff chains if the number of Kuhn segments $n_K$ does not exceed a limiting value $n^*_K$ (which depends on the persistence length). For stretched chains, the Pincus blob size enters as a further characteristic length scale. The anisotropy of the scattering is well described by the modified Debye function, if the actual observed chain extension $<X>$ (end-to-end distance in the direction of the force) as well as the corresponding longitudinal and transverse linear dimensions $<X^2> - <X>^2$, $<R_{g,\bot}^2>$ are used.
△ Less
Submitted 17 August, 2012;
originally announced August 2012.
-
Implementation of atomically defined Field Ion Microscopy tips in Scanning Probe Microscopy
Authors:
William Paul,
Yoichi Miyahara,
Peter Grütter
Abstract:
The Field Ion Microscope (FIM) can be used to characterize the atomic configuration of the apex of sharp tips. These tips are well suited for Scanning Probe Microscopy (SPM) since they predetermine SPM resolution and electronic structure for spectroscopy. A protocol is proposed to preserve the atomic structure of the tip apex from etching due to gas impurities during the transfer period from FIM t…
▽ More
The Field Ion Microscope (FIM) can be used to characterize the atomic configuration of the apex of sharp tips. These tips are well suited for Scanning Probe Microscopy (SPM) since they predetermine SPM resolution and electronic structure for spectroscopy. A protocol is proposed to preserve the atomic structure of the tip apex from etching due to gas impurities during the transfer period from FIM to SPM, and estimations are made regarding the time limitations of such an experiment due to contamination by ultra-high vacuum (UHV) rest gases. While avoiding any current setpoint overshoot to preserve the tip integrity, we present results from approaches of atomically defined tungsten tips to the tunneling regime with Au(111), HOPG, and Si(111) surfaces at room temperature. We conclude from these experiments that adatom mobility and physisorbed gas on the sample surface limit the choice of surfaces for which the tip integrity is preserved in tunneling experiments at room temperature. The atomic structure of FIM tip apices is unchanged only after tunneling to the highly reactive Si(111) surface.
△ Less
Submitted 2 July, 2012;
originally announced July 2012.
-
Breakdown of the Kratky-Porod Wormlike Chain Model for Semiflexible Polymers in Two Dimensions
Authors:
Hsiao-Ping Hsu,
Wolfgang Paul,
Kurt Binder
Abstract:
By large-scale Monte Carlo simulations of semiflexible polymers in $d=2$ dimensions the applicability of the Kratky-Porod model is tested. This model is widely used as "standard model" for describing conformations and force versus extension curves of stiff polymers. It is shown that semiflexible polymers in $d=2$ show a crossover from hard rods to self-avoiding walks, the intermediate Gaussian reg…
▽ More
By large-scale Monte Carlo simulations of semiflexible polymers in $d=2$ dimensions the applicability of the Kratky-Porod model is tested. This model is widely used as "standard model" for describing conformations and force versus extension curves of stiff polymers. It is shown that semiflexible polymers in $d=2$ show a crossover from hard rods to self-avoiding walks, the intermediate Gaussian regime (implied by the Kratky-Porod model) is completely absent. Hence the latter can also describe force versus extension curves only if the contour length is only a few times larger than the persistence length. Consequences for experiments on biopolymers at interfaces are briefly discussed.
△ Less
Submitted 14 November, 2011;
originally announced November 2011.
-
Computer simulation of bottle brush polymers with flexible backbone: Good solvent versus Theta solvent conditions
Authors:
Panagiotis E. Theodorakis,
Hsiao-Ping Hsu,
Wolfgang Paul,
Kurt Binder
Abstract:
By Molecular Dynamics simulation of a coarse-grained bead-spring type model for a cylindrical molecular brush with a backbone chain of $N_b$ effective monomers to which with grafting density $σ$ side chains with $N$ effective monomers are tethered, several characteristic length scales are studied for variable solvent quality. Side chain lengths are in the range $5 \le N \le 40$, backbone chain len…
▽ More
By Molecular Dynamics simulation of a coarse-grained bead-spring type model for a cylindrical molecular brush with a backbone chain of $N_b$ effective monomers to which with grafting density $σ$ side chains with $N$ effective monomers are tethered, several characteristic length scales are studied for variable solvent quality. Side chain lengths are in the range $5 \le N \le 40$, backbone chain lengths are in the range $50 \le N_b \le 200$, and we perform a comparison to results for the bond fluctuation model on the simple cubic lattice (for which much longer chains are accessible, $N_b \le 1027$, and which corresponds to an athermal, very good, solvent). We obtain linear dimensions of side chains and the backbone chain and discuss their $N$-dependence in terms of power laws and the associated effective exponents. We show that even at the Theta point the side chains are considerably stretched, their linear dimension depending on the solvent quality only weakly. Effective persistence lengths are extracted both from the orientational correlations and from the backbone end-to-end distance; it is shown that different measures of the persistence length (which would all agree for Gaussian chains) are not mutually consistent with each other, and depend distinctly both on $N_b$ and the solvent quality. A brief discussion of pertinent experiments is given.
△ Less
Submitted 2 November, 2011;
originally announced November 2011.
-
A fast Monte Carlo algorithm for studying bottle-brush polymers
Authors:
Hsiao-Ping Hsu,
Wolfgang Paul
Abstract:
Obtaining reliable estimates of the statistical properties of complex macromolecules by computer simulation is a task that requires high computational effort as well as the development of highly efficient simulation algorithms. We present here an algorithm combining local moves, the pivot algorithm, and an adjustable simulation lattice box for simulating dilute systems of bottle-brush polymers wit…
▽ More
Obtaining reliable estimates of the statistical properties of complex macromolecules by computer simulation is a task that requires high computational effort as well as the development of highly efficient simulation algorithms. We present here an algorithm combining local moves, the pivot algorithm, and an adjustable simulation lattice box for simulating dilute systems of bottle-brush polymers with a flexible backbone and flexible side chains under good solvent conditions. Applying this algorithm to the bond fluctuation model, very precise estimates of the mean square end-to-end distances and gyration radii of the backbone and side chains are obtained, and the conformational properties of such a complex macromolecule are studied. Varying the backbone length (from $N_b=67$ to $N_b=1027$), side chain length (from N=0 to N=24 or 48), the scaling predictions for the backbone behavior as well as the side chain behavior are checked. We are also able to give a direct comparison of the structure factor between experimental data and the simulation results.
△ Less
Submitted 7 July, 2011;
originally announced July 2011.
-
Mesophase formation in two-component cylindrical bottle-brush polymers
Authors:
Igor Erukhimovich,
Panagiotis E. Theodorakis,
Wolfgang Paul,
Kurt Binder
Abstract:
When two types of side chains (A,B) are densely grafted to a (stiff) backbone and the resulting bottle-brush polymer is in a solution under poor solvent conditions, an incompatibility between A and B leads to microphase separation in the resulting cylindrical brush. The possible types of ordering are reminiscent of the ordering of block copolymers in cylindrical confinement. Starting from this ana…
▽ More
When two types of side chains (A,B) are densely grafted to a (stiff) backbone and the resulting bottle-brush polymer is in a solution under poor solvent conditions, an incompatibility between A and B leads to microphase separation in the resulting cylindrical brush. The possible types of ordering are reminiscent of the ordering of block copolymers in cylindrical confinement. Starting from this analogy, Leibler's theory of microphase separation in block copolymer melts is generalized to derive a description of the system in the weak segregation limit. Also molecular dynamics simulation results of a corresponding coarse-grained bead-spring model are presented. Using side chain lengths up to N = 50 effective monomers, the ratio of the Lennard-Jones energy parameter between unlike monomers $(ε_{AB})$ and monomers of the same kind $(ε_{AA} = ε_{BB})$ is varied. Various correlation functions are analyzed to study the conditions when (local) Janus cylinder-type ordering and when (local) microphase separation in the direction along the cylinder axis occurs. Both the analytical theory and the simulations give evidence for short range order due to a tendency towards microphase separation in the axial direction, with a wavelength proportional to the side chain gyration radius, irrespective of temperature and grafting density, for a wide range of these parameters.
△ Less
Submitted 10 May, 2011;
originally announced May 2011.
-
Polymer chain stiffness versus excluded volume: A Monte Carlo study of the crossover towards the wormlike chain model
Authors:
Hsiao-Ping Hsu,
Wolfgang Paul,
Kurt Binder
Abstract:
When the local intrinsic stiffness of a polymer chain varies over a wide range, one can observe both a crossover from rigid-rod-like behavior to (almost) Gaussian random coils and a further crossover towards self-avoiding walks in good solvents. Using the pruned-enriched Rosenbluth method (PERM) to study self-avoiding walks of up to $N_b=50000$ steps and variable flexibility, the applicability of…
▽ More
When the local intrinsic stiffness of a polymer chain varies over a wide range, one can observe both a crossover from rigid-rod-like behavior to (almost) Gaussian random coils and a further crossover towards self-avoiding walks in good solvents. Using the pruned-enriched Rosenbluth method (PERM) to study self-avoiding walks of up to $N_b=50000$ steps and variable flexibility, the applicability of the Kratky-Porod model is tested. Evidence for non-exponential decay of the bond-orientational correlations $<\cos θ(s) >$ for large distances $s$ along the chain contour is presented, irrespective of chain stiffness. For bottle-brush polymers on the other hand, where experimentally stiffness is varied via the length of side-chains, it is shown that these cylindrical brushes (with flexible backbones) are not described by the Kratky-Porod wormlike chain model, since their persistence length is (roughly) proportional to their cross-sectional radius, for all conditions of practical interest.
△ Less
Submitted 4 October, 2010;
originally announced October 2010.
-
Conformational Studies of bottle-brush polymers absorbed on a flat solid surface
Authors:
Hsiao-Ping Hsu,
Wolfgang Paul,
Kurt Binder
Abstract:
The adsorption of a bottle-brush polymer end-grafted with one chain end of its backbone to a flat substrate surface is studied by Monte Carlo simulation of a coarse-grained model, that previously has been characterized in the bulk, assuming a dilute solution under good solvent conditions. Applying the bond fluctuation model on the simple cubic lattice, we vary the backbone chain length $N_b$ from…
▽ More
The adsorption of a bottle-brush polymer end-grafted with one chain end of its backbone to a flat substrate surface is studied by Monte Carlo simulation of a coarse-grained model, that previously has been characterized in the bulk, assuming a dilute solution under good solvent conditions. Applying the bond fluctuation model on the simple cubic lattice, we vary the backbone chain length $N_b$ from $N_b=67$ to $N_b = 259$ effective monomeric units, the side chain length $N$ from N=6 to N=48, and the grafting density $σ=1$, i.e., parameters that correspond well to the experimentally accessible range. When the adsorption energy strength $ε$ is varied, we find that the adsorption transition (which becomes well-defined in the limit $N_b \rightarrow \infty$, for arbitrary finite $N$) roughly occurs at the same value $ε_c$ as for ordinary linear chains (N=0), at least within our statistical errors. Mean square end-to-end distances and gyration radii of the side chains are obtained, as well as the monomer density profile in the direction perpendicular to the adsorbing surface. We show that for longer side chains the adsorption of bottle-brushes is a two step process, the decrease of the perpendicular linear dimension of side chains with adsorption energy strength can even be non-monotonic. Also the behavior of the static structure factor $S(q)$ is analyzed, evidence for a quasi-two-dimensional scaling is presented, and consequences for the interpretation of experiments are discussed.
△ Less
Submitted 17 September, 2010;
originally announced September 2010.
-
A slow process in confined polymer melts: layer exchange dynamics at a polymer solid interface
Authors:
L. Yelash,
P. Virnau,
K. Binder,
W. Paul
Abstract:
Employing Molecular Dynamics simulations of a chemically realistic model of 1,4-polybutadiene between graphite walls we show that the mass exchange between layers close to the walls is a slow process already in the melt state. For the glass transition of confined polymers this process competes with the slowing down due to packing effects and intramolecular rotation barriers.
Employing Molecular Dynamics simulations of a chemically realistic model of 1,4-polybutadiene between graphite walls we show that the mass exchange between layers close to the walls is a slow process already in the melt state. For the glass transition of confined polymers this process competes with the slowing down due to packing effects and intramolecular rotation barriers.
△ Less
Submitted 16 September, 2010;
originally announced September 2010.
-
Two-state protein-like folding of a homopolymer chain
Authors:
Mark P. Taylor,
Wolfgang Paul,
Kurt Binder
Abstract:
Many small proteins fold via a first-order "all-or-none" transition directly from an expanded coil to a compact native state. Here we study an analogous direct freezing transition from an expanded coil to a compact crystallite for a simple flexible homopolymer. Wang-Landau sampling is used to construct the 1D density of states for square-well chains of length 128. Analysis within both the micro-ca…
▽ More
Many small proteins fold via a first-order "all-or-none" transition directly from an expanded coil to a compact native state. Here we study an analogous direct freezing transition from an expanded coil to a compact crystallite for a simple flexible homopolymer. Wang-Landau sampling is used to construct the 1D density of states for square-well chains of length 128. Analysis within both the micro-canonical and canonical ensembles shows that, for a chain with sufficiently short-range interactions, the usual polymer collapse transition is preempted by a direct freezing or "folding" transition. A 2D free-energy landscape, built via subsequent multi-canonical sampling, reveals a dominant folding pathway over a single free-energy barrier. This barrier separates a high entropy ensemble of unfolded states from a low entropy set of crystallite states and the transition proceeds via the formation of a transition-state folding nucleus. Despite the non-unique homopolymer ground state, the thermodynamics of this direct freezing transition are identical to the thermodynamics of two-state protein folding. The model chain satisfies the van't Hoff calorimetric criterion for two-state folding and an Arrhenius analysis of the folding/unfolding free energy barrier yields a Chevron plot characteristic of small proteins.
△ Less
Submitted 27 May, 2010;
originally announced May 2010.
-
How to Define Variation of Physical Properties Normal to an Undulating One-Dimensional Object
Authors:
Hsiao-Ping Hsu,
Kurt Binder,
Wolfgang Paul
Abstract:
One-dimensional flexible objects are abundant in physics, from polymers to vortex lines to defect lines and many more. These objects structure their environment and it is natural to assume that the influence these objects exert on their environment depends on the distance from the line-object. But how should this be defined? We argue here that there is an intrinsic length scale along the undulat…
▽ More
One-dimensional flexible objects are abundant in physics, from polymers to vortex lines to defect lines and many more. These objects structure their environment and it is natural to assume that the influence these objects exert on their environment depends on the distance from the line-object. But how should this be defined? We argue here that there is an intrinsic length scale along the undulating line that is a measure of its "stiffness" (i.e., orientational persistence), which yields a natural way of defining the variation of physical properties normal to the undulating line. We exemplify how this normal variation can be determined from a computer simulation for the case of a so-called bottle-brush polymer, where side chains are grafted onto a flexible backbone.
△ Less
Submitted 14 November, 2009;
originally announced November 2009.
-
Phase Transitions and Relaxation Processes in Macromolecular Systems: The Case of Bottle-brush Polymers
Authors:
Hsiao-Ping Hsu,
Wolfgang Paul,
Panagiotis E. Theodorakis,
Kurt Binder
Abstract:
As an example for the interplay of structure, dynamics, and phase behavior of macromolecular systems, this article focuses on the problem of bottle-brush polymers with either rigid or flexible backbones. On a polymer with chain length $N_b$, side-chains with chain length $N$ are endgrafted with grafting density $σ$. Due to the multitude of characteristic length scales and the size of these polym…
▽ More
As an example for the interplay of structure, dynamics, and phase behavior of macromolecular systems, this article focuses on the problem of bottle-brush polymers with either rigid or flexible backbones. On a polymer with chain length $N_b$, side-chains with chain length $N$ are endgrafted with grafting density $σ$. Due to the multitude of characteristic length scales and the size of these polymers (typically these cylindrical macromolecules contain of the order of 10000 effective monomeric units) understanding of the structure is a challenge for experiment. But due to excessively large relaxation times (particularly under poor solvent conditions) such macromolecules also are a challenge for simulation studies. Simulation strategies to deal with this challenge, both using Monte Carlo and Molecular Dynamics Methods, will be briefly discussed, and typical results will be used to illustrate the insight that can be gained.
△ Less
Submitted 22 October, 2009;
originally announced October 2009.
-
Coarse-grained models for fluids and their mixtures: Comparison of Monte Carlo studies of their phase behavior with perturbation theory and experiment
Authors:
B. M. Mognetti,
P. Virnau,
L. Yelash,
W. Paul,
K. Binder,
M. Mueller,
L. G. MacDowell
Abstract:
The prediction of the equation of state and the phase behavior of simple fluids (noble gases, carbon dioxide, benzene, methane, short alkane chains) and their mixtures by Monte Carlo computer simulation and analytic approximations based on thermodynamic perturbation theory is discussed. Molecules are described by coarse grained (CG) models, where either the whole molecule (carbon dioxide, benzen…
▽ More
The prediction of the equation of state and the phase behavior of simple fluids (noble gases, carbon dioxide, benzene, methane, short alkane chains) and their mixtures by Monte Carlo computer simulation and analytic approximations based on thermodynamic perturbation theory is discussed. Molecules are described by coarse grained (CG) models, where either the whole molecule (carbon dioxide, benzene, methane) or a group of a few successive CH_2 groups (in the case of alkanes) are lumped into an effective point particle. Interactions among these point particles are fitted by Lennard-Jones (LJ) potentials such that the vapor-liquid critical point of the fluid is reproduced in agreement with experiment; in the case of quadrupolar molecules a quadrupole-quadrupole interaction is included. These models are shown to provide a satisfactory description of the liquid-vapour phase diagram of these pure fluids. Investigations of mixtures, using the Lorentz-Berthelot (LB) combining rule, also produce satisfactory results if compared with experiment, while in some previous attempts (in which polar solvents were modelled without explicitly taking into account quadrupolar interaction), strong violations of the LB rules were required. For this reason, the present investigation is a step towards predictive modelling of polar mixtures at low computational cost. These very simple coarse-grained models of small molecules developed here should be useful e.g. for simulations of polymer solutions with such molecules as solvent.
△ Less
Submitted 1 December, 2008;
originally announced December 2008.
-
Structure of bottle-brush polymers in solution: A Monte Carlo test of models for the scattering function
Authors:
Hsiao-Ping Hsu,
Wolfgang Paul,
Kurt Binder
Abstract:
Extensive Monte Carlo results are presented for a lattice model of a bottle-brush polymer under good solvent or Theta solvent conditions. Varying the side chain length, backbone length, and the grafting density for a rigid straight backbone, both radial density profiles of monomers and side chain ends are obtained, as well as structure factors describing the scattering from a single side chain a…
▽ More
Extensive Monte Carlo results are presented for a lattice model of a bottle-brush polymer under good solvent or Theta solvent conditions. Varying the side chain length, backbone length, and the grafting density for a rigid straight backbone, both radial density profiles of monomers and side chain ends are obtained, as well as structure factors describing the scattering from a single side chain and from the total bottle-brush polymer. To describe the structure in the interior of a very long bottle-brush, a periodic boundary condition in the direction along the backbone is used, and to describe effects due to the finiteness of the backbone length, a second set of simulations with free ends of the backbone is performed. In the latter case, the inhomogeneity of the structure in the direction along the backbone is carefully investigated. We use these results to test various phenomenological models that have been proposed to interpret experimental scattering data for bottle-brush macromolecules. These models aim to extract information on the radial density profile of a bottle-brush from the total scattering via suitable convolution approximations. Possibilities to improve such models, guided by our simulation results, are discussed.
△ Less
Submitted 28 April, 2009; v1 submitted 11 August, 2008;
originally announced August 2008.
-
Spherically averaged versus angle-dependent interactions in quadrupolar fluids
Authors:
Bortolo Matteo Mognetti,
Martin Oettel,
Leonid Yelash,
Peter Virnau,
Wolfgang Paul,
Kurt Binder
Abstract:
Employing simplified models in computer simulation is on the one hand often enforced by computer time limitations but on the other hand it offers insights into the molecular properties determining a given physical phenomenon. We employ this strategy to the determination of the phase behaviour of quadrupolar fluids, where we study the influence of omitting angular degrees of freedom of molecules…
▽ More
Employing simplified models in computer simulation is on the one hand often enforced by computer time limitations but on the other hand it offers insights into the molecular properties determining a given physical phenomenon. We employ this strategy to the determination of the phase behaviour of quadrupolar fluids, where we study the influence of omitting angular degrees of freedom of molecules via an effective spherically symmetric potential obtained from a perturbative expansion. Comparing the liquid-vapor coexistence curve, vapor pressure at coexistence, interfacial tension between the coexisting phases, etc., as obtained from both the models with the full quadrupolar interactions and the (approximate) isotropic interactions, we find discrepancies in the critical region to be typically (such as in the case of carbon dioxide) of the order of 4%. However, when the Lennard-Jones parameters are rescaled such that critical temperatures and critical densities of both models coincide with the experimental results, almost perfect agreement between the above-mentioned properties of both models is obtained. This result justifies the use of isotropic quadrupolar potentials. We present also a detailed comparison of our simulations with a combined integral equation/density functional approach and show that the latter provides an accurate description except for the vicinity of the critical point.
△ Less
Submitted 25 March, 2008;
originally announced March 2008.
-
Efficient prediction of thermodynamic properties of quadrupolar fluids from simulation of a coarse-grained model: The case of carbon dioxide
Authors:
Bortolo Matteo Mognetti,
Leonid Yelash,
Peter Virnau,
Wolfgang Paul,
Kurt Binder,
Marcus Mueller,
Luis Gonzalez MacDowell
Abstract:
Monte Carlo simulations are presented for a coarse-grained model of real quadrupolar fluids. Molecules are represented by particles interacting with Lennard-Jones forces plus the thermally averaged quadrupole-quadrupole interaction. The properties discussed include the vapor-liquid coexistence curve, the vapor pressure along coexistence, and the surface tension. The full isotherms are also acces…
▽ More
Monte Carlo simulations are presented for a coarse-grained model of real quadrupolar fluids. Molecules are represented by particles interacting with Lennard-Jones forces plus the thermally averaged quadrupole-quadrupole interaction. The properties discussed include the vapor-liquid coexistence curve, the vapor pressure along coexistence, and the surface tension. The full isotherms are also accessible over a wide range of temperatures and densities. It is shown that the critical parameters (critical temperature, density, and pressure) depend almost linearly on a quadrupolar parameter $q=Q^{*4} /T^*$, $Q^*$ is the reduced quadrupole moment of the molecule and $T^*$ the reduced temperature.
The model can be applied to a variety of small quadrupolar molecules. We focus on carbon dioxide as a test case, but consider nitrogen and benzene, too. Experimental critical temperature, density and quadrupolar moment are sufficient to fix the parameters of the model. The resulting agreement with experiments is excellent and marks a significant improvement over approaches which neglect quadrupolar effects. The same coarse-grained model was also applied in the framework of Perturbation Theory (PT) in the Mean Spherical Approximation (MSA). As expected, the latter deviates from the Monte Carlo results in the critical region, but is reasonably accurate at lower temperatures.
△ Less
Submitted 8 January, 2008;
originally announced January 2008.
-
Transitions of tethered polymer chains: A simulation study with the bond fluctuation lattice model
Authors:
Jutta Luettmer-Strathmann,
Federica Rampf,
Wolfgang Paul,
Kurt Binder
Abstract:
A polymer chain tethered to a surface may be compact or extended, adsorbed or desorbed, depending on interactions with the surface and the surrounding solvent. This leads to a rich phase diagram with a variety of transitions. To investigate these transitions we have performed Monte Carlo simulations of a bond-fluctuation model with Wang-Landau and umbrella sampling algorithms in a two-dimensiona…
▽ More
A polymer chain tethered to a surface may be compact or extended, adsorbed or desorbed, depending on interactions with the surface and the surrounding solvent. This leads to a rich phase diagram with a variety of transitions. To investigate these transitions we have performed Monte Carlo simulations of a bond-fluctuation model with Wang-Landau and umbrella sampling algorithms in a two-dimensional state space. The simulations' density of states results have been evaluated for interaction parameters spanning the range from good to poor solvent conditions and from repulsive to strongly attractive surfaces. In this work, we describe the simulation method and present results for the overall phase behavior and for some of the transitions. For adsorption in good solvent, we compare with Metropolis Monte Carlo data for the same model and find good agreement between the results. For the collapse transition, which occurs when the solvent quality changes from good to poor, we consider two situations corresponding to three-dimensional (hard surface) and two-dimensional (very attractive surface) chain conformations, respectively. For the hard surface, we compare tethered chains with free chains and find very similar behavior for both types of chains. For the very attractive surface, we find the two-dimensional chain collapse to be a two-step transition with the same sequence of transitions that is observed for three-dimensional chains: a coil-globule transition that changes the overall chain size is followed by a local rearrangement of chain segments.
△ Less
Submitted 3 January, 2008;
originally announced January 2008.
-
Equation of State for Macromolecules of Variable Flexibility in Good Solvents: A Comparison of Techniques for Monte Carlo Simulations of Lattice Models
Authors:
V. A. Ivanov,
E. A. An,
L. A. Spirin,
M. R. Stukan,
M. Mueller,
W. Paul,
K. Binder
Abstract:
The osmotic equation of state for the athermal bond fluctuation model on the simple cubic lattice is obtained from extensive Monte Carlo simulations. For short macromolecules (chain length N=20) we study the influence of various choices for the chain stiffness on the equation of state. Three techniques are applied and compared in order to critically assess their efficiency and accuracy: the repu…
▽ More
The osmotic equation of state for the athermal bond fluctuation model on the simple cubic lattice is obtained from extensive Monte Carlo simulations. For short macromolecules (chain length N=20) we study the influence of various choices for the chain stiffness on the equation of state. Three techniques are applied and compared in order to critically assess their efficiency and accuracy: the repulsive wall method, the thermodynamic integration method (which rests on the feasibility of simulations in the grand canonical ensemble), and the recently advocated sedimentation equilibrium method, which records the density profile in an external (e.g. gravitation-like) field and infers, via a local density approximation, the equation of state from the hydrostatic equilibrium condition. We confirm the conclusion that the latter technique is far more efficient than the repulsive wall method, but we find that the thermodynamic integration method is similarly efficient as the sedimentation equilibrium method. For very stiff chains the onset of nematic order enforces the formation of isotropic-nematic interface in the sedimentation equilibrium method leading to strong rounding effects and deviations from the true equation of state in the transition regime.
△ Less
Submitted 23 July, 2007; v1 submitted 28 June, 2007;
originally announced June 2007.
-
One- and two-component bottle-brush polymers: simulations compared to theoretical predictions
Authors:
Hsiao-Ping Hsu,
Wolfgang Paul,
Kurt Binder
Abstract:
Scaling predictions and results from self-consistent field calculations for bottle-brush polymers with a rigid backbone and flexible side chains under good solvent conditions are summarized and their validity and applicability is assessed by a comparison with Monte Carlo simulations of a simple lattice model. It is shown that under typical conditions, as they are also present in experiments, onl…
▽ More
Scaling predictions and results from self-consistent field calculations for bottle-brush polymers with a rigid backbone and flexible side chains under good solvent conditions are summarized and their validity and applicability is assessed by a comparison with Monte Carlo simulations of a simple lattice model. It is shown that under typical conditions, as they are also present in experiments, only a rather weak stretching of the side chains is realized, and then the scaling predictions based on the extension of the Daoud-Cotton blob picture are not applicable.
Also two-component bottle brush polymers are considered, where two types (A,B) of side chains are grafted, assuming that monomers of different kind repel each other. In this case, variable solvent quality is allowed for, such that for poor solvent conditions rather dense cylinder-like structures result. Theories predict ``Janus Cylinder''-type phase separation along the backbone in this case. The Monte Carlo simulations, using the pruned-enriched Rosenbluth method (PERM) then are restricted to rather short side chain length. Nevertheless, evidence is obtained that the phase separation between an A-rich part of the cylindrical molecule and a B-rich part can only occur locally. The correlation length of this microphase separation can be controlled by the solvent quality. This lack of a phase transition is interpreted by an analogy with models for ferromagnets in one space dimension.
△ Less
Submitted 8 May, 2007;
originally announced May 2007.
-
Intramolecular phase separation of copolymer "bottle brushes": No sharp phase transition but a tunable length scale
Authors:
Hsiao-Ping Hsu,
Wolfgang Paul,
Kurt Binder
Abstract:
A lattice model for a symmetrical copolymer "bottle brush" molecule, where two types (A,B) of flexible side chains are grafted with one chain end to a rigid backbone, is studied by a variant of the pruned-enriched Rosenbluth method (PERM), allowing for simultaneous growth of all side chains in the Monte Carlo sampling. Choosing repulsive binary interactions between unlike monomers and varying th…
▽ More
A lattice model for a symmetrical copolymer "bottle brush" molecule, where two types (A,B) of flexible side chains are grafted with one chain end to a rigid backbone, is studied by a variant of the pruned-enriched Rosenbluth method (PERM), allowing for simultaneous growth of all side chains in the Monte Carlo sampling. Choosing repulsive binary interactions between unlike monomers and varying the solvent quality, it is found that phase separation into an $A$-rich part of the cylindrical molecule and a $B$-rich part can occur only locally. Long range order (in the direction of the backbone) does not occur, and hence the transition from the randomly mixed state of the bottle brush to the phase separated structure is strongly rounded, in contrast to corresponding mean field predictions. This lack of a phase transition can be understood from an analogy with spin models in one space dimension. We predict that the range of microphase separation along the bottle brush backbone can be controlled on the nanoscale by varying the solvent quality.
△ Less
Submitted 10 September, 2007; v1 submitted 7 July, 2006;
originally announced July 2006.
-
A global investigation of phase equilibria using the Perturbed-Chain Statistical-Associating-Fluid-Theory (PC-SAFT) approach
Authors:
L. Yelash,
M. Mueller,
W. Paul,
K. Binder
Abstract:
The recently developed Perturbed-Chain Statistical Associating Fluid Theory (PC-SAFT) is investigated for a wide range of model parameters including the parameter m representing the chain length and the thermodynamic temperature T and pressure p. This approach is based upon the first-order thermodynamic perturbation theory for chain molecules developed by Wertheim and Chapman et al. and includes…
▽ More
The recently developed Perturbed-Chain Statistical Associating Fluid Theory (PC-SAFT) is investigated for a wide range of model parameters including the parameter m representing the chain length and the thermodynamic temperature T and pressure p. This approach is based upon the first-order thermodynamic perturbation theory for chain molecules developed by Wertheim and Chapman et al. and includes dispersion interactions via the second-order perturbation theory of Barker and Henderson. We systematically study a hierarchy of models which are based on the PC-SAFT approach using analytical model calculations and Monte Carlo simulations. For one-component systems we find that the analytical model in contrast to the simulation results exhibits two phase-separation regions in addition to the common gas-liquid coexistence region: One phase separation occurs at high density and low temperature. The second demixing takes place at low density and high temperature where usually the ideal gas phase is expected in the phase diagram. These phenomena, which are referred to as "liquid-liquid" and "gas-gas" equilibria, give rise to multiple critical points in one-component systems, as well as to critical end points (CEP) and equilibria of three fluid phases, which can usually be found in multicomponent mixtures only. Furthermore, it is shown that the "liquid-liquid" demixing in this model is not a consequence of a "softened" repulsive interaction as assumed in the theoretical derivation of the model. Experimental data for the melt density of polybutadiene with molecular mass Mw=45000g/mol are correlated here using the PC-SAFT equation. It is shown that the discrepancies in modeling the polymer density at ambient temperature and high pressure can be traced back to ...
△ Less
Submitted 24 May, 2005;
originally announced May 2005.
-
Analysis of the Phase Transitions in Alkyl- Mica by Density and Pressure Profiles
Authors:
Hendrik Heinz,
Wolfgang Paul,
Ulrich W. Suter,
Kurt Binder
Abstract:
In a previous work [H. Heinz, H. J. Castelijns, U. W. Suter, J. Am. Chem. Soc. 115, 9500 (2003)], we developed an accurate force field and simulated the phase transitions in C18-mica (octadecyltrimethylammonium-mica) as well as the absence of such transitions in 2C18-mica (dioctadecyldimethylammonium-mica) between room temperature and 100 deg C. Here we analyze (i) average z coordinates of the c…
▽ More
In a previous work [H. Heinz, H. J. Castelijns, U. W. Suter, J. Am. Chem. Soc. 115, 9500 (2003)], we developed an accurate force field and simulated the phase transitions in C18-mica (octadecyltrimethylammonium-mica) as well as the absence of such transitions in 2C18-mica (dioctadecyldimethylammonium-mica) between room temperature and 100 deg C. Here we analyze (i) average z coordinates of the carbon atoms and interdigitation of the hydrocarbon bilayers, (ii) density profiles and (iii) pressure profiles of the structures along all Cartesian axes. In C18-mica, the standard deviation in the z coordinate for the chain atoms is high and more than doubles in the disordered phase. The order-disorder transition is accompanied by a change in the orientation of the ammonium headgroup, as well as decreasing tensile and shear stress in the disordered phase. In 2C18-mica, the standard deviation in the z coordinate for the chain atoms is low and does not increase remarkably on heating. The backbones display a highly regular structure, which is slightly obscured by rotations in the C18 backbones and minor headgroup displacements at 100 deg C. Close contacts between the bulky headgroups with sidearms cause significant local pressure which is in part not relieved at 100 deg C. An increase of the basal-plane spacing at higher temperature is found in both systems due to larger separation between the two hydrocarbon layers and an increased z spacing between adjacent chain atoms (=decreased tilt of the chains relative to the surface normal), and, in C18-mica only, a stronger upward orientation of the C18 chain at the ammonium headgroup. The likelihood for chain interdigitation between the two hydrocarbon layers is 0.24 to 0.30 for C18-mica,and 0.65 to 0.26 for 2C18-mica (for 20 deg C to 100 deg C).
△ Less
Submitted 23 November, 2003;
originally announced November 2003.
-
Calculation of local pressure tensors in systems with many-body interactions
Authors:
Hendrik Heinz,
Wolfgang Paul,
Kurt Binder
Abstract:
Local pressures are important in the calculation of interface tensions and in analyzing micromechanical behavior. The calculation of local pressures in computer simulations has been limited to systems with pairwise interactions between the particles, which is not sufficient for chemically detailed systems with many-body potentials such as angles and torsions. We introduce a method to calculate l…
▽ More
Local pressures are important in the calculation of interface tensions and in analyzing micromechanical behavior. The calculation of local pressures in computer simulations has been limited to systems with pairwise interactions between the particles, which is not sufficient for chemically detailed systems with many-body potentials such as angles and torsions. We introduce a method to calculate local pressures in systems with n-body interactions (n=2,3,4, . . .) based on a micromechanical definition of the pressure tensor. The local pressure consists of a kinetic contribution from the linear momentum of the particles and an internal contribution from dissected many-body interactions by infinitesimal areas. To define dissection by a small area, respective n-body interactions are divided into two geometric centers, effectively reducing them to two-body interactions. Consistency with hydrodynamics-derived formulas for systems with two-body interactions (J. H. Irving and J. G. Kirkwood, J. Chem. Phys. 18, 817 (1950)), for average cross-sectional pressures (B. D. Todd, D. J. Evans, and P. J. Daivis, Phys. Rev. E 52, 1627 (1995)), and for volume averaged pressures (virial formula) is shown. As a simple numerical example, we discuss liquid propane in a cubic box. Local, crosssectional,and volume-averaged pressures as well as relative contributions from two-body and three-body forces are analyzed with the proposed method, showing full numerical equivalence with the existing approaches. The method allows computing local pressures in the presence of many-body interactions in atomistic simulations of complex materials and biological systems.
△ Less
Submitted 24 February, 2006; v1 submitted 31 August, 2003;
originally announced September 2003.
-
Molecular Dynamics Simulations
Authors:
Kurt Binder,
Jurgen Horbach,
Walter Kob,
Wolfgang Paul,
Fathollah Varnik
Abstract:
A tutorial introduction to the technique of Molecular Dynamics (MD) is given, and some characteristic examples of applications are described. The purpose and scope of these simulations and the relation to other simulation methods is discussed, and the basic MD algorithms are described. The sampling of intensive variables (temperature T, pressure p) in runs carried out in the microcanonical (NVE)…
▽ More
A tutorial introduction to the technique of Molecular Dynamics (MD) is given, and some characteristic examples of applications are described. The purpose and scope of these simulations and the relation to other simulation methods is discussed, and the basic MD algorithms are described. The sampling of intensive variables (temperature T, pressure p) in runs carried out in the microcanonical (NVE) ensemble (N= particle number, V = volume, E = energy) is discussed, as well as the realization of other ensembles (e.g. the NVT ensemble). For a typical application example, molten SiO2, the estimation of various transport coefficients (self-diffusion constants, viscosity, thermal conductivity) is discussed. As an example of Non-Equilibrium Molecular Dynamics (NEMD), a study of a glass-forming polymer melt under shear is mentioned.
△ Less
Submitted 7 August, 2003;
originally announced August 2003.
-
Studies of YBCO Strip Lines under Voltage Pulses: Optimisation of the Design of Fault Current Limiters
Authors:
M. Decroux,
L. Antognazza,
S. Reymond,
W. Paul,
M. Chen,
O. Fischer
Abstract:
We present experimental results on the behaviour of a superconducting YBCO/Au meander of length L submitted to short circuit tests with constant voltage pulses. The meander, at the beginning of the short-circuit, is divided in two regions; one, with a length L1 proportional to the applied voltage, which first switches into a highly dissipative state (HDS) while the rest remains superconducting.…
▽ More
We present experimental results on the behaviour of a superconducting YBCO/Au meander of length L submitted to short circuit tests with constant voltage pulses. The meander, at the beginning of the short-circuit, is divided in two regions; one, with a length L1 proportional to the applied voltage, which first switches into a highly dissipative state (HDS) while the rest remains superconducting. Then the rest of the meander will progressively switch into the normal state due to the propagation of this HDS (few m/s) from both ends. The part L1 has to initially support a power density proportional to r.Jp^2 (r is the resistivity of the bilayer and Jp the peak current density). To avoid local excessive dissipation of power and over heating on one part of the wafer in the initial period, we have developed a novel design in order to distribute the dissipating section of the meander into many separated small dissipative zones. Furthermore the apparent propagation velocity of these dissipative zones is increased by the number of propagation fronts. We will show results obtained on 3kW (300V, 10A) FCL on a 2" wafer which confirm the benefits of this new design.
△ Less
Submitted 21 January, 2003;
originally announced January 2003.
-
Simulation of models for the glass transition: Is there progress?
Authors:
Kurt Binder,
Jörg Baschnagel,
Walter Kob,
Wolfgang Paul
Abstract:
The glass transition of supercooled fluids is a particular challenge for computer simulation, because the (longest) relaxation times increase by about 15 decades upon approaching the transition temperature T_g. Brute-force molecular dynamics simulations, as presented here for molten SiO_2 and coarse-grained bead-spring models of polymer chains, can yield very useful insight about the first few d…
▽ More
The glass transition of supercooled fluids is a particular challenge for computer simulation, because the (longest) relaxation times increase by about 15 decades upon approaching the transition temperature T_g. Brute-force molecular dynamics simulations, as presented here for molten SiO_2 and coarse-grained bead-spring models of polymer chains, can yield very useful insight about the first few decades of this slowing down. Hence this allows to access the temperature range around T_c of the so-called mode coupling theory, whereas the dynamics around the experimental glass transition is completely out of reach. While methods such as ``parallel tempering'' improve the situation somewhat, a method that allows to span a significant part of the region T_g\leq T\leq T_c is still lacking. Only for abstract models such as the infinite range 10-state Potts glass with a few hundred spins this region can be explored. However this model suffers from very strong finite size effects thus making it difficult to extrapolate the results obtained for the finite system sizes to the thermodynamic limit. For the case of polymer melts, two different strategies to use lattice models instead of continuum models are discussed: In the first approach, a mapping of an atomistically realistic model of polyethylene to the bond fluctuation model with suitable effective potentials and a temperature-dependent time rescaling factor is attempted. In the second approach, devoted to a test of the entropy theory, moves that are artificial but which lead to a faster relaxation (``slithering snake'' algorithm) are used, to get at least static properties at somewhat lower temperatures than possible with a ``realistic'' dynamics. The merits and shortcomings of all these approaches are discussed.
△ Less
Submitted 20 February, 2002;
originally announced February 2002.
-
Molecular-Dynamics Simulation of a Glassy Polymer Melt: Rouse Model and Cage Effect
Authors:
C. Bennemann,
J. Baschnagel,
W. Paul,
K. Binder
Abstract:
We report results of molecular-dynamics simulations for a glassy polymer melt consisting of short, linear bead-spring chains. It was shown in previous work that this onset of the glassy slowing down is compatible with the predictions of the mode coupling theory. The physical process of `caging' of a monomer by its spatial neighbors leads to a distinct two step behavior in the particle mean squar…
▽ More
We report results of molecular-dynamics simulations for a glassy polymer melt consisting of short, linear bead-spring chains. It was shown in previous work that this onset of the glassy slowing down is compatible with the predictions of the mode coupling theory. The physical process of `caging' of a monomer by its spatial neighbors leads to a distinct two step behavior in the particle mean square displacements. In this work we analyze the effects of this caging process on the Rouse description of the melt's dynamics. We show that the Rouse theory is applicable for length and time scales above the typical scales for the caging process. Futhermore, the monomer displacement is compared with simulation results for a binary Lennard-Jones mixture to point out the differences which are introduced by the connectivity of the particles.
△ Less
Submitted 26 February, 1999;
originally announced February 1999.
-
Investigating the influence of different thermodynamic paths on the structural relaxation in a glass forming polymer melt
Authors:
Christoph Bennemann,
Wolfgang Paul,
J"org Baschnagel,
Kurt Binder
Abstract:
We present results from Molecular Dynamics simulations of the thermal glass transition in a dense polymer melt. In previous work we compared the simulation data with the idealized version of mode coupling theory (MCT) and found that the theory provides a good description of the dynamics above the critical temperature. In order to investigate the influence of different thermodynamic paths on the…
▽ More
We present results from Molecular Dynamics simulations of the thermal glass transition in a dense polymer melt. In previous work we compared the simulation data with the idealized version of mode coupling theory (MCT) and found that the theory provides a good description of the dynamics above the critical temperature. In order to investigate the influence of different thermodynamic paths on the structural relaxation (alpha-process), we performed simulations for three different pressures and are thus able to give a sketch of the critical line of MCT in the pressure-temperature-plane [(p,T)-plane], where, according to the idealised version of MCT, an ergodic-nonergodic transition should occur. Furthermore, by cooling our system along two different paths (an isobar and an isochor), with the same impact point on the critical line, we demonstrate that neither the critical temperature nor the exponent gamma depend on the chosen path.
△ Less
Submitted 2 October, 1998;
originally announced October 1998.
-
Monte Carlo and Molecular Dynamics Simulation of the Glass Transition of Polymers
Authors:
K. Binder,
J. Baschnagel,
C. Bennemann,
W. Paul
Abstract:
Two coarse-grained models for polymer chains in dense glass-forming polymer melts are studied by computer simulation: the bond-fluctuation model on a simple cubic lattice, where a bond-length potential favors long bonds, is treated by dynamic Monte Carlo methods, and a bead-spring model in the continuum with a Lennard-Jones potential between the beads is treated by Molecular Dynamics. While the…
▽ More
Two coarse-grained models for polymer chains in dense glass-forming polymer melts are studied by computer simulation: the bond-fluctuation model on a simple cubic lattice, where a bond-length potential favors long bonds, is treated by dynamic Monte Carlo methods, and a bead-spring model in the continuum with a Lennard-Jones potential between the beads is treated by Molecular Dynamics. While the dynamics of both models differ for short length scales and associated time scales, on mesoscopic spatial and temporal scales both models behave similarly. In particular, the mode coupling theory of the glass transition can be used to interpret the slowing down of the undercooled polymer melt. For the off-lattice model, the approach to the critical point of mode coupling is both studied for constant pressure and for constant volume. The lattice model allows a test of the Gibbs-Di Marzio entropy theory of the glass transition, and our finding is that although the entropy does decrease significantly, there is no ``entropy catastrophe'', where the fluid entropy would turn negative. Finally, an outlook on confinement effects on the glass transition in thin film geometry is given.
△ Less
Submitted 29 September, 1998;
originally announced September 1998.
-
Molecular-Dynamics Simulation of a Glassy Polymer Melt: Incoherent Scattering Function
Authors:
C. Bennemann,
J. Baschnagel,
W. Paul
Abstract:
We present simulation results for a model polymer melt, consisting of short, nonentangled chains, in the supercooled state. The analysis focuses on the monomer dynamics, which is monitored by the incoherent intermediate scattering function. The scattering function is recorded over six decades in time and for many different wave-vectors. The lowest temperatures studied are slightly above the crit…
▽ More
We present simulation results for a model polymer melt, consisting of short, nonentangled chains, in the supercooled state. The analysis focuses on the monomer dynamics, which is monitored by the incoherent intermediate scattering function. The scattering function is recorded over six decades in time and for many different wave-vectors. The lowest temperatures studied are slightly above the critical temperature of mode-coupling theory (MCT), which was determined from a quantitative analysis of the beta- and alpha-relaxations. We find evidence for the space-time factorization theorem in the beta-relaxation regime, and for the time-temperature superposition principle in the alpha-regime, if the temperature is not too close to the critical temperature. The wave-vector dependence of the nonergodicity parameter, of the critical amplitude, and the alpha-relaxation time are in qualitative agreement with calculations for hard spheres. For wave-vectors larger than the maximum of the structure factor the alpha-relaxation time already agrees fairly well with the asymptotic MCT-prediction. The behavior of the relaxation time at small wave-vectors can be rationalized by the validity of the Gaussian approximation and the value of the Kohlrausch stretching exponent.
△ Less
Submitted 24 September, 1998;
originally announced September 1998.
-
Molecular Dynamics Simulations of the Thermal Glass Transition in Polymer Melts: Alpha-Relaxation Behaviour
Authors:
C. Bennemann,
W. Paul,
K. Binder,
B. Duenweg
Abstract:
We present Molecular Dynamics simulations of the thermal glass transition in a dense model polymer liquid. We performed a comparative study of both constant volume and constant pressure cooling of the polymer melt. Great emphasis was laid on a careful equilibration of the dense polymer melt at all studied temperatures. Our model introduces competing length scales in the interaction to prevent an…
▽ More
We present Molecular Dynamics simulations of the thermal glass transition in a dense model polymer liquid. We performed a comparative study of both constant volume and constant pressure cooling of the polymer melt. Great emphasis was laid on a careful equilibration of the dense polymer melt at all studied temperatures. Our model introduces competing length scales in the interaction to prevent any cristallisation tendency. In this first manuscript we we analyse the structural properties as a function of temperature and the long alpha-relaxation behaviour as observed in the dynamic structure factor and the self-diffusion of the polymer chains. The relaxation can be consistently analysed in terms of the mode coupling theory (MCT) of the glass transition. The mode coupling critical temperature and the exponent gamma, defining the power law divergence of the alpha-relaxation timescale both depend on the thermodynamic ensemble employed in the simulation.
△ Less
Submitted 14 October, 1997;
originally announced October 1997.
-
Orthorhombic Phase of Crystalline Polyethylene: A Constant Pressure Path Integral Monte Carlo Study
Authors:
R. Martonak,
W. Paul,
K. Binder
Abstract:
In this paper we present a Path Integral Monte Carlo (PIMC) simulation of the orthorhombic phase of crystalline polyethylene, using an explicit atom force field with unconstrained bond lengths and angles. This work represents a quantum extension of our recent classical simulation (J. Chem. Phys. 106, 8918 (1997)). It is aimed both at exploring the applicability of the PIMC method on such polymer…
▽ More
In this paper we present a Path Integral Monte Carlo (PIMC) simulation of the orthorhombic phase of crystalline polyethylene, using an explicit atom force field with unconstrained bond lengths and angles. This work represents a quantum extension of our recent classical simulation (J. Chem. Phys. 106, 8918 (1997)). It is aimed both at exploring the applicability of the PIMC method on such polymer crystal systems, as well as on a detailed assessment of the importance of quantum effects on different quantities. We used the $NpT$ ensemble and simulated the system at zero pressure in the temperature range 25 - 300 K, using Trotter numbers between 12 and 144. In order to investigate finite-size effects, we used chains of two different lengths, C_12 and C_24, corresponding to the total number of atoms in the super-cell being 432 and 864, respectively. We show here the results for structural parameters, like the orthorhombic lattice constants a,b,c, and also fluctuations of internal parameters of the chains, such as bond lengths and bond and torsional angles. We have also determined the internal energy and diagonal elastic constants c_11, c_22 and c_33. We discuss the temperature dependence of the measured quantities and compare to that obtained from the classical simulation. For some quantities, we discuss the way they are related to the torsional angle fluctuation. In case of the lattice parameters we compare our results to those obtained from other theoretical approaches as well as to some available experimental data. In order to study isotope effects, we simulated also a deuterated polyethylene crystal at a low temperature. We also suggest possible ways of extending this study and present some general considerations concerning modeling of polymer crystals.
△ Less
Submitted 18 August, 1997;
originally announced August 1997.
-
Small-Angle Excess Scattering: Glassy Freezing or Local Orientational Ordering?
Authors:
H. Weber,
W. Paul,
W. Kob,
K. Binder
Abstract:
We present Monte Carlo simulations of a dense polymer melt which shows glass-transition-like slowing-down upon cooling, as well as a build up of nematic order. At small wave vectors q this model system shows excess scattering similar to that recently reported for light-scattering experiments on some polymeric and molecular glass-forming liquids. For our model system we can provide clear evidence…
▽ More
We present Monte Carlo simulations of a dense polymer melt which shows glass-transition-like slowing-down upon cooling, as well as a build up of nematic order. At small wave vectors q this model system shows excess scattering similar to that recently reported for light-scattering experiments on some polymeric and molecular glass-forming liquids. For our model system we can provide clear evidence that this excess scattering is due to the onset of short-range nematic order and not directly related to the glass transition.
△ Less
Submitted 24 December, 1996;
originally announced December 1996.
-
Orthorhombic Phase of Crystalline Polyethylene: A Monte Carlo Study
Authors:
R. Martonak,
W. Paul,
K. Binder
Abstract:
In this paper we present a classical Monte Carlo simulation of the orthorhombic phase of crystalline polyethylene, using an explicit atom force field with unconstrained bond lengths and angles and periodic boundary conditions. We used a recently developed algorithm which apart from standard Metropolis local moves employs also global moves consisting of displacements of the center of mass of the…
▽ More
In this paper we present a classical Monte Carlo simulation of the orthorhombic phase of crystalline polyethylene, using an explicit atom force field with unconstrained bond lengths and angles and periodic boundary conditions. We used a recently developed algorithm which apart from standard Metropolis local moves employs also global moves consisting of displacements of the center of mass of the whole chains in all three spatial directions as well as rotations of the chains around an axis parallel to the crystallographic c-direction. Our simulations are performed in the NpT ensemble, at zero pressure, and extend over the whole range of temperatures in which the orthorhombic phase is experimentally known to be stable (10 - 450 K). In order to investigate the finite-size effects in this extremely anisotropic crystal, we used different system sizes and different chain lengths, ranging from C_12 to C_96 chains, the total number of atoms in the super-cell being between 432 and 3456. We show here the results for structural parameters, such as the orthorhombic cell parameters a,b,c, and the setting angle of the chains, as well as internal parameters of the chains, such as the bond lengths and angles. Among thermodynamic quantities, we present results for thermal expansion coefficients, elastic constants and specific heat. We discuss the temperature dependence of the measured quantities as well as the related finite-size effects. In case of lattice parameters and thermal expansion coefficients, we compare our results to those obtained from other theoretical approaches as well as to some available experimental data. We also suggest some possible ways of extending this study.
△ Less
Submitted 10 December, 1996;
originally announced December 1996.