multiscale molecular simulations of pv encapsulants · 2018-08-08 · multiscale molecular...

1
Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-AC04-94AL85000. SAND No. 2011-XXXXP Multiscale Molecular Simulations of PV Encapsulants Lauren Abbott 1 , Ross Larsen 2 , and Stan Moore 1 1 Sandia National Laboratories, Albuquerque, NM, 2 National Renewable Energy Laboratory, Golden, CO Automated generationof candidate molecules for materials discovery Automated generationofinitial conditions for molecular simulations Automated extractionofsub-regions for analysis and calculationofproperties Integration with analysis and visualizationalgorithms Object-oriented framework for easy extension of the functionality to specificneeds Simulation Toolkit for Renewable Energy and Advanced Materials Modeling streamm.nrel.gov Open source classical molecular dynamics code Particle simulations at the atomic, meso, and continuum scales Inclusions of common potentials for hard and soft materials Efficient parallel simulations usingspatial decomposition of domain space Modular framework for easy extension with newfeatures and functionality Large-scale Atomic/Molecular Massively Parallel Simulator lammps.sandia.gov Polymer/Clay Nanocomposites Sandia and Texas A&M are workingon inexpensive,transparentpolymer/clay nanocomposites with layers of polymers and oriented clay platelets Tailoring the clay particles and polymer matrix controls the barrier properties, composite integrity, and fire retardancy Chemically-specificcoarse-grainedmodels of polymer/claynanocomposites can yield insight into material processes (e.g., encapsulation,intercalation,and exfoliation) and resultingmaterialproperties Clay platelets, ~50-100 nm diam *Sandia: Margaret Gordon, Erik Spoerke, Eric Schindelholtz, Ken Armijo, Rob Sorensen, Texas A&M: Jamie Grunlan, Kevin Holder, Shuang Qin neat polymer composite Stress-strain curve Left: Coarse-grained model; Middle: Example of polymer intercalation; Right: Clay dispersion in polymer matrix with polymer not shown *JL Suter, D Groen, PV Coveney, Adv. Mater., 2015, 27, 966-984; Nano Lett., 2015, 15, 8108-8113 Atomistic simulations elucidated possible effects of morphology on electron-transport mechanisms in radical polymer electrodes Molecular models and MD simulations of PTMA were setup using STREAMM Two primary distances (4.5 and 6.5 Å) were identified as contributing to an effective electron transfer distance of 5.5 Å Inter-site couplings were computed between >10,000 pairs of sites with QM calculations, showingthat electron transfer mostly takes place between sites on different chains Radical Polymer Electrodes Top: N-N radial distribution functions; Bottom: Average inter-site electronic coupling MD QM PTMA STREAMM’s automated structure generation from monomer to oligomer to bulk material, and automated extraction of sub- regions for QM calculations *TW Kemper, RE Larsen, T Gennett, J. Phys. Chem. C, 2014, 118, 17213-17220 Polymer Photovoltaic Active Layers Atomistic MD simulations were setup using STREAMM and performed in LAMMPS for three organicPV copolymers (BDT-TPD,PTB7, and PTB7-Th), with variations to the backbone and side-chain structure Aligned parallel chains with pi stacking were observed in BDT-TPD and PTB7, but not in PTB7-Th due to steric hindrance from the side-chain,which contradicts assumption made in the literature Transport in BDT-TPD and PTB7 likely occurs between parallel pi stacks, while transportin PTB7-Th likely occurs between orthogonal pi stacks BDT-TPD PTB7 PTB7-Th Molecular structures (top) and simulation snapshots (bottom) for organic PV copolymers *NREL: Travis Kemper and Ross Larsen Ion-Containing Polymers Sulfonated polyphenylenes(SDAPP) show promisefor proton exchange membrane fuel cells (PEMFCs) and vanadium flow batteries Little is known aboutthe nanoscalestructureof these amorphous polymers, which is not easily characterized with currentexperimental techniques Atomistic simulations yielded unique insight into the morphology of the ionic domains formed by the aggregation of water and ionic groups,specifically the size, shape, and connectivity, as well as resulting implications for ion transport As more water is added into the system, the ionic domains become more fully percolated and the domains become slightly larger and more spherical in shape, which would improveion transport increasing water Left: Molecular structure of SDAPP; Right: Ionic domains with disparate clusters shown in different colors SDAPP *Sandia: Lauren Abbott, Amalie Frischknecht, Cy Fujimoto, Todd Alam, Eric Sorte Molecular simulations can provideimportantinformationaboutnanoscale structure and phenomena in PV encapsulantmaterials with atomic-level detail not currently accessible to experimental techniques The atomic-level detail of molecular simulations are usefulfor studyinglocal packing of the polymer, as well as the interactions and dynamics of additives (e.g., UV absorbers)and contaminants (e.g., H 2 O or O 2 ) within the polymer Molecular simulations are a cost-effective route for efficient high-throughput screening of potentialencapsulantmaterials A multiscaleapproach combines quantummechanical(QM) calculations and atomistic and coarse-grainedmolecular dynamics (MD) simulations to capture a broad range of length and time scales Results from the molecular simulations can be passed to simulation models at higher lengths scales, such as meso or continuum techniques This capability leverages open source tools like Sandia’s LAMMPS MD package and NREL’s STREAMM toolkit to setup, run, and analyze molecular simulations Overview MD techniques follow classical dynamics using Newton’s law (F = ma), by updatingthe positions of particles using the net forces on the particles Potentials are described using force fields with terms for bonded interactions (e.g., bonds and angles) and nonbondedinteractions (e.g., van der Waals) Molecular Modeling QM techniques approximate the wave function followingSchrödinger’s equation ( = ψ/∂t ) to consider quantum effects We can move between different scales using fine- grainingor coarse-grainingtechniques to capture differentfeatures and phenomena Chemically-specificcoarse-grainedmodels can be derived from atomistic models with techniques like iterative Boltzmann inversionand forcematching U = U bond + U angle + U dihedral + U nonbond

Upload: others

Post on 20-Jun-2020

1 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Multiscale Molecular Simulations of PV Encapsulants · 2018-08-08 · Multiscale Molecular Simulations of PV Encapsulants Lauren Abbott1, Ross Larsen2, and Stan ... • Object-oriented

Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-AC04-94AL85000.

SAND No. 2011-XXXXP

MultiscaleMolecularSimulationsofPVEncapsulantsLaurenAbbott1,RossLarsen2,andStanMoore1

1SandiaNationalLaboratories,Albuquerque,NM,2NationalRenewableEnergyLaboratory,Golden,CO

• Automatedgenerationofcandidatemoleculesformaterialsdiscovery• Automatedgenerationofinitialconditionsformolecularsimulations• Automatedextractionofsub-regionsforanalysisandcalculationofproperties• Integrationwithanalysisandvisualizationalgorithms• Object-orientedframeworkforeasyextensionofthefunctionalitytospecificneeds

SimulationToolkitforRenewableEnergyandAdvancedMaterialsModeling

streamm.nrel.gov

• Opensourceclassicalmoleculardynamicscode• Particlesimulationsattheatomic,meso,andcontinuumscales• Inclusionsofcommonpotentialsforhardandsoftmaterials• Efficientparallelsimulationsusingspatialdecompositionofdomainspace• Modularframeworkforeasyextensionwithnewfeaturesandfunctionality

Large-scaleAtomic/MolecularMassivelyParallelSimulator

lammps.sandia.gov

Polymer/ClayNanocomposites• SandiaandTexasA&Mareworkingon

inexpensive,transparentpolymer/claynanocompositeswithlayersofpolymersandorientedclayplatelets

• Tailoringtheclayparticlesandpolymermatrixcontrolsthebarrierproperties,compositeintegrity,andfireretardancy

• Chemically-specificcoarse-grainedmodelsofpolymer/claynanocompositescanyieldinsightintomaterialprocesses(e.g.,encapsulation,intercalation,andexfoliation)andresultingmaterialproperties

Clayplatelets,~50-100nmdiam

*Sandia:MargaretGordon,ErikSpoerke,EricSchindelholtz,KenArmijo,RobSorensen,TexasA&M:JamieGrunlan,KevinHolder,Shuang Qin

neatpolymer

composite

Stress-straincurveLeft:Coarse-grainedmodel;Middle:Exampleofpolymer intercalation;Right:Claydispersion inpolymermatrixwithpolymer notshown

*JLSuter,DGroen,PVCoveney,Adv.Mater.,2015,27,966-984;NanoLett.,2015,15,8108-8113

• Atomisticsimulationselucidatedpossibleeffectsofmorphologyonelectron-transportmechanismsinradicalpolymerelectrodes

• MolecularmodelsandMDsimulationsofPTMAweresetupusingSTREAMM

• Twoprimarydistances(4.5and6.5Å)wereidentifiedascontributingtoaneffectiveelectrontransferdistanceof5.5Å

• Inter-sitecouplingswerecomputedbetween>10,000pairsofsiteswithQMcalculations,showingthatelectrontransfermostlytakesplacebetweensitesondifferentchains

RadicalPolymerElectrodes

Top:N-Nradialdistribution functions;Bottom:Averageinter-siteelectroniccoupling

Figure 3: PMMA and PTMA oligomers 12, 24 and 48 monomerunits, (PMMA38PTMA10), (PMMA2PTMA8)5, (PMMA30PTMA18),(PMMA6PTMA4)5, (PMMA10PTMA38), (PMMA8PTMA8)5

10

MD

QM

PTMA

STREAMM’sautomatedstructuregenerationfrommonomertooligomertobulkmaterial,andautomatedextractionofsub-regionsforQMcalculations

*TWKemper,RELarsen,TGennett,J.Phys.Chem.C,2014,118,17213-17220

PolymerPhotovoltaicActiveLayers• AtomisticMDsimulationsweresetupusingSTREAMMandperformedin

LAMMPSforthreeorganicPVcopolymers(BDT-TPD,PTB7,andPTB7-Th),withvariationstothebackboneandside-chainstructure

• AlignedparallelchainswithpistackingwereobservedinBDT-TPDandPTB7,butnotinPTB7-Thduetosterichindrancefromtheside-chain,whichcontradictsassumptionmadeintheliterature

• TransportinBDT-TPDandPTB7likelyoccursbetweenparallelpistacks,whiletransportinPTB7-Thlikelyoccursbetweenorthogonalpistacks

BDT-TPDPTB7

PTB7-Th

Molecularstructures(top)andsimulation snapshots (bottom)fororganicPVcopolymers

*NREL:TravisKemperandRossLarsen

Ion-ContainingPolymers• Sulfonatedpolyphenylenes(SDAPP)showpromiseforprotonexchange

membranefuelcells(PEMFCs)andvanadiumflowbatteries• Littleisknownaboutthenanoscalestructureoftheseamorphouspolymers,

whichisnoteasilycharacterizedwithcurrentexperimentaltechniques• Atomisticsimulationsyieldeduniqueinsightintothemorphologyoftheionic

domainsformedbytheaggregationofwaterandionicgroups,specificallythesize,shape,andconnectivity,aswellasresultingimplicationsforiontransport

• Asmorewaterisaddedintothesystem,theionicdomainsbecomemorefullypercolatedandthedomainsbecomeslightlylargerandmoresphericalinshape,whichwouldimproveiontransport

increasingwater

Left:MolecularstructureofSDAPP;Right:Ionicdomainswithdisparateclustersshown indifferentcolors

impedance composed of a real component, Z′, and an imaginarycomponent, Z′′. An example of the raw data is shown in Figure3. To compute the membrane proton conductivity from thecomplex impedance response, the impedance line is extrapo-lated to the x-axis. The extrapolated value of the real imped-ance where the imaginary response is zero (Z′ at Z′′ ) 0) isthen taken as the resistance of the membrane, and eq 3 is usedto compute the membrane proton conductivity

where L is the length between the sense electrodes, Z′ is thereal part of the impedance response (extrapolated to Z′′ ) 0),and A is the area available for proton conduction (width ×thickness). All proton conductivities reported here were mea-sured with the film immersed in liquid water at 30 °C duringthe measurement time.

The tensile stress-strain properties of the SDAPP3 mem-branes were compared with Nafion 117 in order to evaluatethe mechanical strength of a typical sulfonated poly(phenylene)vs the state-of-the-art membrane. The samples were testedusing a Com-Ten Industries 95T series load frame equippedwith a 200 lbf load cell and computerized data acquisitionsoftware. Samples of 9 mm width were deformed at a cross-head speed of 5 mm/min with gauge length of 30 mm. Reportedstress-strain measurements are the average of at least threetests.

Results and DiscussionSynthesis and Characterization. Preparation of

Diels-Alder poly(phenylene)s has been described else-where16 and briefly described above (vide supra). Sul-fonation of SDAPP0 has not been reported, but substi-

tuted 1,4-poly(phenylene) has been sulfonated withconcentrated sulfuric acid with levels of sulfonationcontrolled by varying the reaction time.11 Our approachto sulfonating SDAPP0 was the creation of a 6 wt %polymer solution of SDAPP0 in methylene chloride thatwas subsequently homogeneously sulfonated in situwith chlorosulfonic acid.

The repeat unit of SDAPP0 with its six pendentphenyl groups (Scheme 2) provides a number of possiblesites for sulfonation. However, sulfonation is thoughtto occur predominantly at the para positions of thependant phenyl groups due to their positioning17 aboutthe sterically congested, concoplanar, rigid-rod backbone(Figure 1). Therefore, the limits of sulfonation areprojected to be between 0 and 6 sulfonic acid groups perrepeat unit. By varying the ratio of moles of chlorosul-fonic acid to moles of polymer repeat unit charged tothe sulfonation reaction, 0.8-2.1 sulfonic acid groupsper repeat unit were achieved as measured by titration(see below).

Introduction of increasing numbers of sulfonic acidgroups onto polymers improves their ionic conductivity

Figure 2. Schematic of four-point membrane proton conductivity cell.

Figure 3. Impedance response of a typical proton conductingmembrane between 100 kHz and 100 Hz.

κ ) LZ′A

(3)

Scheme 2. Sulfonation of Diels-AlderPoly(phenylene)

5012 Fujimoto et al. Macromolecules, Vol. 38, No. 12, 2005

SDAPP

*Sandia:LaurenAbbott, AmalieFrischknecht,CyFujimoto,ToddAlam,EricSorte

• MolecularsimulationscanprovideimportantinformationaboutnanoscalestructureandphenomenainPVencapsulantmaterialswithatomic-leveldetailnotcurrentlyaccessibletoexperimentaltechniques

• Theatomic-leveldetailofmolecularsimulationsareusefulforstudyinglocalpackingofthepolymer,aswellastheinteractionsanddynamicsofadditives(e.g.,UVabsorbers)andcontaminants(e.g.,H2OorO2)withinthepolymer

• Molecularsimulationsareacost-effectiverouteforefficienthigh-throughputscreeningofpotentialencapsulantmaterials

• Amultiscaleapproachcombinesquantummechanical(QM)calculationsandatomisticandcoarse-grainedmoleculardynamics(MD)simulationstocaptureabroadrangeoflengthandtimescales

• Resultsfromthemolecularsimulationscanbepassedtosimulationmodelsathigherlengthsscales,suchasmeso orcontinuumtechniques

• ThiscapabilityleveragesopensourcetoolslikeSandia’sLAMMPSMDpackageandNREL’sSTREAMMtoolkittosetup,run,andanalyzemolecularsimulations

Overview• MDtechniquesfollowclassicaldynamicsusingNewton’slaw(F =ma),by

updatingthepositionsofparticlesusingthenetforcesontheparticles• Potentialsaredescribedusingforcefieldswithtermsforbondedinteractions

(e.g.,bondsandangles)andnonbondedinteractions(e.g.,vanderWaals)

MolecularModeling

• QMtechniquesapproximatethewavefunctionfollowingSchrödinger’sequation(Hψ =iħ ∂ψ/∂t)toconsiderquantumeffects

• Wecanmovebetweendifferentscalesusingfine-grainingorcoarse-grainingtechniquestocapturedifferentfeaturesandphenomena

• Chemically-specificcoarse-grainedmodelscanbederivedfromatomisticmodelswithtechniqueslikeiterativeBoltzmanninversionandforcematching

U =Ubond +Uangle +Udihedral +Unonbond