Modeling of Electrochemical Cells: HYD Lecture 10. Atomistic and Coarse-Grained Models of Polymers
|
|
- Kimberly Hunt
- 6 years ago
- Views:
Transcription
1 Modeling of Electrochemical Cells: Proton Exchange Membrane Fuel Cells HYD Lecture 10. Atomistic and Coarse-Grained Models of Polymers Dept. of Chemical & iomolecular Engineering Yonsei University Spring, 2011 Prof. David Keffer Prof. David Keffer
2 Lecture Outline Atomistic simulations method sample results (Wang, Q., Keffer, D.J., Petrovan, S., Thomas, J.., Molecular Dynamics Simulation of Polyethylene Terephthalate Oligomers, J. Phys. Chem.. 114(2) 2010 pp ) Coarse-graining Procedure (Wang, Q., Keffer, D.J., Nicholson, D.M., Thomas, J.., Use of the Ornstein Zernike Percus Yevick Equation to Extract Interaction Potentials from Pair Correlation Functions, Phys. Rev. E 81(6) 2010 article # ) Coarse-grained simulations method sample results (Wang, Q., Keffer, D.J., Nicholson, D.M., Thomas, J.., Coarse grained Molecular Dynamics Simulation of Polyethylene Terephthalate (PET), Macromolecules 43(24) 2010 pp )
3 Motivation The simulations presented thus far in the course have focused on the structure and dynamics within the aqueous phase. The movement of water and ions occurs on a much faster timescale than the movement of the polymer. Thus, in those analyses, the polymer is considered to be essentially rigid. If we are interested in understanding processes that require polymer mobility, such as the response of the membrane to a sustained applied electric field, then we must implement simulations that are capable of much longer timescales. This is one motivation for moving to coarse-grained simulations of polymers.
4 Molecular l Level lsimulation i : Atomistic MD simulation: Small system, short time. number of oligomers 125 to 216 number of atoms 4320 to simulation duration 1 to 100 ns system size ~ 6 nm The use ofperiodic boundary conditions allows bulk properties to be extracted from these simulations.
5 Can t Simulate long chains hi with ihatomistic i Model Why? Longest relaxational time of polymers scales with chain length to the 3 4 power. PET Relaxation Times DP time (ns) For long chains (DP > 10), the relaxation time is longer than the maximum simulation time available on modern computational resources What is the solution? Coarse grained models
6 Molecular l vs. Coarse Grained dmodels: Example Polymer: PET Molecular Model Coarse Grained Model A atomistic detail particles are molecular l groups structure is known grouping is based on intuition interaction potentials known interaction potentials unknown many degrees of freedom computationally expensive limited to short chains fewer degrees of freedom computationally modest works for long chains use for DP 1 to 10 use for DP 25 to 50
7 Procedure 1. Run molecular level simulation for DP = 1 to Parameterize coarse grained potential based on molecularlevel simulations 3. Run coarse grained simulations for DP = 25 to 50
8 Molecular l Dynamics Simulation: i ased on solving Newton s equation of motion To simulate N atoms in 3 D, you must solve a set of 3N coupled nonlinear ordinary differential equations. F ma The force is completely determined by an interaction potential. F U The ODE for particle i in dimension is thus d x i, 1 U dt m x 2 i, We must provide an interaction potential from either theory, quantum mechanical calculations or experiment.
9 Molecular l Dynamics Simulation: i Intramolecular l potential il We use the anisotropic united atom (AUA) potential model developed by Hedenqvist, haradwaj and oyd for PET to describe the intramolecular and intermolecular potential of HET molecules. Namely, all atoms of HET molecule are explicitly represented expect for the hydrogen connected to carbons θ 5 4 1: bond stretching 2: bond bending 3: bond torsion 4: out of plane bending 5: intramolecular L J and electrostatic potential (for atoms over four bonds) Note: not all pairs are listed O: red; C: grey; H: white Hedenqvist et al. Macromolecules (neglected OH interactions) OH interactions from in Chen et al. J. Phys. Chem., 2001
10 Collaboration with Oak Ridge National Laboratory To solve the thousands of ODEs (largest system thus far is 55,000), we use the massively parallel supercomputers at ORNL, both Cray and IM. National Center for Computational Science Today the computing resources of the NCCS are among the fastest in the world, able to perform more than 119 trillion calculations per second. These resources are available to researchers at UT through discretionary accounts of the program directors.
11 Molecular Dynamic Simulation Snapshots all molecules shown selected molecules shown Snapshots of PET hexamer at T = 563 K and p = 0.13 kpa.
12 Molecular Dynamic Simulation Snapshots dimer trimer tetramer hexamer octamer decamer All snaphots at T = 563 K and p = 0.13 kpa.
13 Structural properties: Relaxation Times monomer Two traditional models: both fit well <R(t).R(0)>/<R R 2 > MD simulation fit to KWW model fit to Rouse model 0.4 exp(-1/τ R ) β exp(-(t/τ KWW )) time (ns) Tsolou, V.G. et al. Macromolecules, 2005 MD simulations provide characteristic ti relaxations times through h the analysis of the decay of structural auto correlation functions. For polymers, longest relaxation time is due to rotation.
14 Structural properties: Relaxation Times DP time (ns) end to end dista ance auto correla ation function monomer dimer trimer tetramer hexamer octamer time (ns) Characteristic relaxations times increase dramatically with DP, creating difficulties for MD simulation of longer chain molecules.
15 Structural properties: End to end Pair Correlation Function end to end pair corre elation function function end to end pair correlation f monomer dimer trimer DP <R > (Å) ete tetramer hexamer octamer trimer tetramer hexamer octamer distance (Angstroms) ± ± ± ± ± ± distance (Angstroms) The long tail of end to end pair correlation function is typical for polymer molecules, indicating the complexity of chain molecules.
16 Structural properties: Hydrogen onding Hydrogen onding as a function of DP monomer dimer trimer tetramer hexamer octamer g OH (r) distance (Angstrom) While the number of alcohol end groups decreases with DP While the number of alcohol end groups decreases with DP, the probability of finding the end groups forming hydrogen bonds increase with DP.
17 Thermodynamic properties: Density Melt density as a function of DP DP < ρ> > (g/cm 3 ) ±0.006 density (g/cm^3 3) monomer PET melt ± ± ± Need experimental references. Is bulk the right word? Report different state points ± ± ±0.034 DP The density does not change much as DPincreases, the melt density of PET is around 1.2 (g/cm 3 ) (ASTM D1238,T=285 ), our simulation results are comparable to the experimental data.
18 Thermodynamic properties: Enthalpy Enthalpy per molecule as a function of DP H (aj/molecu ule) enthalpy per molecule enthalpy per DP DP <H> DP (kj /mole) ± ± ± ± ± ± ± The enthalpy per molecule increases with DP. The enthalpy per repeat unit also increases with DP.
19 Thermodynamic properties: Heatcapacity (C p ), isothermal compressibility thermal expansion coefficient Fluctuations method Finite difference method: (one simulation): (multiple simulations) 1 1 V p+ε,t -V p-ε,t 2 V Vk T 1 (V)(H+PV) 2 Vk T C P (H+PV) Vk T β= - V(p,T) 2ε 1 V p,t+ε -V p,t-ε α= V(p,T) 2ε C= p H p,t+ε -H p,t-ε 2ε Give a small perturbation of density and temperature. Calculate changes of pressure, volume, and enthalpy. Derivative is the corresponding thermodynamic property. The fluctuation method for these thermodynamic properties is based on the variance of volume, enthalpy and covariance of volume and enthalpy of the system, which can be obtained from one MD simulation.
20 Thermodynamic properties: Heatcapacity (C p ), isothermal compressibility(β), thermal expansion coefficient (α) Finite difference : very close to the experimental data Fluctuation: ti too much uncertainty. Note: experimental ldt data is tk taken from Polymer Handbook written by J. randrup, E.H. Immergut, and E.A. Grulke 1999.
21 Transport properties: Diffusivity Diffusion coefficient is calculated from the center of mass mean square displacement DP D (m 2 /sec))x10 11 log (<(R (t)-r cm (0))2 >( 2 cm )) nomomer(slope:1.00) dimer (slope: 1.00) trimer (slope:1.01) tetramer (slope:1.00) hexamer (slope:1.01) octamer (slope: 1.01) ± ± ± ± 4.01 ± ± ± ±0.60 log (simulation time (fs)) The slope of the curve equal to 1 means our simulations are long enough to guarantee a good diffusivity. The self diffusion coefficient decreases with DP as the diffusion process becomes harder for longer molecules.
22 Transport properties: Viscosity Zero shear rate viscosity is obtained by integrating momentum autocorrelation function over time, by following Green Kubo relations. 5e-9 5e-9 n function tum auto correlatio moment 4e-9 3e-9 2e-9 1e-9 momentum auto correlation functio on 4e-9 3e-9 2e-9 1e-9 0-1e-9-2e-9 time (ps) monomer dimer trimer tetramer hexamer octamer -3e The decay to zero ensures convergence of the integration. 0-1e time (ps)
23 Transport properties:viscosity Zero shear viscosity as a function of DP. viscosity (pa.s) monomer dimer trimer tetramer hexamer octamer tetramer trimer octamer hexamer DP η(pa.s) x ± ± 5.45 ± ± ± 6.51 ± monomer dimer time (ps) ± ± ±8.0 The viscosity increases with DP, this is true based on polymer physics. We measured the viscosity of the monomer using a rheometer at UT. The experimental viscosity of the monomer at T=533 K and 1 atm is 3.8 x 10 3 Pa.s. Compare to a simulated viscosity at T=563 K and 1 atm of 3.83 x10 3 Pa.s.
24 Transport properties: Stokes Einstein equation Relationship between diffusion coefficient and viscosity R kt 4 D D kt 1 4 R ln (Rg) & ln (R) radius of gyration Rg(MD) R(MD) fit (R(MD)) fit (Rg(MD)) Once we know R, either D or η as a function of DP, we can calculate the other by using above equation ln (DP) For monomer and dimer, R agrees with R g, R increases with DP, fluctuating around its trend line with slope equals to 1. This is governed by scaling behavior of D and η.
25 Transport properties: Thermal conductivity ased on integration of heat flux auto correlation function over time. 2e-4 heat flux au uto correlation fu unction 2e-4 1e-4 5e-55 0 monomer dimer trimer tetramer hexemer octamer elation function heat flux auto corr 2.0e-4 1.5e-4 1.0e-4 5.0e e-5-1.0e-4-1.5e time (ps) -5e time (ps) The heat flux auto correlation function has a faster decay compared with the momentum auto correlation function.
26 Transport Properties: Thermal Conductivity thermal con nductivity (W/m m.k*10^-7) 4e-8 3e-8 2e-8 1e-8 monomer dimer trimer tetramer hexamer octamer octamer trimer monomer tetramer hexamer dimer DP λ(w/m.k) ± ± ± ± ± ± time ( ps) ±0.01 The thermal conductivity obtained from simulation is between 0.1 and 0.3(W/m.K), consistent with the experimental data (0.24 from polymer handbook).
27 Scaling Analysis Objectives: The structural, thermodynamic and transport properties all scale with DP according to polymer dynamic theory, these scaling behaviors enable us to verify our simulation results and give predictions to the properties of longer chains, such as DP=10,which isdifficultto to obtain frommdsimulation simulation. Principles: Study how the physical quantities change when DP is changed. Properties such as R g, R ete, kww, C p, D are all a(dp) b, The scaling exponent b is the constant we are interested in and different from property to property D a DP b Method: y fitting our simulation data to a(dp) b in log log plots, the slope of the curve is b. The literature reported values of b for polymer melts are as follows: b = 2.1 for D (reptation model) b = for R g, R ete b =1 for (below l critical molecular l weight) iht) M. Doi and S.F. Edwards, The Theory of Polymer Dynamics, G. Tsolou, V.G. Mavrantzas, D.N. Theodourou, Macromolecules, 2005.
28 ScalingAnalysis: Chain radius of gyration (R g ) and average end to end distance (R ete ) 5 ) and ln (R g ) 4 3 R g (fit), slope=0.594 R ete (fit),slope=0.571 R g (MD) R ete (MD) Estimated value for decamer: <R g > =15.19 (Å) <R ete > =35.65 (Å) ln (R ete ) ln (DP) The scaling exponent for R g and R ete is and separately, very close to theoretical prediction of Our result is valid and we also give the prediction of the value of R g and R ete for decamer.
29 ScalingAnalysis: Heat Capacity (C p ) MD fit (slope=-0.081, R2=0.99) b=0.081 from MD ln ( Cp) 0.9 no reported literature value Estimated value for decamer: C p =2250 (J K 1 kg 1 ) ln (DP)
30 ScalingAnalysis: Diffusion Coefficient (D) MD fit (slope=-2.01, R 2 =0.99) MD simulation b = 2.01 ln (D) theoretical value: b = ln (DP) Estimated value for decamer: D =6.40x ( m 2 /sec)
31 ScalingAnalysis: Zero Shear Rate Viscosity (η) MD fit ( slope=0.96, R 2 =0.85) MD value b = 0.96 ln (viscosity) Literature reported -5 value is 1 (for PET with molecular weight -6 under critical point) ln (DP) Estimated value for decamer: η =2.24x10 2 (pa.s)
32 Coarse Grained model: based on atomisticmd simulation: Atomistic MD simulation is limited to DP up to 10, CG MD simulation will allow us to get to DP from 25 to 50 or even longer chains by neglecting gsome local degrees of freedom. A A A A A Atomistic model and com of CG beads CG model A CG model on top of atomistic model for hexamer (DP=6) CG chain for hexamer, total interaction centers is 13 for each molecule ( original:90).
33 Coarse Grained Potentials: parameterization We must figure out all interaction modes according to the CG model and parameterize the CG potentials for each mode (stretching, bending, torsion, nonbonded interactions) to match the structures of atomistic simulation. A AA bend A A A stretch CG CG CG CG CG U U U U U AA torsion A bend CG chain model and intramolecular interaction for hexamer, 1: bond stretching (A); 2: bond bending (A and AA); 3: bond torsion (AA) A A A stretch bend torsion non bond CG CG U (, r T) k Tln P (, r T) C stretch stretch r CG CG U (, T) k Tln P (, T) C bend bend CG CG Utorsion (, T ) k T ln Ptorsion (, T) C To get the potential of each mode, we have to obtain the distribution function first, this can be extracted by analyzing atomistic simulation configurations of short chains.
34 Coarse Grained potentials: distributions of CG beads We obtained P(r), P( ), P( ), g(r) based on atomistic MD simulation configurations. P(r) A stretching tetramer hexamer octamer decamer tetramer hexamer octamer decamer P( A-tetramer AA-tetramer A-hexamer AA-hexamer A-octamer AA-octamer A-decamer AA-decamer bending AA A separation r (Angostroms) (Angstroms) bending angle (degree AA tetramer hexamer octamer decamer non bonded P( ) torsion g(r) tetramer hexamer octamer decamer torsion angle degree separation (Angstroms)
35 Coarse Grained potentials: OZPY 1 method for CG polymers OZ integral equation: exact relationship between pair correlation function (PCF) and interaction potential g r 1 c ( r ) n c s h t d V A A A 1 A A 3 2 A OZPY: given U, find PCF OZPY 1 : given PCF, find U need to know: PCFs allowable combinations (A) n linear chain various combinations of stretching, bending and torsion interactions For this model, there are 34 allowable combinations that contribute to the indirect portion of the correlation.
36 Coarse Grained potentials: bonded and non bonded.
37 Coarse Grained MD simulation: comparison of structures. We compared the P(r), P( ) and P( ) from CG and atomistic MD simulation. istance distribution stretching di tion to orsion angle disbribu A A target CGMD separation (Angstroms) AA-target AA-CGMD AA torsion angle (degrees) istribution bending angle d A-target A-CGMD AA-target AA-CGMD AA bending angle (degrees) A The distributions of beads with bonded interactions from MD simulations match relatively well.
38 Coarse Grained MDsimulation: comparison ofstructures comparison of non bonded PCFs (, A and AA) from CGMD and atomistic MD simulations nonbonded distance disyribution target -CGMD nonbonded distance distribution AA AA-target AA-CGMD distance (Angstroms) distance (Angstroms) ribution nonb bonded distance disbr A distance (Angstroms) A-target A-CGMD The non bonded PCFs from CGMD and atomistic i MD simulations i match fil fairly well. Solving the integral equations provides us a way to get a reliable non bonded interaction potential for CGMD simulation.
39 Coarse Grained MD simulation: chain end to end distance distribution bution pr robability distri DP = 10 DP = DP = DP = 40 DP = Decamer: Comparison of chain end to end distance distributions from CGMD and atomistic MD match fairly well distance (Angstroms) The distributions are not typical Gaussian distributions, with chain length increases, they become more Gaussian like.
40 CGMD simulation: comparison of dynamic property chain end to end distance auto correlation function (ACF) Note: We must scale time in CGMD simulation (Harmandaris et al. Macromol. Chem. Phys. (2007)) (based on ratio of the values of selfdiffusivities) Atomistic MD: ACFs for long chain systmems can not reach 0 in short times For Atomistic MD: ACFs for long chain systmems can not reach 0 in short times. For decamer, it took roughly 6 months to finish a run of 30ns CGMD: simulations of the same systems. Apparent speed up is observed. All ACFs can reach 0 in short times. For decamer, it took just 2 weeks to finish a run of 180 ns.
41 Scalingfactor and Scalingexponent (b) Scaling factors: based on the ratio of the values of self diffusivities from atomistic and CG simulations. The value 5.38 is used to scale the dynamic properties back to atomistic scope. Scaling exponents: b from polymer physics. DP 1~10 X a Simulation method Atomistic MD DP b D KWW <R ete > <R g > ~10 CGMD ~50 CGMD Rouse Model (1, 2) Reptation Model (1, 2) N/A N/A Tzoumanekas et al. Macromolecules Lahmar et al. Macromolecules 2009
42 Coarse Grained MD simulation: Scaling exponents of R ete and R g 6 ln (R ete ) and ln (R g ) for Long chains: obtained b is close to 0.5 for short chains: obtained b is in the range of 0.57~ ln (DP) The values for both short chains hi and long chains hi agree with the Rouse and Reptation theory respectively.
43 Coarse Grained MD simulation: Scaling exponent of self diffusivity i it Rouse model: for short chains, b = 1 ln(d(m 2 /s)) Reptation model: for long chains, b = ln(dp) Obtained b deviatesfrom Rousetheoryfor shortchains, agrees with reptation theory for long chains.
44 Coarse Grained dmd simulation: entanglement tanalysis To further understand the reptation behavior, we can do entanglement analysisby extracting entanglement information directly from configurations of the chains. The Z code: A common algorithm to study the entanglements in polymeric systems (Kroger, M. Comput. Phys. Commun ) What can we get from Z code? mean contour length of primitive path (<Lpp>), tube diameter (d), number of monomers between entanglement points (N e ), interentanglement strand length (N ES ), defined as: N ES Z N ( ( N N 1) 1) N Kamio et al. Macromolecules 2007, 40, 710.
45 Coarse Grained MD simulation: entanglement analysis DP <L pp > (Å) d (Å) N e Z N ES rheology models N/A 35 a, 30.2a, b 24.2 b, N/A 14.9 d 25.0 c For DP = 10, unentangled system. for DP = 20 to 50, (d), (Ne) and (N ES ) are very close to the reported values for entangled PET melts. a Fetters et al InPhysical Properties ofpolymers Handbook; James E Mark 2007 a Fetters et al. In Physical Properties of Polymers Handbook; James E. Mark, 2007 b Fetters et al. Macromolecules 1994, 27, c Lorentz, G.; Tassin, J. F. Polymer 1994, 35, d Kamio et al. Macromolecules 2007, 40, 710.
46 Conclusions: Atomistic simulations of PET from DP = 1 to 10 yield a variety of structural, thermodynamic and transport properties that compare well with experimental data. CG model of PET was implemented in CGMD simulations of PET chains with DP up to 50. This allows simulation up to 1000 ns. The CG potential is parameterized to structural distribution functions obtained from atomistic simulations using OZPY 1. The CGMD simulation of PET chains satisfactorily reproduces the structural and dynamic properties from atomistic MD simulation of the same systems. For the long gchains, we find the escaling exponents epoetsof 0.51, 05,050a 0.50 and 2.00 for average chain end to end distance, radius of gyration and self diffusivity respectively. Th t l t l i h th t t b di t (d) (N ) d (N ) The entanglement analysis shows that tube diameter (d), (Ne) and (N ES ) of long chain systems are very close to the reported values for entangled PET melts.
Multiscale Materials Modeling
Multiscale Materials Modeling Lecture 06 Polymers These notes created by David Keffer, University of Tennessee, Knoxville, 2012. Outline Multiscale Modeling of Polymers I. Introduction II. Atomistic Simulation
More informationMultiscale Materials Modeling
Multiscale Materials Modeling Lecture 02 Capabilities of Classical Molecular Simulation These notes created by David Keffer, University of Tennessee, Knoxville, 2009. Outline Capabilities of Classical
More informationPart III. Polymer Dynamics molecular models
Part III. Polymer Dynamics molecular models I. Unentangled polymer dynamics I.1 Diffusion of a small colloidal particle I.2 Diffusion of an unentangled polymer chain II. Entangled polymer dynamics II.1.
More informationMultiscale Materials Modeling At the University of Tennessee
Multiscale Materials Modeling At the University of Tennessee David Keffer Dept. of Materials Science & Engineering The University of Tennessee Knoxville, TN 37996-2100 dkeffer@utk.edu http://clausius.engr.utk.edu/
More informationModeling of Electrochemical Cells: HYD Lecture 06. Multiscale Modeling of Water & Charge Transport
Modeling of Electrochemical Cells: Proton Exchange Membrane Fuel Cells HYD7007 01 Lecture 06. Multiscale Modeling of Water & Charge Transport Dept. of Chemical & Biomolecular Engineering Yonsei University
More informationGlass-Transition and Side-Chain Dynamics in Thin Films: Explaining. Dissimilar Free Surface Effects for Polystyrene and Poly(methyl methacrylate)
Supporting Information for Glass-Transition and Side-Chain Dynamics in Thin Films: Explaining Dissimilar Free Surface Effects for Polystyrene and Poly(methyl methacrylate) David D. Hsu, Wenjie Xia, Jake
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 1 David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Time/s Multi-Scale Modeling Based on SDSC Blue Horizon (SP3) 1.728 Tflops
More informationSupporting Information
Projection of atomistic simulation data for the dynamics of entangled polymers onto the tube theory: Calculation of the segment survival probability function and comparison with modern tube models Pavlos
More informationPolymer dynamics. Course M6 Lecture 5 26/1/2004 (JAE) 5.1 Introduction. Diffusion of polymers in melts and dilute solution.
Course M6 Lecture 5 6//004 Polymer dynamics Diffusion of polymers in melts and dilute solution Dr James Elliott 5. Introduction So far, we have considered the static configurations and morphologies of
More informationUsing Molecular Dynamics to Compute Properties CHEM 430
Using Molecular Dynamics to Compute Properties CHEM 43 Heat Capacity and Energy Fluctuations Running an MD Simulation Equilibration Phase Before data-collection and results can be analyzed the system
More informationStructural Properties from Molecular Dynamics Simulation
Structural Properties from Molecular Dynamics Simulation David Keffer Department of Materials Science & Engineering University of Tennessee, Knoxville date begun: April 7, 2016 date last updated: April
More informationCoarse-Grained Models!
Coarse-Grained Models! Large and complex molecules (e.g. long polymers) can not be simulated on the all-atom level! Requires coarse-graining of the model! Coarse-grained models are usually also particles
More informationMultiscale Modeling of Materials
Multiscale Modeling of Materials David Keffer Dept. of Materials Science & Engineering The University of Tennessee Knoxville, TN 37996-2100 dkeffer@utk.edu http://clausius.engr.utk.edu/ ASM Summer Materials
More informationOn the Dynamics and Disentanglement in Thin and Two-Dimensional Polymer Films
J. Phys. IV France 1 (006) Pr1-1 c EDP Sciences, Les Ulis On the Dynamics and Disentanglement in Thin and Two-Dimensional Polymer Films H. Meyer, T. Kreer, A. Cavallo, J. P. Wittmer and J. Baschnagel 1
More informationChapter 7. Entanglements
Chapter 7. Entanglements The upturn in zero shear rate viscosity versus molecular weight that is prominent on a log-log plot is attributed to the onset of entanglements between chains since it usually
More informationViscoelastic Flows in Abrupt Contraction-Expansions
Viscoelastic Flows in Abrupt Contraction-Expansions I. Fluid Rheology extension. In this note (I of IV) we summarize the rheological properties of the test fluid in shear and The viscoelastic fluid consists
More informationModeling of Electrochemical Cells: HYD Lecture 04. Overview of transport processes in PEMs
Modeling of Electrochemical Cells: Proton Exchange Membrane Fuel Cells HYD7007 01 Lecture 04. Overview of transport processes in PEMs Dept. of Chemical & Biomolecular Engineering Yonsei University Spring,
More informationChap. 2. Polymers Introduction. - Polymers: synthetic materials <--> natural materials
Chap. 2. Polymers 2.1. Introduction - Polymers: synthetic materials natural materials no gas phase, not simple liquid (much more viscous), not perfectly crystalline, etc 2.3. Polymer Chain Conformation
More informationQENS in the Energy Domain: Backscattering and Time-of
QENS in the Energy Domain: Backscattering and Time-of of-flight Alexei Sokolov Department of Polymer Science, The University of Akron Outline Soft Matter and Neutron Spectroscopy Using elastic scattering
More informationPolymer Dynamics. Tom McLeish. (see Adv. Phys., 51, , (2002)) Durham University, UK
Polymer Dynamics Tom McLeish Durham University, UK (see Adv. Phys., 51, 1379-1527, (2002)) Boulder Summer School 2012: Polymers in Soft and Biological Matter Schedule Coarse-grained polymer physics Experimental
More informationChapter 5: Molecular Scale Models for Macroscopic Dynamic Response. Fluctuation-Dissipation Theorem:
G. R. Strobl, Chapter 6 "The Physics of Polymers, 2'nd Ed." Springer, NY, (1997). R. B. Bird, R. C. Armstrong, O. Hassager, "Dynamics of Polymeric Liquids", Vol. 2, John Wiley and Sons (1977). M. Doi,
More informationA MOLECULAR DYNAMICS STUDY OF POLYMER/GRAPHENE NANOCOMPOSITES
A MOLECULAR DYNAMICS STUDY OF POLYMER/GRAPHENE NANOCOMPOSITES Anastassia N. Rissanou b,c*, Vagelis Harmandaris a,b,c* a Department of Applied Mathematics, University of Crete, GR-79, Heraklion, Crete,
More information(Polymer rheology Analyzer with Sliplink. Tatsuya Shoji JCII, Doi Project
Rheology Simulator PASTA (Polymer rheology Analyzer with Sliplink model of entanglement) Tatsuya Shoji JCII, Doi Project 0 sec -3 msec -6 sec -9 nsec -12 psec -15 fsec GOURMET SUSHI PASTA COGNAC MUFFIN
More informationSupplementary Information for Atomistic Simulation of Spinodal Phase Separation Preceding Polymer Crystallization
Supplementary Information for Atomistic Simulation of Spinodal Phase Separation Preceding Polymer Crystallization Richard H. Gee * Naida Lacevic and Laurence E. Fried University of California Lawrence
More informationEntanglements. M < M e. M > M e. Rouse. Zero-shear viscosity vs. M (note change of slope) Edwards degennes Doi. Berry + Fox, slope 3.4.
Entanglements Zero-shear viscosity vs. M (note change of slope) M < M e Rouse slope 3.4 M > M e Edwards degennes Doi slope 1 Berry + Fox, 1968 Question: Which factors affect the Me: T, P, M, flexibility,
More informationAnalytical Rheology Linear Viscoelasticity of Model and Commercial Long-Chain-Branched Polymer Melts
Analytical Rheology Linear Viscoelasticity of Model and Commercial Long-Chain-Branched Polymer Melts Sachin Shanbhag, Seung Joon Park, Qiang Zhou and Ronald G. Larson Chemical Engineering, University of
More informationNotes. Prediction of the Linear Viscoelastic Shear Modulus of an Entangled Polybutadiene Melt from Simulation and Theory (1) 3π 2 k B T D(T)N (2)
134 Macromolecules 2001, 34, 134-139 Notes Prediction of the Linear Viscoelastic Shear Modulus of an Entangled Polybutadiene Melt from Simulation and Theory Oleksiy Byutner and Grant D. Smith* Department
More informationMP10: Process Modelling
MP10: Process Modelling MPhil Materials Modelling Dr James Elliott 0.1 MP10 overview 6 lectures on process modelling of metals and polymers First three lectures by JAE Introduction to polymer rheology
More informationPart III. Polymer Dynamics molecular models
Part III. Polymer Dynamics molecular models I. Unentangled polymer dynamics I.1 Diffusion of a small colloidal particle I.2 Diffusion of an unentangled polymer chain II. Entangled polymer dynamics II.1.
More information1.3 Molecular Level Presentation
1.3.1 Introduction A molecule is the smallest chemical unit of a substance that is capable of stable, independent existence. Not all substances are composed of molecules. Some substances are composed of
More informationChapter 6: The Rouse Model. The Bead (friction factor) and Spring (Gaussian entropy) Molecular Model:
G. R. Strobl, Chapter 6 "The Physics of Polymers, 2'nd Ed." Springer, NY, (1997). R. B. Bird, R. C. Armstrong, O. Hassager, "Dynamics of Polymeric Liquids", Vol. 2, John Wiley and Sons (1977). M. Doi,
More information2.1 Traditional and modern applications of polymers. Soft and light materials good heat and electrical insulators
. Polymers.1. Traditional and modern applications.. From chemistry to statistical description.3. Polymer solutions and polymer blends.4. Amorphous polymers.5. The glass transition.6. Crystalline polymers.7.
More informationMolecular-Level Description of Reaction Equilibrium
Molecular-Level Description of Reaction Equilibrium CBE 450: Reactor Engineering Fundamentals Fall, 2009 Prof. David Keffer dkeffer@utk.edu Thermodynamic Properties Quantum Mechanics can be combined with
More informationPolymer Dynamics and Rheology
Polymer Dynamics and Rheology 1 Polymer Dynamics and Rheology Brownian motion Harmonic Oscillator Damped harmonic oscillator Elastic dumbbell model Boltzmann superposition principle Rubber elasticity and
More informationCoarse-grained Models for Oligomer-grafted Silica Nanoparticles
Modeling So+ Ma-er: Linking Mul3ple Length and Time Scales KITP Conference, Santa Barbara, June 4-8, 2012 Coarse-grained Models for Oligomer-grafted Silica Nanoparticles Bingbing Hong Alexandros Chremos
More informationMolecular Theories of Linear Viscoelasticity THE ROUSE MODEL (P. 1)
THE ROUSE ODEL (P. 1) odel polymer dynamics by a system of N + 1 beads connected by N springs. Figure 1: apping the Polymer Chain onto a Chain of Beads Connected by Springs. ROUSE SCALING Recall that a
More informationViscosity and Diffusion Constants Calculation of n-alkanes by Molecular Dynamics Simulations
Viscosity and Diffusion Constants Calculation of n-alkanes by Molecular Dynamics Simulations Song Hi Lee * and Taihyun Chang Department of Chemistry, Kyungsung University, Busan 608-736, Korea Department
More informationRouse chains, unentangled. entangled. Low Molecular Weight (M < M e ) chains shown moving past one another.
Physical Picture for Diffusion of Polymers Low Molecular Weight (M < M e ) chains shown moving past one another. Z X Y Rouse chains, unentangled Figure by MIT OCW. High Molecular weight (M > M e ) Entanglements
More informationHyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea
The Scaled Effective Solvent Method for Predicting the Equilibrium Ensemble of Structures with Analysis of Thermodynamic Properties of Amorphous Polyethylene Glycol-Water Mixtures Hyeyoung Shin a, Tod
More informationAppendix C - Persistence length 183. Consider an ideal chain with N segments each of length a, such that the contour length L c is
Appendix C - Persistence length 183 APPENDIX C - PERSISTENCE LENGTH Consider an ideal chain with N segments each of length a, such that the contour length L c is L c = Na. (C.1) If the orientation of each
More informationRHEOLOGY OF BRANCHED POLYMERS
RHEOLOGY OF BRANCHED POLYMERS Overview: The Tube Model Shear and elongational viscosity Albena Lederer Leibniz-Institute of Polymer Research Dresden Member of Gottfried Wilhelm Leibniz Society WGL Hohe
More informationSupplementary Information for: Controlling Cellular Uptake of Nanoparticles with ph-sensitive Polymers
Supplementary Information for: Controlling Cellular Uptake of Nanoparticles with ph-sensitive Polymers Hong-ming Ding 1 & Yu-qiang Ma 1,2, 1 National Laboratory of Solid State Microstructures and Department
More informationChapter 1 Introduction
Chapter 1 Introduction This thesis is concerned with the behaviour of polymers in flow. Both polymers in solutions and polymer melts will be discussed. The field of research that studies the flow behaviour
More informationA Phenomenological Model for Linear Viscoelasticity of Monodisperse Linear Polymers
Macromolecular Research, Vol. 10, No. 5, pp 266-272 (2002) A Phenomenological Model for Linear Viscoelasticity of Monodisperse Linear Polymers Kwang Soo Cho*, Woo Sik Kim, Dong-ho Lee, Lee Soon Park, Kyung
More informationMultiscale Materials Modeling
Multiscale Materials Modeling Lecture 09 Quantum Mechanics/Molecular Mechanics (QM/MM) Techniques Fundamentals of Sustainable Technology These notes created by David Keffer, University of Tennessee, Knoxville,
More informationUniversty of Patras, Patras, GR 26504, Greece
8 th GRACM International Congress on Computational Mechanics Volos, 12 July 15 July 2015 STRUCTURAL, CONFORMATIONAL, DYNAMIC AND TOPOLOGICAL PROPERTIES OF RING POLY(ETHYLE OXIDE) MELTS FROM MOLECULAR DYNAMICS
More informationMechanical Properties of Tetra-Polyethylene and Tetra-Polyethylene Oxide Diamond Networks via Molecular Dynamics Simulations
Supplemental Information Mechanical Properties of Tetra-Polyethylene and Tetra-Polyethylene Oxide Diamond Networks via Molecular Dynamics Simulations Endian Wang and Fernando A. Escobedo Table S1 Lennard-Jones
More informationHydrodynamics, Thermodynamics, and Mathematics
Hydrodynamics, Thermodynamics, and Mathematics Hans Christian Öttinger Department of Mat..., ETH Zürich, Switzerland Thermodynamic admissibility and mathematical well-posedness 1. structure of equations
More informationMultiscale Materials Modeling At the University of Tennessee
Multiscale Materials Modeling At the University of Tennessee David Keffer Dept. of Materials Science & Engineering The University of Tennessee Knoxville, TN 37996-2100 dkeffer@utk.edu http://clausius.engr.utk.edu/
More informationNAPLES New Algorism for Polymeric Liquids Entangled and Strained. Yuichi Masubuchi ICR, Kyoto University
NAPLES New Algorism for Polymeric Liquids Entangled and Strained Yuichi Masubuchi ICR, Kyoto University The ul:mate goal of polymer rheology Predict viscoelas:c behaviors from Mw Mw/Mn LCB Blends Copolymers
More informationEnergy Barriers and Rates - Transition State Theory for Physicists
Energy Barriers and Rates - Transition State Theory for Physicists Daniel C. Elton October 12, 2013 Useful relations 1 cal = 4.184 J 1 kcal mole 1 = 0.0434 ev per particle 1 kj mole 1 = 0.0104 ev per particle
More informationPolymers Dynamics by Dielectric Spectroscopy
Polymers Dynamics by Dielectric Spectroscopy What s a polymer bulk? A condensed matter system where the structural units are macromolecules Polymers Shape of a Macromolecule in the Bulk Flory's prediction
More informationSystematic Coarse-Graining and Concurrent Multiresolution Simulation of Molecular Liquids
Systematic Coarse-Graining and Concurrent Multiresolution Simulation of Molecular Liquids Cameron F. Abrams Department of Chemical and Biological Engineering Drexel University Philadelphia, PA USA 9 June
More informationMolecular Dynamics Simulation of Methanol-Water Mixture
Molecular Dynamics Simulation of Methanol-Water Mixture Palazzo Mancini, Mara Cantoni University of Urbino Carlo Bo Abstract In this study some properties of the methanol-water mixture such as diffusivity,
More informationModeling of Electrochemical Cells: HYD Lecture 08. Composite Membranes
Modeling of Electrochemical Cells: Proton Exchange Membrane Fuel Cells HYD7007 01 Lecture 08. Composite Membranes Dept. of Chemical & Biomolecular Engineering Yonsei University Spring, 2011 Prof. David
More informationThermodynamic behaviour of mixtures containing CO 2. A molecular simulation study
Thermodynamic behaviour of mixtures containing. A molecular simulation study V. Lachet, C. Nieto-Draghi, B. Creton (IFPEN) Å. Ervik, G. Skaugen, Ø. Wilhelmsen, M. Hammer (SINTEF) Introduction quality issues
More informationKirkwood-Buff Integrals for Aqueous Urea Solutions Based upon the Quantum Chemical Electrostatic Potential and Interaction Energies
Supporting Information for Kirkwood-Buff Integrals for Aqueous Urea Solutions Based upon the Quantum Chemical Electrostatic Potential and Interaction Energies Shuntaro Chiba, 1* Tadaomi Furuta, 2 and Seishi
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationDynamic force matching: Construction of dynamic coarse-grained models with realistic short time dynamics and accurate long time dynamics
for resubmission Dynamic force matching: Construction of dynamic coarse-grained models with realistic short time dynamics and accurate long time dynamics Aram Davtyan, 1 Gregory A. Voth, 1 2, a) and Hans
More informationNon equilibrium thermodynamics: foundations, scope, and extension to the meso scale. Miguel Rubi
Non equilibrium thermodynamics: foundations, scope, and extension to the meso scale Miguel Rubi References S.R. de Groot and P. Mazur, Non equilibrium Thermodynamics, Dover, New York, 1984 J.M. Vilar and
More informationVagelis A. Harmandaris, Kostas Ch. Daoulas, and Vlasis G. Mavrantzas*
5796 Macromolecules 2005, 38, 5796-5809 Molecular Dynamics Simulation of a Polymer Melt/Solid Interface: Local Dynamics and Chain Mobility in a Thin Film of Polyethylene Melt Adsorbed on Graphite Vagelis
More informationV. A. Harmandaris,, V. G. Mavrantzas,*, D. N. Theodorou,, M. Kro1 ger,, J. Ramírez, H. C. O2 ttinger, and D. Vlassopoulos
1376 Macromolecules 2003, 36, 1376-1387 Crossover from the Rouse to the Entangled Polymer Melt Regime: Signals from Long, Detailed Atomistic Molecular Dynamics Simulations, Supported by Rheological Experiments
More informationMultiscale Materials Modeling
Multiscale Materials Modeling Lecture 10 Adaptive Resolution Fundamentals of Sustainable Technology These notes created by David Keffer, University of Tennessee, Knoxville, 2012. Reminder In the previous
More informationBioengineering 215. An Introduction to Molecular Dynamics for Biomolecules
Bioengineering 215 An Introduction to Molecular Dynamics for Biomolecules David Parker May 18, 2007 ntroduction A principal tool to study biological molecules is molecular dynamics simulations (MD). MD
More informationSupporting Information for. Dynamics of Architecturally Engineered All- Polymer Nanocomposites
Supporting Information for Dynamics of Architecturally Engineered All- Polymer Nanocomposites Erkan Senses,,,,* Madhusudan Tyagi,, Madeleine Pasco, Antonio Faraone,* NIST Center for Neutron Research, National
More informationK ex. Conformational equilibrium. equilibrium K B
Effects of Chemical Exchange on NMR Spectra Chemical exchange refers to any yprocess in which a nucleus exchanges between two or more environments in which its NMR parameters (e.g. chemical shift, scalar
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy
More informationModeling the combined effect of surface roughness and shear rate on slip flow of simple fluids
Modeling the combined effect of surface roughness and shear rate on slip flow of simple fluids Anoosheh Niavarani and Nikolai Priezjev www.egr.msu.edu/~niavaran November 2009 A. Niavarani and N.V. Priezjev,
More informationShear Thinning Near the Rough Boundary in a Viscoelastic Flow
Advanced Studies in Theoretical Physics Vol. 10, 2016, no. 8, 351-359 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/astp.2016.6624 Shear Thinning Near the Rough Boundary in a Viscoelastic Flow
More informationGEM4 Summer School OpenCourseWare
GEM4 Summer School OpenCourseWare http://gem4.educommons.net/ http://www.gem4.org/ Lecture: Polymer Chains by Ju Li. Given August 16, 2006 during the GEM4 session at MIT in Cambridge, MA. Please use the
More informationThe Effect of Model Internal Flexibility Upon NEMD Simulations of Viscosity
Draft: September 29, 1999 The Effect of Model Internal Flexibility Upon NEMD Simulations of Viscosity N. G. Fuller 1 and R. L. Rowley 1,2 Abstract The influence of model flexibility upon simulated viscosity
More informationOrigins of Mechanical and Rheological Properties of Polymer Nanocomposites. Venkat Ganesan
Department of Chemical Engineering University of Texas@Austin Origins of Mechanical and Rheological Properties of Polymer Nanocomposites Venkat Ganesan $$$: NSF DMR, Welch Foundation Megha Surve, Victor
More informationHierarchical Modeling of Polystyrene: From Atomistic to Coarse-Grained Simulations
6708 Macromolecules 2006, 39, 6708-6719 Hierarchical Modeling of Polystyrene: From Atomistic to Coarse-Grained Simulations V. A. Harmandaris, N. P. Adhikari, N. F. A. van der Vegt, and K. Kremer* Max Planck
More informationStatistical and Thermal Physics. Problem Set 5
Statistical and Thermal Physics xford hysics Second year physics course Dr A. A. Schekochihin and Prof. A. T. Boothroyd (with thanks to Prof. S. J. Blundell Problem Set 5 Some useful constants Boltzmann
More informationAccuracy, Transferability, and Efficiency of Coarse-Grained Models of Molecular Liquids
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes. Cite This: pubs.acs.org/jpcb
More informationScientific Computing II
Scientific Computing II Molecular Dynamics Simulation Michael Bader SCCS Summer Term 2015 Molecular Dynamics Simulation, Summer Term 2015 1 Continuum Mechanics for Fluid Mechanics? Molecular Dynamics the
More informationMolecular Dynamics Simulation Study of Transport Properties of Diatomic Gases
MD Simulation of Diatomic Gases Bull. Korean Chem. Soc. 14, Vol. 35, No. 1 357 http://dx.doi.org/1.51/bkcs.14.35.1.357 Molecular Dynamics Simulation Study of Transport Properties of Diatomic Gases Song
More informationMichael W. Mahoney Department of Physics, Yale University, New Haven, Connecticut 06520
JOURNAL OF CHEMICAL PHYSICS VOLUME 115, NUMBER 23 15 DECEMBER 2001 Quantum, intramolecular flexibility, and polarizability effects on the reproduction of the density anomaly of liquid water by simple potential
More informationCase study: molecular dynamics of solvent diffusion in polymers
Course MP3 Lecture 11 29/11/2006 Case study: molecular dynamics of solvent diffusion in polymers A real-life research example to illustrate the use of molecular dynamics Dr James Elliott 11.1 Research
More informationFile ISM02. Dynamics of Soft Matter
File ISM02 Dynamics of Soft Matter 1 Modes of dynamics Quantum Dynamics t: fs-ps, x: 0.1 nm (seldom of importance for soft matter) Molecular Dynamics t: ps µs, x: 1 10 nm Brownian Dynamics t: ns>ps, x:
More informationEXPERIMENT 3. HEAT-CAPACITY RATIOS FOR GASES
EXERIMENT 3. HEAT-CAACITY RATIOS FOR GASES The ratio Cp/Cv of the heat capacity of a gas at constant pressure to that at constant volume will be determined by either the method of adiabatic expansion.
More informationMutual diffusion in the ternary mixture of water + methanol + ethanol: Experiments and Molecular Simulation
- 1 - Mutual diffusion in the ternary mixture of water + methanol + ethanol: Experiments and Molecular Simulation Tatjana Janzen, Gabriela Guevara-Carrión, Jadran Vrabec University of Paderborn, Germany
More informationSimulation of Coarse-Grained Equilibrium Polymers
Simulation of Coarse-Grained Equilibrium Polymers J. P. Wittmer, Institut Charles Sadron, CNRS, Strasbourg, France Collaboration with: M.E. Cates (Edinburgh), P. van der Schoot (Eindhoven) A. Milchev (Sofia),
More informationsimulations Extracting coarse-grained elastic behavior of capsid proteins from molecular dynamics Outline Feb. 2009
Extracting coarse-grained elastic behavior of capsid proteins from molecular dynamics simulations Stephen D. Hicks Christopher L. Henley Cornell University Gordon Conf. on Physical Virology Feb. 2009 Outline
More informationChemical Engineering 160/260 Polymer Science and Engineering. Lecture 14: Amorphous State February 14, 2001
Chemical Engineering 160/260 Polymer Science and Engineering Lecture 14: Amorphous State February 14, 2001 Objectives! To provide guidance toward understanding why an amorphous polymer glass may be considered
More information3.1 Hydrodynamic interactions in a Gaussian chain
Chapter 3 The Zimm model 3. Hydrodynamic interactions in a Gaussian chain In the previous chapter we have focused on the Rouse chain, which gives a good description of the dynamics of unentangled concentrated
More informationLecture 5: Macromolecules, polymers and DNA
1, polymers and DNA Introduction In this lecture, we focus on a subfield of soft matter: macromolecules and more particularly on polymers. As for the previous chapter about surfactants and electro kinetics,
More informationAnalysis of the simulation
Analysis of the simulation Marcus Elstner and Tomáš Kubař January 7, 2014 Thermodynamic properties time averages of thermodynamic quantites correspond to ensemble averages (ergodic theorem) some quantities
More informationUNIVERSITY OF KWAZULU-NATAL WESTVILLE CAMPUS DEGREE/DIPLOMA EXAMINATIONS: NOVEMBER 2006 CHEMISTRY CHEM230W: PHYSICAL CHEMISTRY 2
UNIVERSITY OF KWAZULU-NATAL WESTVILLE CAMPUS DEGREE/DIPLOMA EXAMINATIONS: NOVEMBER 006 CHEMISTRY CHEM30W: PHYSICAL CHEMISTRY TIME: 180 MINUTES MARKS: 100 EXAMINER: PROF S.B. JONNALAGADDA ANSWER FIVE QUESTIONS.
More informationElements of Polymer Structure and Viscoelasticity. David M. Parks Mechanics and Materials II February 18, 2004
Elements of Polymer Structure and Viscoelasticity David M. Parks Mechanics and Materials II 2.002 February 18, 2004 Outline Elements of polymer structure Linear vs. branched; Vinyl polymers and substitutions
More informationThis is a repository copy of Molecular observation of contour-length fluctuations limiting topological confinement in polymer melts.
This is a repository copy of Molecular observation of contour-length fluctuations limiting topological confinement in polymer melts. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/57/
More informationMathematical and Computational Modelling of Molecular Systems
Mathematical and Computational Modelling of Molecular Systems Vagelis Harmandaris Institute of Applied and Computational Mathematics (IACM/FORTH), Heraklion, Greece 11 η Επιστημονική Διημερίδα ΙΤΕ MS:
More informationPart 8. Special Topic: Light Scattering
Part 8. Special Topic: Light Scattering Light scattering occurs when polarizable particles in a sample are placed in the oscillating electric field of a beam of light. The varying field induces oscillating
More informationENAS 606 : Polymer Physics
ENAS 606 : Polymer Physics Professor Description Course Topics TA Prerequisite Class Office Hours Chinedum Osuji 302 Mason Lab, 432-4357, chinedum.osuji@yale.edu This course covers the static and dynamic
More informationCitation for published version (APA): Bulacu, M. I. (2008). Molecular dynamics studies of entangled polymer chains s.n.
University of Groningen Molecular dynamics studies of entangled polymer chains Bulacu, Monica Iulia IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite
More informationEQUATION OF STATE DEVELOPMENT
EQUATION OF STATE DEVELOPMENT I. Nieuwoudt* & M du Rand Institute for Thermal Separation Technology, Department of Chemical Engineering, University of Stellenbosch, Private bag X1, Matieland, 760, South
More informationComputational Modeling of Protein-Ligand Interactions
Computational Modeling of Protein-Ligand Interactions Steven R. Gwaltney Department of Chemistry Mississippi State University Mississippi State, MS 39762 Auguste Comte, 1830 Every attempt to refer chemical
More informationEffect of equilibration on primitive path analyses of entangled polymers
PHYSICAL REVIEW E 72, 061802 2005 Effect of equilibration on primitive path analyses of entangled polymers Robert S. Hoy* and Mark O. Robbins Department of Physics and Astronomy, Johns Hopkins University,
More informationSubject of the Lecture:
Subject of the Lecture: Conceptual basis for the development of force fields. Implementation/validation Water - a worked example Extensions - combining molecular mechanics and quantum mechanics (QM/MM)
More informationElectronic Supplementary Information
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2014 Electronic Supplementary Information Computational investigation of structural
More information