Angular overlap model parameterisation of Anderson''s superexchange theory. I. A quantification of Goodenough-Kanamori rules

A parameterisation of Anderson's exchange formulas on the basis of an extension of the angular overlap model (AOM) is proposed. Transfer integrals are expressed in terms of the metal-metal and metal-ligand bonding parameters, which can be estimated from independent spectroscopic studies, or calculated using solid state expressions. Analytical expressions for the transfer integrals between various d-orbitals, appropriate for cubic crystal lattices, comprising octahedra sharing common vertices and common edges serve as a quantification of the Goodenough-Kanamori rules. On the basis of the present parameterisation we give an explanation of the “exchange integral versus bond-distance” dependence. Some potential applications of the model are briefly discussed.     

Time-dependent quantum theory III. Model for vibrational relaxation in crystals

A one-dimensional model consisting of a “diatomic” spring attached on one side to a rigid wall and on the other side to a linear array of mass-springs is proposed as a model for the vibrational relaxation of small solute molecules in host lattices. A modification allowing a change in the equilibrium internuclear extension of the diatomic spring is also incorporated. The Hamiltonian divides naturally into pure diatomic, pure linear crystal, and the two mixed perturbation terms, one giving rise to stepwise vibrational cascade damping accompanied by phonon emission, and the other process, lattice relaxation, giving rise to phonon emission without any change of the quantum number of the diatomic spring. The cascade damping rate for a diatomic spring with a frequency less than the the maximum frequency of the linear crystal is calculated to second-order, and it is shown that the perturbation series converges in this range. An upper bound to the cascade damping rate for a diatomic spring with a frequency greater (i.e., 4.5 ×) than the cut-off frequency of the linear crystal is determined to be very small, λ ≦ 10⁴ sec ^?;1. The rate for the lattice relaxation process corresponds to a line-width λ = 6 cm ^?1 at ⁰K. An explanation for the thermal quenching of the low-temperature luminescence of SO₂ is based upon induced cascade-phonon emission.     

Derivation of quantum langevin equation from an explicit molecule-medium treatment in interaction picture

A quantum mechanical form of the Langevin equation is derived from an explicit consideration of the molecule-medium interaction, as advocated by Simons in 1978, and by using two identities in the interaction picture. This can be easily reduced to the classical regime, and further simplified to the macroscopic Langevin equation by considering the stochastic Langevin force autocorrelation function. One of the so-called Einstein relations appears as a byproduct. By following the methodology proposed by Simons, an exact expression for the momentum autocorrelation function is obtained. The latter can be used to calculate the zero-frequency macroscopic diffusion coefficient that is observed to satisfy the second Einstein relation. The formalism described above gives rise to the possibility of explicitly computing the transport characteristics such as friction constant and diffusion coefficient from the corresponding quantum statistical mechanical expressions. A discussion on the Langevin equation becomes complete only when the corresponding Fokker-Planck equation is obtained. Therefore, the probability of the evolution of states with a particular absolute magnitude of linear momentum from those of another momentum eigenvalue is quantum mechanically defined. This probability appears as a special average value of a projection operator and as a special projection operator correlation function. A classical identity is introduced that is shown to be valid also for the quantum mechanically defined probability function. By using this identity, the so-called Fokker-Planck equation for the evolution probability is easily established.     

Modeling vibrational dephasing and energy relaxation of intramolecular anharmonic modes for multidimensional infrared spectroscopies

Starting from a system-bath Hamiltonian in a molecular coordinate representation, we examine an applicability of a stochastic multilevel model for vibrational dephasing and energy relaxation in multidimensional infrared spectroscopy. We consider an intramolecular anharmonic mode nonlinearly coupled to a colored noise bath at finite temperature. The system-bath interaction is assumed linear plus square in the system coordinate, but linear in the bath coordinates. The square-linear system-bath interaction leads to dephasing due to the frequency fluctuation of system vibration, while the linear-linear interaction contributes to energy relaxation and a part of dephasing arises from anharmonicity. To clarify the role and origin of vibrational dephasing and energy relaxation in the stochastic model, the system part is then transformed into an energy eigenstate representation without using the rotating wave approximation. Two-dimensional (2D) infrared spectra are then calculated by solving a low-temperature corrected quantum Fokker-Planck (LTC-QFP) equation for a colored noise bath and by the stochastic theory. In motional narrowing regime, the spectra from the stochastic model are quite different from those from the LTC-QFP. In spectral diffusion regime, however, the 2D line shapes from the stochastic model resemble those from the LTC-QFP besides the blueshifts caused by the dissipation from the colored noise bath. The preconditions for validity of the stochastic theory for molecular vibrational motion are also discussed.     

Force constants and many-body electrostatic interactions in two-dimensional colloidal crystal

A model of a two-dimensional colloidal crystal with a hexagonal lattice, the electrostatic interactions in which are described by the nonlinear Poisson-Boltzmann equation, is considered. The calculation procedure for force constants of this crystal is treated in detail. Properties of system symmetry, which make it possible to significantly decrease the volume of calculations and to classify force constants, are analyzed. Numerical data for force constants of a crystal as functions of lattice parameters at different particle sizes are reported. A method that allows us to disclose the presence of many-body interactions in a system by the behavior of force constants at some interval of the values of lattice parameters is proposed. The application of this method to the system under consideration demonstrated that electrostatic interparticle interactions in the system cannot be reduced to simply a pair interaction of any kind; the introduction of many-body potentials is required for the adequate representation of the elastic properties of a crystal.     

Applications of the RISM equation to diatomic fluids: the liquids nitrogen,oxygen and bromine

The reference interaction site model (RISM) integral equation is used to study the equilibrium pair correlation for one-component liquids composed of homonuclear diatomic molecules. The integral equation is first tested by comparing the results obtained from it with those of computer simulation calculations. An internal consistency test is developed which seems to provide an a priori measure of the accuracy of the RISM equation. Then the theory is used to interpret the neutron scattering structure factor data taken on the liquids nitrogen and oxygen. For both of these liquids, a satisfactory explanation of the data is obtained by assuming that the short ranged repulsive forces between molecules are mimicked by a two-site hard core model. However, this simple model does not provide a satisfactory explanation of the data taken from liquid bromine. But, it is shown that a slightly more sophisticated model for the short-ranged repulsion between Br₂ molecules does provide an adequate explanation. With the molecular models determined by fitting the neutron scattering data, the RISM equation provides a method for determining the atom, atom to center-of-mass, and center-of-mass to center-of-mass intermolecular distribution functions in the diatomic liquids. From these functions, the local structures in the three liquids are analyzed. While orientational pair correlations are nearly negligible in both the liquids nitrogen and oxygen, these correlations are fairly substantial in liquid bromine. Furthermore, even when the orientation of a molecule does not greatly influence the orientation of its neighbors, it is found that the orientation of a molecule does have an important effect on the location of its neighbors. Thus, the coupling of translational and orientational coordinates is significant, even in the liquids nitrogen and oxygen.     

Crystal Memory near Discontinuous Triacylglycerol Phase Transitions: Models,Metastable Regimes,and Critical Points

It is proposed that “crystal memory”, observed in a discontinuous solid-liquid phase transition of saturated triacylglycerol (TAG) molecules, is due to the coexistence of solid TAG crystalline phases and a liquid TAG phase, in a superheated metastable regime. Such a coexistence has been detected. Solid crystals can act as heterogeneous nuclei onto which molecules can condense as the temperature is lowered. We outlined a mathematical model, with a single phase transition, that shows how the time-temperature observations can be explained, makes predictions, and relates them to recent experimental data. A modified Vogel-Fulcher-Tammann (VFT) equation is used to predict time-temperature relations for the observation of “crystal memory” and to show boundaries beyond which “crystal memory” is not observed. A plot of the lifetime of a metastable state versus temperature, using the modified VFT equation, agrees with recent time-temperature data. The model can be falsified through its predictions: the model possesses a critical point and we outline a procedure describing how it could be observed by changing the hydrocarbon chain length. We make predictions about how thermodynamic functions will change as the critical point is reached and as the system enters a crossover regime. The model predicts that the phenomenon of “crystal memory” will not be observed unless the system is cooled from a superheated metastable regime associated with a discontinuous phase transition.     

Electronic energy transfer in an impurity band

In this paper we explore the dynamics of triplet electronic energy transfer in an impurity band of a substitutionally-disordered material. The general characteristics of the solutions of the master equation in the strong scattering regime have been studied, leading to a generalized diffusion equation that contains a memory term, which results in a time-dependent diffusion coefficient. Approximate results for the average density of excitation were derived within the framework of the pair-approximation, resulting in explicit expressions for the time-resolved spectral diffusion, the mean-square displacement, as well as the time-dependent diffusion coefficient for electronic transfer processes induced by exchange interactions.     

Crystal lattice photochemistry often proceeds in discrete stages. Mechanistic and exploratory organic photochemistry

[formula: see text] While the solid-state photochemistry of crystals has been assumed to proceed to some critical point beyond which the crystal is destroyed, it has now been found that crystal lattice photochemistry often occurs in stages or regimes. In each regime different chemistry generally results. The reaction in different stages may be followed kinetically.     

Accordian-type Laser-Raman scattering by polymers

All the improvements of the independent-rod model of longitudinal accordian-type acoustic mode (LAM) oscillations have assumed that the oscillation energy is retained either on the isolated macromolecule oscillating in a vacuum or in a narrow cylinder containing the straight sections of the macromolecules in the crystal lattice and their straight continuations through the amorphous layers. According to such concepts, concentration of the oscillation energy in gauche defects or amorphous layers occurs, respectively, whenever the axial elastic modulus of the straight sections (crystal lattice) is very much larger than that of the kinked sections (amorphous layers). The effect is enhanced by low crystallinity. Actually such behavior has never been observed. To agree with experimental data the model has to be modified in such a manner that the oscillation amplitude in the amorphous layer steadily decreases with increasing distance from the boundary between the two phases. The necessary large damping of the LAM oscillation in the kinked sections results from true damping in the viscoelastic amorphous component and energy transfer to adjacent chains which turns out to be just as easy as energy conduction along the kinked chain. Such a transfer is equivalent to radiation of the oscillation energy in all directions in the kinked phase. As a consequence of damping, the coupling of chains in adjacent crystals becomes so small that it may be completely neglected. Such a model explains in a satisfactory manner the observed accordion Laser-Raman spectra of the semicrystalline polymers and the infrared absorption of paraffins in the liquid state.     

Comparison between theoretical values and simulation results of viscosity for the dissipative particle dynamics method

In order to investigate the validity of the dissipative particle dynamics method, which is a mesoscopic simulation technique, we have derived an expression for viscosity from the equation of motion of dissipative particles. In the concrete, we have shown the Fokker-Planck equation in phase space, and macroscopic conservation equations such as the equation of continuity and the equation of momentum conservation. The basic equations of the single-particle and pair distribution functions have been derived using the Fokker-Planck equation. The solutions of these distribution functions have approximately been solved by the perturbation method under the assumption of molecular chaos. The expressions of the viscosity due to momentum and dissipative forces have been obtained using the approximate solutions of the distribution functions. Also, we have conducted nonequilibrium dynamics simulations to investigate the influence of the parameters, which have appeared in defining the equation of motion in the dissipative particle dynamics method. The theoretical values of the viscosity due to dissipative forces in the Hoogerbrugge-Koelman theory are in good agreement with the simulation results obtained by the nonequilibrium dynamics method, except in the range of small number densities. There are restriction conditions for taking appropriate values of the number density, number of particles, time interval, shear rate, etc., to obtain physically reasonable results by means of dissipative particle dynamics simulations.     

How accurate are the nonlinear chemical Fokker-Planck and chemical Langevin equations?

The chemical Fokker-Planck equation and the corresponding chemical Langevin equation are commonly used approximations of the chemical master equation. These equations are derived from an uncontrolled, second-order truncation of the Kramers-Moyal expansion of the chemical master equation and hence their accuracy remains to be clarified. We use the system-size expansion to show that chemical Fokker-Planck estimates of the mean concentrations and of the variance of the concentration fluctuations about the mean are accurate to order Ω(-3∕2) for reaction systems which do not obey detailed balance and at least accurate to order Ω(-2) for systems obeying detailed balance, where Ω is the characteristic size of the system. Hence, the chemical Fokker-Planck equation turns out to be more accurate than the linear-noise approximation of the chemical master equation (the linear Fokker-Planck equation) which leads to mean concentration estimates accurate to order Ω(-1∕2) and variance estimates accurate to order Ω(-3∕2). This higher accuracy is particularly conspicuous for chemical systems realized in small volumes such as biochemical reactions inside cells. A formula is also obtained for the approximate size of the relative errors in the concentration and variance predictions of the chemical Fokker-Planck equation, where the relative error is defined as the difference between the predictions of the chemical Fokker-Planck equation and the master equation divided by the prediction of the master equation. For dimerization and enzyme-catalyzed reactions, the errors are typically less than few percent even when the steady-state is characterized by merely few tens of molecules.     

Non-Markovian theory of open systems in classical limit

A fully classical limit of the recently published quantum-classical approximation [A. A. Neufeld, J. Chem. Phys. 119, 2488 (2003)] is obtained and analyzed. The resulting kinetic equations are capable of describing the evolution of an open system on the entire time axis, including the short-time non-Markovian stage, and are valid beyond linear response regime. We have shown, that proceeding to the classical mechanics limit we restrict the class of allowed correlations between an open system and a canonical bath, so that the initial conditions and the relaxation operator has to be appropriately modified (projected). Disregard of the projection may lead to unphysical behavior, since mechanism of the decay of some correlations is essentially of quantum-mechanical nature, and is not correctly described by classical mechanics. The projection (quantum correction to the kinetics) is particularly important for the non-Markovian regime of relaxation towards canonical equilibrium. The conformity of the developed method to the conventional approaches is demonstrated using a model of Brownian motion (heavy particle in the bath of light ones), for which the obtained non-Markovian equations are reduced to the standard Fokker-Planck equation in phase space.     

Kinetic theory of binary nucleation based on a first passage time analysis

The binary classical nucleation theory (BCNT) is based on the Gibbsian thermodynamics and applies the macroscopic concept of surface tension to nanosize clusters. This leads to severe inconsistencies and large discrepancies between theoretical predictions and experimental results regarding the nucleation rate. We present an alternative approach to the kinetics of binary nucleation which avoids the use of classical thermodynamics for clusters. The new approach is an extension to binary mixtures of the kinetic theory previously developed by Narsimhan and Ruckenstein and Ruckenstein and Nowakowski [J. Colloid Interface Sci. 128, 549 (1989); 137, 583 (1990)] for unary nucleation which is based on molecular interactions and in which the rate of emission of molecules from a cluster is determined via a mean first passage time analysis. This time is calculated by solving the single-molecule master equation for the probability distribution of a "surface" molecule moving in a potential field created by the cluster. The starting master equation is a Fokker-Planck equation for the probability distribution of a surface molecule with respect to its phase coordinates. Owing to the hierarchy of characteristic time scales in the evolution of the molecule, this equation can be reduced to the Smoluchowski equation for the distribution function involving only the spatial coordinates. The new theory is combined with density functional theory methods to determine the density profiles. This is essential for nucleation in binary systems particularly when one of the components is surface active. Knowing these profiles, one can determine the potential fields created by the cluster, its rate of emission of molecules, and the nucleation rate more accurately than by using the uniform density approximation. The new theory is illustrated by numerical calculations for a model binary mixture of Lennard-Jones monomers and rigidly bonded dimers of Lennard-Jones atoms. The amphiphilic character of the dimer component (i.e., its surface activity) is induced by the asymmetry in the interaction between a monomer and the two different sites of a dimer. The inconsistencies of the BCNT are avoided in the new theory.     

Recent developments in the kinetic theory of nucleation

A review of recent progress in the kinetics of nucleation is presented. In the conventional approach to the kinetic theory of nucleation, it is necessary to know the free energy of formation of a new-phase particle as a function of its independent variables at least for near-critical particles. Thus the conventional kinetic theory of nucleation is based on the thermodynamics of the process. The thermodynamics of nucleation can be examined by using various approaches, such as the capillarity approximation, density functional theory, and molecular simulation, each of which has its own advantages and drawbacks. Relatively recently a new approach to the kinetics of nucleation was proposed [Ruckenstein E, Nowakowski B. J Colloid Interface Sci 1990;137:583; Nowakowski B, Ruckenstein E. J Chem Phys 1991;94:8487], which is based on molecular interactions and does not employ the traditional thermodynamics, thus avoiding such a controversial notion as the surface tension of tiny clusters involved in nucleation. In the new kinetic theory the rate of emission of molecules by a new-phase particle is determined with the help of a mean first passage time analysis. This time is calculated by solving the single-molecule master equation for the probability distribution function of a surface layer molecule moving in a potential field created by the rest of the cluster. The new theory was developed for both liquid-to-solid and vapor-to-liquid phase transitions. In the former case the single-molecule master equation is the Fokker-Planck equation in the phase space which can be reduced to the Smoluchowski equation owing to the hierarchy of characteristic time scales. In the latter case, the starting master equation is a Fokker-Planck equation for the probability distribution function of a surface layer molecule with respect to both its energy and phase coordinates. Unlike the case of liquid-to-solid nucleation, this Fokker-Planck equation cannot be reduced to the Smoluchowski equation, but the hierarchy of time scales does allow one to reduce it to the Fokker-Plank equation in the energy space. The new theory provides an equation for the critical radius of a new-phase particle which in the limit of large clusters (low supersaturations) yields the Kelvin equation and hence an expression for the macroscopic surface tension. The theory was illustrated with numerical calculations for a molecular pair interaction potential combining the dispersive attraction with the hard-sphere repulsion. The results for the liquid-to-solid nucleation clearly show that at given supersaturation the nucleation rate depends on the cluster structure (for three cluster structures considered-amorphous, fcc, and icosahedral). For both the liquid-to-solid and vapor-to-liquid nucleation, the predictions of the theory are consistent with the results of classical nucleation theory (CNT) in the limit of large critical clusters (low supersaturations). For small critical clusters the new theory provides higher nucleation rates than CNT. This can be accounted for by the fact that CNT uses the macroscopic interfacial tension which presumably overpredicts the surface tension of small clusters, and hence underpredicts nucleation rates.     

Geometry optimization of periodic systems using internal coordinates

An algorithm is proposed for the structural optimization of periodic systems in internal (chemical) coordinates. Internal coordinates may include in addition to the usual bond lengths, bond angles, out-of-plane and dihedral angles, various "lattice internal coordinates" such as cell edge lengths, cell angles, cell volume, etc. The coordinate transformations between Cartesian (or fractional) and internal coordinates are performed by a generalized Wilson B-matrix, which in contrast to the previous formulation by Kudin et al. [J. Chem. Phys. 114, 2919 (2001)] includes the explicit dependence of the lattice parameters on the positions of all unit cell atoms. The performance of the method, including constrained optimizations, is demonstrated on several examples, such as layered and microporous materials (gibbsite and chabazite) as well as the urea molecular crystal. The calculations used energies and forces from the ab initio density functional theory plane wave method in the projector-augmented wave formalism.     

Thermal Destiny of an Ionic Crystal

Recently it has been shown that defect‐defect interactions leading to anomalous rise in defect concentrations and eventually to phase transitions into a disordered state can be successfully described by a quasi defect lattice model [1]. Using this approach the fact is highlighted that the evolution of an ionic crystal from the perfect state at T = 0 K to a weakly defective ideal crystal and eventually via passing a pre‐melting regime to a superionic crystal (or melt) is unavoidable (given sufficient stability) and can be roughly phenomenologically calculated by using a few basic parameters only. Semi‐empirical relations (e. g. Tammann rule) can be quantified in this way.     

Thermally activated escape rate for the Brownian motion of a fixed axis rotator in an asymmetrical double-well potential for all values of the dissipation

The Kramers theory of the escape rate of a Brownian particle from a potential well as extended by Mel'nikov and Meshkov is used to evaluate the relaxation times and the dynamic susceptibility for the rotational Brownian motion of fixed axis rotators in an asymmetric double-well potential. An expression for the escape rate valid for all values of the dissipation including the very low damping (VLD), very high damping (VHD), and crossover regimes is derived. It is shown that this expression provides a good asymptotic estimate of the inverse of the smallest nonvanishing eigenvalue lambda(1) of the underlying Fokker-Planck operator calculated by using the matrix-continued fraction method. For low barriers, where the Mel'nikov and Meshkov approach is not applicable, analytic equations for the correlation time tau( parallel) of the longitudinal dipole correlation function in the VLD and VHD limits are derived and a simple extrapolating equation valid for all values of the damping is proposed.