Search

Recent theories suggest that low frequency dynamic intramolecular and intermolecular motions in organic semiconductors (OSCs) are critical to determining the hole mobility. So far, however, it has not been possible to probe these motions directly experimentally and therefore no unequivocal and quantitative link exists between molecular-scale thermal disorder and macroscale hole mobility in OSCs. Here we use inelastic neutron scattering to probe thermal disorder directly by measuring the phonon spectrum in six different small molecule OSCs, which we accurately reproduce with first principles simulations. We use the simulated phonons to generate a set of electron-phonon coupling parameters. Using these parameters, the theoretical mobility is in excellent agreement with macroscopic measurements. Comparison of mobility between different materials reveals routes to improve mobility by engineering phonon and electron-phonon coupling.

When an organic semiconductor (OSC) is blended with an electron acceptor molecule that can act as a p-type dopant, there should ideally be complete (integer) transfer of charge from the OSC to the dopant. However, some dopant-OSC blends instead form charge transfer complexes (CTCs), characterized by fractional charge transfer (CT) and strong orbital hybridization between the two molecules. Fractional CT doping does not efficiently generate free charge carriers, but it is unclear what conditions lead to incomplete charge transfer. Here we show that by modifying film processing conditions in the semiconductor-dopant couple poly(3-hexylthiophene):2,3,5,6-tetrafluoro-7,7,8,8,-tetracyanoquinodimethane (P3HT:F4TCNQ), we can selectively obtain nearly pure integer or fractional CT phases. Fractional CT films show electrical conductivities approximately 2 orders of magnitude lower than corresponding integer CT films, and remarkably different optical absorption spectra. Grazing incidence wide-angle X-ray diffraction (GIXD) reveals that fractional CT films display an unusually dense and well-ordered crystal structure. These films show lower paracrystallinity and shorter lamellar and π-stacking distances than undoped films processed under similar conditions. Using plane-wave DFT we obtain a structure with unit cell parameters closely matching those observed by GIXD. This first-ever observation of both fractional and integer CT in a single OSC-dopant system demonstrates the importance of structural effects on OSC doping and opens the door to further studies.

Carbon nanoscrolls (CNSs), with controllable interlayer distances, have been systematically studied with respect to noble gas adsorption, separation, and transport properties. Grand canonical Monte Carlo simulations were performed on pure noble gases and noble gas mixtures interacting with CNSs at different temperatures. There is separation of gases across a series of CNS interlayer distances at different pressures. The optimal interlayer distance for each noble gas was determined. Selective adsorption was characterized for noble gas mixtures, and the corresponding selectivity was explored. A molecular sieving effect was observed for a certain interlayer gap range, in which the energetically unfavorable smaller gas atoms were better stored in the CNS. In addition, molecular dynamics simulations were conducted to understand the transport properties of noble gases in carbon nanoscrolls in both the cross-sectional and axial directions of the CNS. Gas molecules are found to be captured most strongly when the CNS has the optimal interlayer distance for adsorption.

Conformational and energetic disorder in organic semiconductors reduces charge and exciton transport because of the structural defects, thus reducing the efficiency in devices such as organic photovoltaics and organic light-emitting diodes. The main structural heterogeneity is because of the twisting of the polymer backbone that occurs even in polymers that are mostly crystalline. Here, we explore the relationship between polymer backbone twisting and exciton delocalization by means of transient absorption spectroscopy and density functional theory calculations. We study the PffBT4T-2DT polymer which has exhibited even higher device efficiency with nonfullerene acceptors than the current record breaking PCE11 polymer. We determine the driving force for planarization of a polymer chain caused by excitation. The methodology is generally applicable and demonstrates a higher penalty for nonplanar structures in the excited state than in the ground state. This study highlights the morphological and electronic changes in conjugated polymers that are brought about by excitation.

Helium gas accumulation from alpha decay during extended storage of spent fuel has potential to compromise the structural integrity the fuel. Here we report results obtained with surrogate nickel particles which suggest that alumina formed by atomic layer deposition can serve as a low volume-fraction, uniformly-distributed phase for retention of helium generated in fuel particles such as uranium oxide. Thin alumina layers may also form transport paths for helium in the fuel rod, which would otherwise be impermeable. Micron-scale nickel particles, representative of uranium oxide particles in their low helium solubility and compatibility with the alumina synthesis process, were homogeneously coated with alumina approximately 3–20 nm by particle atomic layer deposition (ALD) using a fluidized bed reactor. Particles were then loaded with helium at 800 °C in a tube furnace. Subsequent helium spectroscopy measurements showed that the alumina phase, or more likely a related nickel/alumina interface structure, retains helium at a density of at least 1017 atoms/cm3. High resolution transmission electron microscopy revealed that the thermal treatment increased the alumina thickness and generated additional porosity. Results from Monte Carlo simulations on amorphous alumina predict the helium retention concentration at room temperature could reach 1021 atoms/cm3 at 400 MPa, a pressure predicted by others to be developed in uranium oxide without an alumina secondary phase. This concentration is sufficient to eliminate bubble formation in the nuclear fuel for long-term storage scenarios, for example. Measurements by others of the diffusion coefficient in polycrystalline alumina indicate values several orders of magnitude higher than in uranium oxide, which then can also allow for helium transport out of the spent fuel.

Organic electronics is a popular and rapidly growing field of research. The optical, electrical and mechanical properties of organic molecules and materials can be tailored using increasingly well controlled synthetic methods. The challenge and fascination with this field of research is derived from the fact that not only the chemical identity, but also the spatial arrangement of the molecules critically affects the performance of the material. Thus synthetic, fabrication, characterisation and computational scientists need to work closely to relate a materials performance in a device to the molecular details that cause and optimise that performance. For computational scientists in particular, the need to relate macroscopic device performance to details of molecular electronic structure brings challenges in methodology due to the need to bridge many orders of time and length scales. This article provides a survey of computational methods applied to multiple-length and time scale problems in organic electronic materials. Here we seek to highlight a few particular approaches that expand the simulation toolbox.

The greatest advantage of organic materials is the ability to synthetically tune desired properties. However, structural heterogeneity often obfuscates the relationship between chemical structure and functional properties. Inelastic neutron scattering (INS) is sensitive to both local structure and chemical environment and provides atomic level details that cannot be obtained through other spectroscopic or diffraction methods. INS data are composed of a density of vibrational states with no selection rules, which means that every structural configuration is equally weighted in the spectrum. This allows the INS spectrum to be quantitatively decomposed into different structural motifs. We present INS measurements of the semiconducting polymer P3HT doped with F4TCNQ supported by density functional theory calculations to identify two dominant families of undoped crystalline structures and one dominant doped structural motif, in spite of considerable heterogeneity. The differences between the undoped and doped structures indicate that P3HT side chains flatten upon doping.

The supplementary material, in the form of one pdf file, is available from the journal web page at http://dx.doi.org/10.1140/epjst/e2016-02676-9 The Publishers regret this mistake and apologize for any inconvenience caused.

The study of protein mass transport via atomistic simulation requires time and length scales beyond the computational capabilities of modern computer systems. The raspberry model for colloidal particles in combination with the mesoscopic hydrodynamic method of lattice Boltzmann facilitates coarse-grained simulations that are on the order of microseconds and hundreds of nanometers for the study of diffusive transport of protein-like colloid particles. The raspberry model reproduces linearity in resistance to motion versus particle size and correct enhanced drag within cylindrical pores at off-center coordinates for spherical particles. Owing to the high aspect ratio of many proteins, ellipsoidal raspberry colloid particles were constructed and reproduced the geometric resistance factors of Perrin and of Happel and Brenner in the laboratory-frame and in the moving body-frame. Accurate body-frame rotations during diffusive motion have been captured for the first time using projections of displacements. The spatial discretization of the fluid leads to a renormalization of the hydrodynamic radius, however, the data describes a self-consistent hydrodynamic frame within this renormalized system.

A coarse-grained (CG) molecular simulation model has been refined for poly(2,6-dimethyl-1,4-phenylene ether) (PPE). This was successfully validated against atomistic simulation and experimental data. Particularly, the glass transition temperature (Tg) of PPE was studied using both atomistic and CG models and compared favourably to experimental data. In addition, we used the CG model together with an existing Martini CG model of polystyrene (PS) to study the blending behaviour of these two polymers. We solved the problem to mix the different potentials and molecular dynamics of high-molecular-weight blends of PPE/PS was performed in detail.

We use Molecular Dynamics and Monte Carlo simulations to examine molecular transport phenomena in nanochannels, explaining four orders of magnitude difference in wheat germ agglutinin (WGA) protein diffusion rates observed by fluorescence correlation spectroscopy (FCS) and by direct imaging of fluorescently-labeled proteins. We first use the ESPResSo Molecular Dynamics code to estimate the surface transport distance for neutral and charged proteins. We then employ a Monte Carlo model to calculate the paths of protein molecules on surfaces and in the bulk liquid transport medium. Our results show that the transport characteristics depend strongly on the degree of molecular surface coverage. Atomic force microscope characterization of surfaces exposed to WGA proteins for 1000 s show large protein aggregates consistent with the predicted coverage. These calculations and experiments provide useful insight into the details of molecular motion in confined geometries.

To meet the huge demand for renewable energy, significant research effort focuses on creating efficient organic photovoltaic (OPV) devices. In comparison to silicon-based semiconductors, OPV materials have many superior properties such as cost effectiveness, being lightweight, and flexibility, which lead to a high potential for the replacement of silicon-based semiconductors. Recently, a large number of new alternating copolymer materials have demonstrated high power conversion efficiency (PCE). These successful polymers typically have low long-range order but a high hole mobility which directly affects the PCE which depends on polymer structure. In this study, a solution molecular model for poly[N-9″-hepta-decanyl-2,7-carbazole-alt-5,5-(4′,7′-di-2-thienyl-2′,1′,3′-benzothiadiazole)]:[6,6]-phenyl (PCDTBT) is developed and subsequently a molecular dynamics simulation conducted in order to understand the structure of the polymer solution. The simulation results are consistent with a low-solubility polymer that requires long equilibration times to planarize. The structural addition of side chains to inhibit rotation of thiophene rings could improve the conjugation and processability of PDCTBT leading to further improvements in OPV efficiency or hole mobility.

We performed molecular dynamics simulations of the low-molecular weight organic glass former ortho-terphenyl in bulk and freestanding films. The main motivation is to provide molecular insight into the confinement effect without explicit interfaces. Based on earlier models of ortho-terphenyl we developed an atomistic model for bulk simulations. The model reproduces literature data both from simulations and experiments starting from specific volume and diffusivity to mean square displacement and radial distribution functions. After characterizing the bulk model we form freestanding films by the elongation and expansion method. These films give us the opportunity to study the dynamical heterogeneity near the glass transition through in-plane mobility and reorientation dynamics. We finally compare the model in bulk and under confinement. We found qualitatively a lower glass transition temperature for the freestanding film compared to the bulk. [ABSTRACT FROM AUTHOR]

The dynamics of designed short polyisoprene (PI) chains in the melt is investigated on a wide temperature window using dielectric relaxation spectroscopy and pulsed field gradient nuclear magnetic resonance (NMR). At high temperatures, molecular dynamics (MD) simulations performed using two different models (an explicit atom model and a united atom one) capture very well the dynamic properties documented experimentally. Structures pre-equilibrated with end-bridging Monte Carlo are used as initial configurations for MD runs at different temperatures, providing predictions for the temperature dependence of the dynamics of this bulk PI. Local dynamics is unique, independently of the probe (dielectric relaxation, dynamic light scattering, nuclear magnetic resonance, neutron scattering), although mean correlation times are significantly affected, to different extents, by librations. Chain dynamics over the molecular weight and temperature range studied can be described well by the Rouse model, as shown by both experimental data and a normal mode analysis on simulation trajectories. Deviations from the Rouse model emerge for the high modes at short times; still, this model offers a rather accurate picture. © 2003 American Institute of Physics. [ABSTRACT FROM AUTHOR]

In cariology and research involving the interactions between restorative materials and dental hard tissues, it is important for small changes in tooth mineral content to be measurable. Currently, transverse microradiography (TMR) is the most accepted tool for the above purpose. Electron Probe Microanalysis (EPMA) can yield both qualitative identification of elements and quantitative compositional information. The purpose of this study was to compare the mineral distribution in well-defined artificial lesions, in dentin and enamel, by the use of both TMR and EPMA on the same sample. The good correlation between the two sets of data validates EPMA as a technique and helps in the interpretation of its results. The data from TMR analysis are expressed as vol% of mineral, while EPMA gives the levels of Ca and phosphate in wt%. The conversion between the two sets of data is complicated by the fact that local density is, as yet, unknown. [ABSTRACT FROM AUTHOR]