-
Graph Random Walk for Time-of-Flight Charge Mobilities
Authors:
Zhongquan Chen,
Pim van der Hoorn,
Björn Baumeier
Abstract:
We present a graph random walk (GRW) method for the study of charge transport properties of complex molecular materials in the time-of-flight regime. The molecules forming the material are represented by the vertices of a directed weighted graph, and the charge carriers are random walkers. The edge weights are rates for elementary jumping processes for a charge carrier to move along the edge and a…
▽ More
We present a graph random walk (GRW) method for the study of charge transport properties of complex molecular materials in the time-of-flight regime. The molecules forming the material are represented by the vertices of a directed weighted graph, and the charge carriers are random walkers. The edge weights are rates for elementary jumping processes for a charge carrier to move along the edge and are determined from a combination of the energies of the involved vertices and an interaction strength. Exclusions are built into the random walk to account for the Pauli exclusion principle. In time-of-flight experiments, charge carriers are injected into the material and the time until they reach a collecting electrode is recorded. In this setting, our GRW approach allows direct evaluation of the expected hitting time of the collecting nodes in the graph in terms of a typically sparse, linear system, thereby avoiding numerically cumbersome and potentially fluctuations-prone methods based on explicit time evolution from solutions of a high-dimensional system of coupled ordinary differential equations (the Master Equation) or from kinetic Monte Carlo (KMC). We validate the GRW approach by conducting numerical studies of charge dynamics of single and multiple carriers in diffusive and drift-diffusive regimes using a surrogate lattice model. The surrogate model allows varying types and strengths of energetic disorder from the reference baseline. Comparison with results from the Master Equation confirms the theoretical equivalence of both approaches also in numerical implementations. We further show that KMC results show substantial deviations due to inadequate sampling. All in all, we find that the GRW method provides a powerful alternative to the more commonly used methods without sampling issues and with the benefit of making use of sparse matrix methods.
△ Less
Submitted 24 May, 2024;
originally announced May 2024.
-
Glassy dynamics from generalized mode-coupling theory: existence and uniqueness of solutions for hierarchically coupled integro-differential equations
Authors:
Rutger A. Biezemans,
Simone Ciarella,
Onur Çaylak,
Björn Baumeier,
Liesbeth M. C. Janssen
Abstract:
Generalized mode-coupling theory (GMCT) is a first-principles-based and systematically correctable framework to predict the complex relaxation dynamics of glass-forming materials. The formal theory amounts to a hierarchy of infinitely many coupled integro-differential equations, which may be approximated using a suitable finite-order closure relation. Although previous studies have suggested that…
▽ More
Generalized mode-coupling theory (GMCT) is a first-principles-based and systematically correctable framework to predict the complex relaxation dynamics of glass-forming materials. The formal theory amounts to a hierarchy of infinitely many coupled integro-differential equations, which may be approximated using a suitable finite-order closure relation. Although previous studies have suggested that finite-order GMCT leads to well-defined solutions, and that the hierarchy converges as the closure level increases, no rigorous and general result in this direction is known. Here we unambiguously establish the existence and uniqueness of solutions to generic, schematic GMCT hierarchies that are closed at arbitrary order. We consider two types of commonly invoked closure approximations, namely mean-field and exponential closures. We also distinguish explicitly between overdamped and underdamped glassy dynamics, corresponding to hierarchies of first-order and second-order integro-differential equations, respectively. We find that truncated GMCT hierarchies closed under an exponential closure conform to previously developed mathematical theories, such that the existence of a unique solution can be readily inferred. Self-consistent mean-field closures, however, of which the well-known standard-MCT closure approximation is a special case, warrant additional arguments for mathematical rigour. We demonstrate that the existence of a priori bounds on the solution is sufficient to also prove that unique solutions exist for such self-consistent hierarchies. To complete our analysis, we present simple arguments to show that these a priori bounds must exist, motivated by the physical interpretation of the GMCT solutions as density correlation functions. Overall, our work contributes to the theoretical justification of GMCT for studies of the glass transition, placing GMCT on a firmer mathematical footing.
△ Less
Submitted 13 August, 2020; v1 submitted 8 June, 2020;
originally announced June 2020.
-
Jamming and force distribution in growing epithelial tissue
Authors:
Pranav Madhikar,
Jan Åström,
Björn Baumeier,
Mikko Karttunen
Abstract:
We investigate morphologies of proliferating cellular tissues using a newly developed numerical simulation model for mechanical cell division and migration in 2D. The model is applied to a bimodal mixture consisting of stiff cells with a low growth potential and soft cells with a high growth potential; cancer cells are typically considered to be softer than healthy cells. In an even mixture, the s…
▽ More
We investigate morphologies of proliferating cellular tissues using a newly developed numerical simulation model for mechanical cell division and migration in 2D. The model is applied to a bimodal mixture consisting of stiff cells with a low growth potential and soft cells with a high growth potential; cancer cells are typically considered to be softer than healthy cells. In an even mixture, the soft cells develop into a tissue matrix and the stiff cells into a dendrite-like network structure. When soft cells are placed inside a tissue consisting of stiff cells (to model cancer growth), the soft cells develop to a fast growing tumor-like structure that gradually evacuates the stiff cell matrix. The model also demonstrates 1) how soft cells orient themselves in the direction of the largest effective stiffness as predicted by the theory of Bischofs and Schwarz (Proc. Natl. Acad. Sci U.S.A., 100, 9274--9279 (2003) and 2) that the orientation and force generation continue a few cell rows behind the soft-stiff interface. With increasing inter-cell friction, tumor growth slows down and cell death occurs. The contact force distribution between cells is demonstrated to be highly sensitive to cell type mixtures and cell-cell interactions, which indicates that local mechanical forces can be useful as a regulator of tissue formation. The results shed new light on established experimental data.
△ Less
Submitted 4 May, 2020;
originally announced May 2020.
-
Quantitative Predictions of Photoelectron Spectra in Amorphous Molecular Solids from Multiscale Quasiparticle Embedding
Authors:
Gianluca Tirimbo,
Xander de Vries,
Christ H. L. Weijtens,
Peter A. Bobbert,
Reinder Coehoorn,
Bjoern Baumeier
Abstract:
We present a first-principles-based multiscale simulation framework for quantitative predictions of the high-energy part of the Ultraviolet Photoelectron Spectroscopy (UPS) spectra of amorphous molecular solids. The approach combines a deposition simulation, many-body Green's Function Theory, polarizable film-embedding, and multimode electron-vibrational coupling and provides a molecular-level vie…
▽ More
We present a first-principles-based multiscale simulation framework for quantitative predictions of the high-energy part of the Ultraviolet Photoelectron Spectroscopy (UPS) spectra of amorphous molecular solids. The approach combines a deposition simulation, many-body Green's Function Theory, polarizable film-embedding, and multimode electron-vibrational coupling and provides a molecular-level view on the interactions and processes giving rise to spectral features. This insight helps bridging the current gap between experimental UPS and theoretical models as accurate analyses are hampered by the energetic disorder, surface-sensitivity of the measurement and the complexity of excitation processes. In particular this is relevant for the unambiguous determination the highest occupied molecular orbital energy (HOMO) of organic semiconductors, a key quantity for tailoring and engineering new opto-electronic devices. We demonstrate the capabilities of the simulation approach studying the spectrum of two isomers of 2-methyl-9,10-bis(naphthalen-2-yl)anthracene (MADN) as archetypical materials showing a clearly separated HOMO peak in experiment. The agreement with experiment is excellent, suggesting that our approach provides a route for determining the HOMO energy with an accuracy better than 0.1eV.
△ Less
Submitted 15 December, 2019; v1 submitted 16 July, 2019;
originally announced July 2019.
-
Morphology of Proliferating Epithelial Cellular Tissue
Authors:
Pranav Madhikar,
Jan Åström,
Björn Baumeier,
Mikko Karttunen
Abstract:
We investigate morphologies of proliferating cellular tissue using a newly developed numerical simulation model for mechanical cell division. The model reproduces structures of simple multi-cellular organisms via simple rules for selective division and division plane orientation. The model is applied to a bimodal mixture of stiff cells with a low growth potential and soft cells with a high growth…
▽ More
We investigate morphologies of proliferating cellular tissue using a newly developed numerical simulation model for mechanical cell division. The model reproduces structures of simple multi-cellular organisms via simple rules for selective division and division plane orientation. The model is applied to a bimodal mixture of stiff cells with a low growth potential and soft cells with a high growth potential. In an even mixture, the soft cells develop into a tissue matrix and the stiff cells into a dendrite-like network structure. For soft cell inclusion in a stiff cellular matrix, the soft cells develop to a fast growing tumour like structure that gradually evacuates the stiff cell matrix. With increasing inter-cell friction, the tumour growth slows down and parts of it is driven to self-inflicted cell death.
△ Less
Submitted 19 November, 2018;
originally announced November 2018.
-
Improved general-purpose five-point model for water: TIP5P/2018
Authors:
Yuriy Khalak,
Bjoern Baumeier,
Mikko Karttunen
Abstract:
A new five point potential for liquid water, TIP5P/2018, is presented along with the techniques used to derive its charges from ab initio per-molecule electrostatic potentials in the liquid phase using the split charge equilibration (SQE) of Nistor et al. [J. Chem. Phys. 125, 094108 (2006)]. By taking the density and diffusion dependence on temperature as target properties, significant improvement…
▽ More
A new five point potential for liquid water, TIP5P/2018, is presented along with the techniques used to derive its charges from ab initio per-molecule electrostatic potentials in the liquid phase using the split charge equilibration (SQE) of Nistor et al. [J. Chem. Phys. 125, 094108 (2006)]. By taking the density and diffusion dependence on temperature as target properties, significant improvements to the behavior of isothermal compressibility were achieved along with improvements to other thermodynamic and rotational properties. While exhibiting a dipole moment close to ab initio values, TIP5P/2018 suffers from a too small quadrupole moment due to the charge assignment procedure and results in an overestimation of the dielectric constant.
△ Less
Submitted 18 October, 2018;
originally announced October 2018.
-
Getting excited: Challenges in quantum-classical studies of excitons in polymeric systems
Authors:
Behnaz Bagheri,
Björn Baumeier,
Mikko Karttunen
Abstract:
A combination of classical molecular dynamics (MM/MD) and quantum chemical calculations based on the density functional theory (DFT) was performed to describe conformational properties of diphenylethyne (DPE), methylated-DPE and poly para phenylene ethynylene (PPE). DFT calculations were employed to improve and develop force field parameters for MM/MD simulations. Many-body Green's functions theor…
▽ More
A combination of classical molecular dynamics (MM/MD) and quantum chemical calculations based on the density functional theory (DFT) was performed to describe conformational properties of diphenylethyne (DPE), methylated-DPE and poly para phenylene ethynylene (PPE). DFT calculations were employed to improve and develop force field parameters for MM/MD simulations. Many-body Green's functions theory within the GW approximation and the Bethe-Salpeter equation were utilized to describe excited states of the systems. Reliability of the excitation energies based on the MM/MD conformations was examined and compared to the excitation energies from DFT conformations. The results show an overall agreement between the optical excitations based on MM/MD conformations and DFT conformations. This allows for calculation of excitation energies based on MM/MD conformations.
△ Less
Submitted 9 May, 2016;
originally announced May 2016.
-
Solvent effects on optical excitations of poly para phenylene ethynylene studied by QM/MM simulations based on Many-Body Green's Functions Theory
Authors:
Behnaz Bagheri,
Mikko Karttunen,
Björn Baumeier
Abstract:
Electronic excitations in dilute solutions of poly para phenylene ethynylene (poly-PPE) are studied using a QM/MM approach combining many-body Green's functions theory within the $GW$ approximation and the Bethe-Salpeter equation with polarizable force field models. Oligomers up to a length of 7.5\,nm (10 repeat units) functionalized with nonyl side chains are solvated in toluene and water, respec…
▽ More
Electronic excitations in dilute solutions of poly para phenylene ethynylene (poly-PPE) are studied using a QM/MM approach combining many-body Green's functions theory within the $GW$ approximation and the Bethe-Salpeter equation with polarizable force field models. Oligomers up to a length of 7.5\,nm (10 repeat units) functionalized with nonyl side chains are solvated in toluene and water, respectively. After equilibration using atomistic molecular dynamics (MD), the system is partitioned into a quantum region (backbone) embedded into a classical (side chains and solvent) environment. Optical absorption properties are calculated solving the coupled QM/MM system self-consistently and special attention is paid to the effects of solvents. The model allows to differentiate the influence of oligomer conformation induced by the solvation from electronic effects related to local electric fields and polarization. It is found that the electronic environment contributions are negligible compared to the conformational dynamics of the conjugated PPE. An analysis of the electron-hole wave function reveals a sensitivity of energy and localization characteristics of the excited states to bends in the global conformation of the oligomer rather than to the relative of phenyl rings along the backbone.
△ Less
Submitted 2 May, 2016;
originally announced May 2016.