Charge-Transfer Mobility Parameters in Photoelectronic Devices: The Advanced Miller-Abrahams Computation.

The local hopping step of the electron transfer (ET) reaction is investigated for a real organic material composed of molecules M (N,N'-di(1-naphthyl)-N,N'-diphenyl-(1,1'-biphenyl)-4,4'-diamine). This material is implemented in light-emitting photoelectronic devices. The conductivity effect is simulated and calculated at a molecular level. We have studied the ET mechanism alternative to that suggested by the usually employed Marcus-like polaron model. The ion-molecular binary complex M(+)M (for hole transfer) is considered as a reaction center. The reaction dynamics is carried through the low-frequency intermolecular vibration coordinate connecting its fragments (the promotion mode). Its coupling to the acoustic phonon bath serves for a dissipation of the reaction energy misfit. The high-frequency intramolecular vibrations (the reorganization modes) modulate the reaction kinetics via Franck-Condon factors induced by their polarization. The ET rate constants are evaluated in terms of the computational algorithm described earlier (Basilevsky, M. V.; et al. J. Chem. Phys. 2013 139, 234102). Standard quantum-chemical and molecular dynamical techniques are used for a calculation of all necessary parameters of this model. The macroscopic charge-carrier mobility of the material is estimated by properly averaging the rate constants over the total simulation cell.

Golden rule kinetics of transfer reactions in condensed phase: the microscopic model of electron transfer reactions in disordered solid matrices.

The algorithm for a theoretical calculation of transfer reaction rates for light quantum particles (i.e., the electron and H-atom transfers) in non-polar solid matrices is formulated and justified. The mechanism postulated involves a local mode (an either intra- or inter-molecular one) serving as a mediator which accomplishes the energy exchange between the reacting high-frequency quantum mode and the phonon modes belonging to the environment. This approach uses as a background the Fermi golden rule beyond the usually applied spin-boson approximation. The dynamical treatment rests on the one-dimensional version of the standard quantum relaxation equation for the reduced density matrix, which describes the frequency fluctuation spectrum for the local mode under consideration. The temperature dependence of a reaction rate is controlled by the dimensionless parameter ξ0 = ℏω0/k(B)T where ω0 is the frequency of the local mode and T is the temperature. The realization of the computational scheme is different for the high/intermediate (ξ0 < 1 - 3) and for low (ξ0 ≫ 1) temperature ranges. For the first (quasi-classical) kinetic regime, the Redfield approximation to the solution of the relaxation equation proved to be sufficient and efficient in practical applications. The study of the essentially quantum-mechanical low-temperature kinetic regime in its asymptotic limit requires the implementation of the exact relaxation equation. The coherent mechanism providing a non-vanishing reaction rate has been revealed when T → 0. An accurate computational methodology for the cross-over kinetic regime needs a further elaboration. The original model of the hopping mechanism for electronic conduction in photosensitive organic materials is considered, based on the above techniques. The electron transfer (ET) in active centers of such systems proceeds via local intra- and intermolecular modes. The active modes, as a rule, operate beyond the kinetic regimes, which are usually postulated in the existing theories of the ET. Our alternative dynamic ET model for local modes immersed in the continuum harmonic medium is formulated for both classical and quantum regimes, and accounts explicitly for the mode∕medium interaction. The kinetics of the energy exchange between the local ET subsystem and the surrounding environment essentially determine the total ET rate. The efficient computer code for rate computations is elaborated on. The computations are available for a wide range of system parameters, such as the temperature, external field, local mode frequency, and characteristics of mode/medium interaction. The relation of the present approach to the Marcus ET theory and to the quantum-statistical reaction rate theory [V. G. Levich and R. R. Dogonadze, Dokl. Akad. Nauk SSSR, Ser. Fiz. Khim. 124, 213 (1959); J. Ulstrup, Charge Transfer in Condensed Media (Springer, Berlin, 1979); M. Bixon and J. Jortner, Adv. Chem. Phys. 106, 35 (1999)] underlying it is discussed and illustrated by the results of computations for practically important target systems.

Association constants and distribution functions for ion pairs in binary solvent mixtures: application to a cyanine dye system.

The computations of the association constants K(ass) were performed at the microscopic level for the ion pair Cy(+)I(-) composed of the complex cyanine dye cation Cy(+) coupled to the negative iodine counterion. The wide array of K(ass) values is arranged by a variation of the composition of the binary solvent mixtures toluene/dimethylsulfoxide with the accompanying change of the solvent polarity. The potentials of mean force (PMFs) are calculated for a set of interionic separations R in the Cy(+)I(-) by a methodology which combines the quantum-chemical techniques for the treatment of the electronic structure of the Cy(+)I(-) system with the recent dielectric continuum approach which accounts for the solvation effects. For a given solute/solvent system the probability function P(R), which describes the distribution of interionic separations, is constructed in terms of the PMFs and implemented for the evaluation of the K(ass).

Decay kinetics of the excited S1 state of the cyanine dye Cy(+)I(-) (thiacarbocyanine iodide): the computation of quantum yields for different pathways.

This work explains the unordinary solvent effect which was observed in the photochemical decay kinetics for the cyanine dye thiacarbocyanine iodide (Cy(+)I(-)) in binary solvent mixtures toluene/dimethylsulfoxide. The interpretation is formulated in terms of the probability density F(R) describing the distribution of interionic distances R in the ion pair Cy(+)I(-) and depending on the solvent composition. The proper normalization of this distribution is expressed via the degree of association α for the ion pair in a given solvent mixture. The α values are, in turn, extracted by means of the mass action law from the ionic association constants computed in a separate publication. The detailed kinetic scheme includes the empirical parametrization of the R-dependent kinetic constants for different decay channels. The multiparameter fitting procedure represents, with the reasonable parameter values, the dependence of the observed quantum yields on the solvent composition.

Specific features of the dielectric continuum solvation model with a position-dependent permittivity function.

We consider a modified formulation for the recently developed new approach in the continuum solvation theory (Basilevsky, M. V.; Grigoriev, F. V.; Nikitina, E. A.; Leszczynski, J. J. Phys. Chem. B 2010, 114, 2457) which is based on the exact solution of the electrostatic Poisson equation with the space-dependent dielectric permittivity. Its present modification ensures the property curl E = 0 for the electric strength field E inherent to this solution, which is the obligatory condition imposed by Maxwell equations. The illustrative computation is made for the model system of the point dipole immersed in a spherical cavity of excluded volume.

Cavitation free energy for organic molecules having various sizes and shapes.

Cavitation free energy DeltaG(cav), corresponding to the formation of an excluded volume cavity in water, is calculated for a large set of organic molecules employing the thermodynamic integration procedure, which is realized as the original two-step algorithm for growing the interaction potential between the hard cavity wall and the water molecules. A large variety of solute systems is considered. Their characteristic radii change in the range 3-7 A; spherical cavities with radii 3-6 A are also studied. The interaction between water molecules is described by the four-site nonpolarizable TIP4P model. The diversity of the trial molecular set is provided by using a specially formulated nonspherical criterion classifying the cavity shapes according to their deviation from a sphere. Molecular objects were partly taken from the data base NCI Diversity with the aid of this criterion. The so-computed free energies are approximated by the linear volume dependence DeltaG(cav)V = XiV, where V is the cavity volume. This relation works fairly well until the cavity size becomes very large (the effective radius larger than 7 A). The volume dependence valid for solutes of arbitrary shapes can be included in a calculation of the nonpolar free energy component as required in the implicit water model.

Non-Markovian modification of the golden rule rate expression.

The reformulation of the standard golden rule approach considered in this paper for treating reactive tunneling reduces the computation of the reaction rate to a derivation of band shapes for energy levels of reactant and product states. This treatment is based on the assumption that the medium environment is actively involved as a partner in the energy exchange with the reactive subsystem but its reorganization effect is negligible. Starting from the quantum relaxation equation for the density matrix, the required band shapes are represented in terms of the spectral density function, exhibiting the continuum spectrum inherent to the interaction between the reactants and the medium in the total reactive system. The simplest Lorentzian spectral bands, obtained under Redfield approximation, proved to be unsatisfactory because they produced a divergent rate expression at low temperature. The problem is resolved by invoking a refined spectral band shape, which behaves as Lorentzian one at the band center but decays exponentially at its tails. The corresponding closed non-Markovian rate expression is derived and investigated taking as an example the photochemical H-transfer reaction between fluorene and acridine proceeding in the fluorene molecular crystal. The kinetics in this reactive system was thoroughly studied experimentally in a wide temperature range [B. Prass et al., Ber. Bunsenges. Phys. Chem. 102, 498 (1998)].

The model of level broadening in condensed phase.

We study a model of non-Markovian kinetics for a harmonic oscillator embedded in a harmonic heat bath. We present a new scheme for approximately solving the quantum relaxation equation for the density matrix to find a distribution of level populations. It is found to be an extended Lorentzian with the width depending on the energy. A more convenient non-Markovian distribution called square root Fourier distribution that was implemented in the preceding paper [M. V. Basilevsky et al., J. Chem. Phys. 125, 194513 (2006)] is closely related to this extended Lorentzian model. Both distributions decay exponentially far away from their centers and reproduce well standard Lorentzian widths in the vicinity of the central region. A conventional Lorentzian model with such widths results when the Redfield approximation is applied in the frame of the present procedure.

Molecular simulations of outersphere reorganization energies in polar and quadrupolar solvents. The case of intramolecular electron and hole transfer.

Outersphere reorganization energies (lambda) for intramolecular electron and hole transfer are studied in anion- and cation-radical forms of complex organic substrates (p-phenylphenyl-spacer-naphthyl) in polar (water, 1,2-dichloroethane, tetrahydrofuran) and quadrupolar (supercritical CO2) solvents. Structure and charge distributions of solute molecules are obtained at the HF/6-31G(d,p) level. Standard Lennard-Jones parameters for solutes and the nonpolarizable simple site-based models of solvents are used in molecular dynamics (MD) simulations. Calculation of lambda is done by means of the original procedure, which treats electrostatic polarization of a solvent in terms of a usual nonpolarizable MD scheme supplemented by scaling of reorganization energies at the final stage. This approach provides a physically relevant background for separating inertial and inertialless polarization responses by means of a single parameter epsilon(infinity), optical dielectric permittivity of the solvent. Absolute lambda values for hole transfer in 1,2-dichloroethane agree with results of previous computations in terms of the different technique (MD/FRCM, Leontyev, I. V.; et al. Chem. Phys. 2005, 319, 4). Computed lambda values for electron transfer in tetrahydrofuran are larger than the experimental values by ca. 2.5 kcal/mol; for the case of hole transfer in 1,2-dichloroethane the discrepancy is of similar magnitude provided the experimental data are properly corrected. The MD approach gives nonzero lambda values for charge-transfer reaction in supercritical CO2, being able to provide a uniform treatment of nonequilibrium solvation phenomena in both quadrupolar and polar solvents.

Excluded volume effect for large and small solutes in water.

The cavitation effect, i.e., the process of the creation of a void of excluded volume in bulk solvent (a cavity), is considered. The cavitation free energy is treated in terms of the information theory (IT) approach [Hummer, G.; Garde, S.; Garcia, A. E.; Paulaitis, M. E.; Pratt, L. R. J. Phys. Chem. B 1998, 102, 10469]. The binomial cell model suggested earlier is applied as the IT default distribution p(m) for the number m of solute (water) particles occupying a cavity of given size and shape. In the present work, this model is extended to cover the entire range of cavity size between small ordinary molecular solutes and bulky biomolecular structures. The resulting distribution consists of two binomial peaks responsible for producing the free energy contributions, which are proportional respectively to the volume and to the surface area of a cavity. The surface peak dominates in the large cavity limit, when the two peaks are well separated. The volume effects become decisive in the opposite limit of small cavities, when the two peaks reduce to a single-peak distribution as considered in our earlier work. With a proper interpolation procedure connecting these two regimes, the MC simulation results for model spherical solutes with radii increasing up to R = 10 A [Huang, D. H.; Geissler, P. L.; Chandler, D. J. Phys. Chem. B 2001, 105, 6704] are well reproduced. The large cavity limit conforms to macroscopic properties of bulk water solvent, such as surface tension, isothermal compressibility and Tolman length. The computations are extended to include nonspherical solutes (hydrocarbons C1-C6).

Momentum representation of the solute/bath interaction in the dynamic theory of chemical processes in condensed phase.

For the system consisting of the chemically reactive solute immersed in the oscillator bath, we consider an approach based on the solute/medium interaction expressed in terms of momenta rather than coordinates. In the adiabatic representation the medium reorganization effects are suppressed, being hidden in the solute renormalized potential and new spectral density function. The advantage proposed by the bilinear interaction in momentum representation is its spatial uniformity important for approximate dynamical treatments. The procedure of explicit transforming a standard spectral density (coordinate representation of interaction) into the spectral density in adiabatic representation (momentum representation of interaction) is the main new result of the present study. Illustrative calculations for several types of spectral functions are performed. Special discussion is devoted to clarifying the nature of the slow diffusion coordinate, to which the present approach is mainly addressed.