Coarse-Grained Dynamically Accurate Simulations of Ionic Liquids: [pyr14][TFSI] and [EMIM][BF4].

In this work, coarse-grained (CG) models for two different sets of ionic liquids were developed from atomistic molecular dynamics (MD) reference systems, expanding their system size and time duration capabilities. The bonded force field of the CG systems was built using harmonic oscillator potential (HOP) fitting, while the nonbonded force field was generated with the multiscale coarse-graining (MS-CG) approach based on force matching. The dynamics of each system were corrected using the probability distribution function-based coarse-grained molecular dynamics (PDF-based CGMD) method. The structure and dynamics of each system were proven to match reference system data at two temperature scales. CG models and force fields for these liquids were developed to exemplify a general purpose methodology for producing MD results of ionic liquids and other fluids with accurate structural as well as dynamic properties. As an application, developed ionic liquids CG models were then applied to study vacuum-interface interaction. Density profile results of vacuum-interface exposure show significant deviation from bulk behavior. At the interface, multilayer ordering of ionic liquids is predicted to be similar to those observed from an experimental work. This ordering is intensified by decreasing temperature and use of the PDF-based CGMD method as opposed to conventional CG methods.

Proton Abstraction from DME n···X+ by OH-, O2-, and XO2-, for X = Li, Na, and K: Implications for Li-O2 Batteries.

The abstraction of a proton by OH-, O2-, and XO2- from DME n···X+, where X is Li, Na, or K, is studied using density functional theory. Both the gas phase and the solution phase are studied. In general, when explicit solvent molecules are added, the difference between the gas-phase and solution results becomes rather small. While the DME n···X+ binding energies differ significantly for various alkali cations, the reaction energies and transition-state energies are far less sensitive to the choice of an alkali cation. XO2- has a lower barrier height than OH-, which, in turn, has a lower barrier height than O2-. The reaction energies follow the same trends.

Polarizable Molecular Dynamics and Experiments of 1,2-Dimethoxyethane Electrolytes with Lithium and Sodium Salts: Structure and Transport Properties.

The structure and transport properties of electrolyte solutions of 1,2-dimethoxyethane (DME) having salts of Li+ with bis(trifluoromethanesulfonyl)imide ([TFSI]-) or Na+ with [TFSI]- are investigated with polarizable molecular dynamics and experiments. Polarizable force fields for Li+ and Na+ with DME and [TFSI]- were developed based on quantum chemistry calculations, ab initio molecular dynamics simulations, and thermodynamic liquid-state properties. Simulation results for density, viscosity, self-diffusion coefficient, and conductivity of the electrolytes all agree well with the trends and magnitudes of available experimental data for a wide range of salt concentrations. As the concentration of salt increases, the electrolytes become more viscous and molecular species become less mobile. Ionic conductivity does not change monotonically with salt concentration and exhibits a maximum between 0.5 and 1.0 M for both Li[TFSI] and Na[TFSI] electrolytes. Comparatively, both cations are solvated by 5-6 DME or [TFSI]- oxygen atoms and exhibit similar diffusivities and conductivities. The solvation shell of Na+ is found to be more weakly bound and to have a lower binding residence time than that of Li+. The transport of Li+ therefore is more vehicular, through the motion of the solvation shell, while the transport of Na+ is based more on exchange, through the replacement of solvating species. The atomistic insight provided by this work can be used as the basis for future rational design of improved electrolyte solvents for lithium-oxygen, sodium-oxygen, and lithium-sulfur batteries.

Phenolic Polymer Solvation in Water and Ethylene Glycol, I: Molecular Dynamics Simulations.

Interactions between pre-cured phenolic polymer chains and a solvent have a significant impact on the structure and properties of the final postcured phenolic resin. Developing an understanding of the nature of these interactions is important and will aid in the selection of the proper solvent that will lead to the desired final product. Here, we investigate the role of the phenolic chain structure and the solvent type on the overall solvation performance of the system through molecular dynamics simulations. Two types of solvents are considered: ethylene glycol (EGL) and H2O. In addition, three phenolic chain structures are considered, including two novolac-type chains with either an ortho-ortho (OON) or an ortho-para (OPN) backbone network and a resole-type (RES) chain with an ortho-ortho network. Each system is characterized through a structural analysis of the solvation shell and the hydrogen-bonding environment as well as through a quantification of the solvation free energy along with partitioned interaction energies between specific molecular species. The combination of simulations and the analyses indicate that EGL provides a higher solvation free energy than H2O due to more energetically favorable hydrophilic interactions as well as favorable hydrophobic interactions between CH element groups. In addition, the phenolic chain structure significantly affects the solvation performance, with OON having limited intermolecular hydrogen-bond formations, while OPN and RES interact more favorably with the solvent molecules. The results suggest that a resole-type phenolic chain with an ortho-para network should have the best solvation performance in EGL, H2O, and other similar solvents.

Phenolic Polymer Solvation in Water and Ethylene Glycol, II: Ab Initio Computations.

Ab initio techniques are used to study the interaction of ethylene glycol and water with a phenolic polymer. The water bonds more strongly with the phenolic OH than with the ring. The phenolic OH groups can form hydrogen bonds between themselves. For more than one water molecule, there is a competition between water-water and water-phenolic interactions. Ethylene glycol shows the same effects as those of water, but the potential energy surface is further complicated by CH2-phenolic interactions, different conformers of ethylene glycol, and two OH groups on each molecule. Thus, the ethylene glycol-phenolic potential is more complicated than the water-phenolic potential. The results of the ab initio calculations are compared to those obtained using a force field. These calibration studies show that the water system is easier to describe than the ethylene glycol system. The calibration studies confirm the reliability of force fields used in our companion molecular dynamics study of a phenolic polymer in water and ethylene solutions.

Finite Temperature Properties of NiTi from First Principles Simulations: Structure, Mechanics, and Thermodynamics.

We present a procedure to determine temperature-dependent thermodynamic properties of crystalline materials from density functional theory molecular dynamics (DFT-MD). Finite temperature properties (structural, thermal, and mechanical properties) of the phases (ground state monoclinic B33, martensitic B19', and austenitic B2) of the shape memory alloy NiTi are investigated. Fluctuation formulas and numerical derivatives are used to evaluate mechanical and thermal properties. A modified version of thermodynamic upsampling is used to enable simulations with lower DFT convergence thresholds. In addition, a thermodynamic integration expression is developed to compute free energies from isobaric DFT-MD simulations that accounts for volume changes. Structural parameters, elastic constants, volume expansion, and specific heats as a function of temperature are evaluated. Phase transitions between B2 and B19' and between B19' and B33 are characterized according to their thermal energy, entropy, and free energy differences as well as their latent heats. Anharmonic effects are shown to play a large role in both stabilizing the austenite B2 phase as well as suppressing the martensitic phase transition. The quasiharmonic approximation to the free energy results in large errors in estimating the martensitic transition temperature by neglecting these large anharmonic components.

Evaluation of molecular dynamics simulation methods for ionic liquid electric double layers.

We investigate how systematically increasing the accuracy of various molecular dynamics modeling techniques influences the structure and capacitance of ionic liquid electric double layers (EDLs). The techniques probed concern long-range electrostatic interactions, electrode charging (constant charge versus constant potential conditions), and electrolyte polarizability. Our simulations are performed on a quasi-two-dimensional, or slab-like, model capacitor, which is composed of a polarizable ionic liquid electrolyte, [EMIM][BF4], interfaced between two graphite electrodes. To ensure an accurate representation of EDL differential capacitance, we derive new fluctuation formulas that resolve the differential capacitance as a function of electrode charge or electrode potential. The magnitude of differential capacitance shows sensitivity to different long-range electrostatic summation techniques, while the shape of differential capacitance is affected by charging technique and the polarizability of the electrolyte. For long-range summation techniques, errors in magnitude can be mitigated by employing two-dimensional or corrected three dimensional electrostatic summations, which led to electric fields that conform to those of a classical electrostatic parallel plate capacitor. With respect to charging, the changes in shape are a result of ions in the Stern layer (i.e., ions at the electrode surface) having a higher electrostatic affinity to constant potential electrodes than to constant charge electrodes. For electrolyte polarizability, shape changes originate from induced dipoles that soften the interaction of Stern layer ions with the electrode. The softening is traced to ion correlations vertical to the electrode surface that induce dipoles that oppose double layer formation. In general, our analysis indicates an accuracy dependent differential capacitance profile that transitions from the characteristic camel shape with coarser representations to a more diffuse profile with finer representations.

Computational and Experimental Study of Li-Doped Ionic Liquids at Electrified Interfaces.

We evaluate the influence of Li-salt doping on the dynamics, capacitance, and structure of three ionic liquid electrolytes, [pyr14][TFSI], [pyr13][FSI], and [EMIM][BF4], using molecular dynamics and polarizable force fields. In this respect, our focus is on the properties of the electric double layer (EDL) formed by the electrolytes at the electrode surface as a function of surface potential (Ψ). The rates of EDL formation are found to be on the order of hundreds of picoseconds and only slightly influenced by the addition of Li-salt. The EDLs of three electrolytes are shown to have different energy storage capacities, which we relate to the EDL formation free energy. The differential capacitance obtained from our computations exhibits asymmetry about the potential of zero charge and is consistent with the camel-like profiles noted from mean field theories and experiments on metallic electrodes. The introduction of Li-salt reduces the noted asymmetry in the differential capacitance profile. Complementary experimental capacitance measurements have been made on our three electrolytes in their neat forms and with Li-salt. The measurements, performed on glassy carbon electrodes, produce U-like profiles, and Li-salt doping is shown to strongly affect capacitance at high magnitudes of Ψ. Differences in the theoretical and experimental shapes and magnitudes of capacitance are rationalized in terms of the electrode surface and pseudocapacitive effects. In both neat and Li-doped liquids, the details of the computational capacitance profile are well described by Ψ-induced changes in the density and molecular orientation of ions in the molecular layer closest to the electrode. Our results suggest that the addition of Li+ induces disorder in the EDL, which originates from the strong binding of anions to Li+. An in-depth analysis of the distribution of Li+ in the EDL reveals that it does not readily enter the molecular layer at the electrode surface, preferring instead to be localized farther away from the surface in the second molecular layer. This behavior is validated through an analysis of the free energy of Li+ solvation as a function of distance from the electrode. Free energy wells are found to coincide with localized concentrations of Li+, the depths of which increase with Ψ and suggest a source of impedance for Li+ to reach the electrode. Finally, we make predictions of the specific energy at ideal graphite utilizing the computed capacitance and previously derived electrochemical windows of the liquids.

Ab Initio Simulations and Electronic Structure of Lithium-Doped Ionic Liquids: Structure, Transport, and Electrochemical Stability.

Density functional theory (DFT), density functional theory molecular dynamics (DFT-MD), and classical molecular dynamics using polarizable force fields (PFF-MD) are employed to evaluate the influence of Li(+) on the structure, transport, and electrochemical stability of three potential ionic liquid electrolytes: N-methyl-N-butylpyrrolidinium bis(trifluoromethanesulfonyl)imide ([pyr14][TFSI]), N-methyl-N-propylpyrrolidinium bis(fluorosulfonyl)imide ([pyr13][FSI]), and 1-ethyl-3-methylimidazolium boron tetrafluoride ([EMIM][BF4]). We characterize the Li(+) solvation shell through DFT computations of [Li(Anion)n]((n-1)-) clusters, DFT-MD simulations of isolated Li(+) in small ionic liquid systems, and PFF-MD simulations with high Li-doping levels in large ionic liquid systems. At low levels of Li-salt doping, highly stable solvation shells having two to three anions are seen in both [pyr14][TFSI] and [pyr13][FSI], whereas solvation shells with four anions dominate in [EMIM][BF4]. At higher levels of doping, we find the formation of complex Li-network structures that increase the frequency of four anion-coordinated solvation shells. A comparison of computational and experimental Raman spectra for a wide range of [Li(Anion)n]((n-1)-) clusters shows that our proposed structures are consistent with experiment. We then compute the ion diffusion coefficients and find measures from small-cell DFT-MD simulations to be the correct order of magnitude, but influenced by small system size and short simulation length. Correcting for these errors with complementary PFF-MD simulations, we find DFT-MD measures to be in close agreement with experiment. Finally, we compute electrochemical windows from DFT computations on isolated ions, interacting cation/anion pairs, and liquid-phase systems with Li-doping. For the molecular-level computations, we generally find the difference between ionization energy and electron affinity from isolated ions and interacting cation/anion pairs to provide upper and lower bounds, respectively, to experiment. In the liquid phase, we find the difference between the lowest unoccupied and highest occupied electronic levels in pure and hybrid functionals to provide lower and upper bounds, respectively, to experiment. Li-doping in the liquid-phase systems results in electrochemical windows little changed from the neat systems.

Structure and energetics of Li(+)-(BF4(-))n, Li(+)-(FSI(-))n, and Li(+)-(TFSI(-))n: ab initio and polarizable force field approaches.

The Li(+)-BF4(-) and BF4(-)-BF4(-) interactions are studied using second order perturbation theory (MP2) and coupled cluster singles and doubles approach, including the effect of connected triples, CCSD(T). The MP2 and CCSD(T) results are in excellent agreement. Using only the MP2 approach, the interactions of Li(+) with bis(trifluoromethane)sulfonimide anion (TFSI) and Li(+) with bis(fluorosulfonyl)imide anion (FSI) are studied. The results of these high level calculations are compared with density functional theory (DFT) calculations for a variety of functionals and with the APPLE&P force field. The B3LYP approach well reproduces the accurate calculations using both a small and large basis set. The M06 and M06L functionals in the larger basis set are in good agreement with the high level calculations. While the APPLE&P force field does not outperform the best functionals, the APPLE&P results agree better with the accurate results than do some of the functionals tested.

Computational and experimental investigation of Li-doped ionic liquid electrolytes: [pyr14][TFSI], [pyr13][FSI], and [EMIM][BF4].

We employ molecular dynamics (MD) simulation and experiment to investigate the structure, thermodynamics, and transport of N-methyl-N-butylpyrrolidinium bis(trifluoromethylsufonyl)imide ([pyr14][TFSI]), N-methyl-N-propylpyrrolidinium bis(fluorosufonyl)imide ([pyr13][FSI]), and 1-ethyl-3-methylimidazolium boron tetrafluoride ([EMIM][BF4]), as a function of Li-salt mole fraction (0.05 ≤ xLi(+) ≤ 0.33) and temperature (298 K ≤ T ≤ 393 K). Structurally, Li(+) is shown to be solvated by three anion neighbors in [pyr14][TFSI] and four anion neighbors in both [pyr13][FSI] and [EMIM][BF4], and at all levels of xLi(+) we find the presence of lithium aggregates. Pulsed field gradient spin-echo NMR measurements of diffusion and electrochemical impedance spectroscopy measurements of ionic conductivity are made for the neat ionic liquids as well as 0.5 molal solutions of Li-salt in the ionic liquids. Bulk ionic liquid properties (density, diffusion, viscosity, and ionic conductivity) are obtained with MD simulations and show excellent agreement with experiment. While the diffusion exhibits a systematic decrease with increasing xLi(+), the contribution of Li(+) to ionic conductivity increases until reaching a saturation doping level of xLi(+) = 0.10. Comparatively, the Li(+) conductivity of [pyr14][TFSI] is an order of magnitude lower than that of the other liquids, which range between 0.1 and 0.3 mS/cm. Our transport results also demonstrate the necessity of long MD simulation runs (∼200 ns) to converge transport properties at room temperature. The differences in Li(+) transport are reflected in the residence times of Li(+) with the anions (τ(Li/-)), which are revealed to be much larger for [pyr14][TFSI] (up to 100 ns at the highest doping levels) than in either [EMIM][BF4] or [pyr13][FSI]. Finally, to comment on the relative kinetics of Li(+) transport in each liquid, we find that while the net motion of Li(+) with its solvation shell (vehicular) significantly contributes to net diffusion in all liquids, the importance of transport through anion exchange increases at high xLi(+) and in liquids with large anions.

Equilibrium limit of thermal conduction and boundary scattering in nanostructures.

Determining the lattice thermal conductivity (κ) of nanostructures is especially challenging in that, aside from the phonon-phonon scattering present in large systems, the scattering of phonons from the system boundary greatly influences heat transport, particularly when system length (L) is less than the average phonon mean free path (MFP). One possible route to modeling κ in these systems is through molecular dynamics (MD) simulations, inherently including both phonon-phonon and phonon-boundary scattering effects in the classical limit. Here, we compare current MD methods for computing κ in nanostructures with both L ⩽ MFP and L ≫ MFP, referred to as mean free path constrained (cMFP) and unconstrained (uMFP), respectively. Using a (10,0) CNT (carbon nanotube) as a benchmark case, we find that while the uMFP limit of κ is well-defined through the use of equilibrium MD and the time-correlation formalism, the standard equilibrium procedure for κ is not appropriate for the treatment of the cMFP limit because of the large influence of boundary scattering. To address this issue, we define an appropriate equilibrium procedure for cMFP systems that, through comparison to high-fidelity non-equilibrium methods, is shown to be the low thermal gradient limit to non-equilibrium results. Further, as a means of predicting κ in systems having L ≫ MFP from cMFP results, we employ an extrapolation procedure based on the phenomenological, boundary scattering inclusive expression of Callaway [Phys. Rev. 113, 1046 (1959)]. Using κ from systems with L ⩽ 3 µm in the extrapolation, we find that the equilibrium uMFP κ of a (10,0) CNT can be predicted within 5%. The equilibrium procedure is then applied to a variety of carbon-based nanostructures, such as graphene flakes (GF), graphene nanoribbons (GNRs), CNTs, and icosahedral fullerenes, to determine the influence of size and environment (suspended versus supported) on κ. Concerning the GF and GNR systems, we find that the supported samples yield consistently lower values of κ and that the phonon-boundary scattering remains dominant at large lengths, with L = 0.4 µm structures exhibiting a third of the periodic result. We finally characterize the effect of shape in CNTs and fullerenes on κ, showing the angular components of conductivity in CNTs and icosahedral fullerenes are similar for a given circumference.

Control of thermal and electronic transport in defect-engineered graphene nanoribbons.

The influence of the structural detail and defects on the thermal and electronic transport properties of graphene nanoribbons (GNRs) is explored by molecular dynamics and non-equilibrium Green's function methods. A variety of randomly oriented and distributed defects, single and double vacancies, Stone-Wales defects, as well as two types of edge form (armchair and zigzag) and different edge roughnesses are studied for model systems similar in sizes to experiments (>100 nm long and >15 nm wide). We observe substantial reduction in thermal conductivity due to all forms of defects, whereas electrical conductance reveals a peculiar defect-type-dependent response. We find that a 0.1% single vacancy concentration and a 0.23% double vacancy or Stone-Wales concentration lead to a drastic reduction in thermal conductivity of GNRs, namely, an 80% reduction from the pristine one of the same width. Edge roughness with an rms value of 7.28 Å leads to a similar reduction in thermal conductivity. Randomly distributed bulk vacancies are also found to strongly suppress the ballistic nature of electrons and reduce the conductance by 2 orders of magnitude. However, we have identified that defects close to the edges and relatively small values of edge roughness preserve the quasi-ballistic nature of electronic transport. This presents a route of independently controlling electrical and thermal transport by judicious engineering of the defect distribution; we discuss the implications of this for thermoelectric performance.