Simulation of polymeric mixed ionic and electronic conductors with a combined classical and quantum mechanical model.

In organic polymeric materials with mixed ionic and electronic conduction (OMIEC), the excess charge in doped polymers is very mobile and the dynamics of the polymer chain cannot be accurately described with a model including only fixed point charges. Ions and polymer are comparatively slower and a methodology to capture the correlated motions of excess charge and ions is currently unavailable. Considering a prototypical interface encountered in this type of materials, we constructed a scheme based on the combination of MD and QM/MM to evaluate the classical dynamics of polymer, water and ions, while allowing the excess charge of the polymer chains to rearrange following the external electrostatic potential. We find that the location of the excess charge varies substantially between chains. The excess charge changes across multiple timescales as a result of fast structural fluctuations and slow rearrangement of the polymeric chains. Our results indicate that such effects are likely important to describe the phenomenology of OMIEC, but additional features should be added to the model to enable the study of processes such as electrochemical doping.

Development of hybrid coarse-grained atomistic models for rapid assessment of local structuring of polymeric semiconductors.

Decades of work in the field of computational study of semiconducting polymers using atomistic models illustrate the challenges of generating equilibrated models for this class of materials. While adopting a coarse-grained model can be helpful, the process of developing a suitable model is particularly non-trivial and time-consuming for semiconducting polymers due to a large number of different interactions with some having an anisotropic nature. This work introduces a procedure for the rapid generation of a hybrid model for semiconducting polymers where atoms of secondary importance (those in the alkyl side chains) are transformed into coarse-grained beads to reduce the computational cost of generating an equilibrated structure. The parameters are determined from easy-to-equilibrate simulations of very short oligomers and the model is constructed to enable a very simple back-mapping procedure to reconstruct geometries with atomistic resolution. The model is illustrated for three related polymers containing DPP (diketopyrrolopyrrole) to evaluate the transferability of the potential across different families of polymers. The accuracy of the model, determined by comparison with the results of fully equilibrated simulations of the same material before and after back-mapping, is fully satisfactory for two out of the three cases considered. We noticed that accuracy can be determined very early in the workflow so that it is easy to assess when the deployment of this method is advantageous. The hybrid representation can be used to evaluate directly the electronic properties of structures sampled by the simulations.

Investigating the quasi-liquid layer on ice surfaces: a comparison of order parameters.

Ice surfaces are characterized by pre-melted quasi-liquid layers (QLLs), which mediate both crystal growth processes and interactions with external agents. Understanding QLLs at the molecular level is necessary to unravel the mechanisms of ice crystal formation. Computational studies of the QLLs heavily rely on the accuracy of the methods employed for identifying the local molecular environment and arrangements, discriminating between solid-like and liquid-like water molecules. Here we compare the results obtained using different order parameters to characterize the QLLs on hexagonal ice (Ih) and cubic ice (Ic) model surfaces investigated with molecular dynamics (MD) simulations in a range of temperatures. For the classification task, in addition to the traditional Steinhardt order parameters in different flavours, we select an entropy fingerprint and a deep learning neural network approach (DeepIce), which are conceptually different methodologies. We find that all the analysis methods give qualitatively similar trends for the behaviours of the QLLs on ice surfaces with temperature, with some subtle differences in the classification sensitivity limited to the solid-liquid interface. The thickness of QLLs on the ice surface increases gradually as the temperature increases. The trends of the QLL size and of the values of the order parameters as a function of temperature for the different facets may be linked to surface growth rates which, in turn, affect crystal morphologies at lower vapour pressure. The choice of the order parameter can be therefore informed by computational convenience except in cases where a very accurate determination of the liquid-solid interface is important.

Local structuring of diketopyrrolopyrrole (DPP)-based oligomers from molecular dynamics simulations.

The microscopic structure of high mobility semiconducting polymers is known to be essential for their performance but it cannot be easily deduced from the available experimental data. A series of short oligomers of diketopyrrolopyrrole (DPP)-based materials that display high charge mobility are studied by molecular dynamics simulations to understand their local structuring at an atomic level. Different analyses are proposed to compare the ability of different oligomers to form large aggregates and their driving force. The simulations show that the tendency for this class of materials to form aggregates is driven by the interaction between DPP fragments, but this is modulated by the other conjugated fragments of the materials which affect the rigidity of the polymer and, ultimately, the size of the aggregates that are formed. The main structural features and the electronic structure of the oligomers are fairly similar above the glass transition temperature and at room temperature.