Condensed-phase relaxation of multilevel quantum systems. II. Comparison of path integral calculations and second-order relaxation theory for a nondegenerate three-level system

An exactly solvable model of multisite condensed-phase vibrational relaxation was studied in Paper I (Peter, S.; Evans, D. G.; Coalson, R. D. J. Phys. Chem. B 2006, 110, 18758.), where it was shown that long-time steady-state site populations of a degenerate N-level system are not equal (hence, they are non-Boltzmann) and depend on the initial preparation of the system and the number of sites that it comprises. Here we consider a generalization of the model to the case of a nondegenerate three-level system coupled to a high-dimensional bath: such a model system has direct relevance to a large class of donor-bridge-acceptor electron transfer processes. Because the quantum dynamics of this system cannot be computed analytically, we compare numerically exact path integral calculations to the predictions of second-order time-local relaxation theory. For modest system-bath coupling strengths, the two sets of results are in excellent agreement. They show that non-Boltzmann long-time steady-state site populations are obtained when the level splitting is small but nonzero, whereas at larger values of the system bias (asymmetry) these populations become Boltzmann distributed.     

An exactly solvable model for the Fermi contact interaction

A model for the Fermi contact interaction is proposed in which the nuclear moment is represented as a magnetized spherical shell of radius r ₀. For a hydrogen-like system thus perturbed, the Schrödinger equation is solvable without perturbation theory by use of the Coulomb Green's function. Approximation formulas are derived in terms of a quantum defect in the Coulombic energy formula. It is shown that the usual Fermi potential cannot be applied beyond first-order perturbation theory.     

An exactly solvable Ogston model of gel electrophoresis. VI. Towards a theory for macromolecules

In this article, we present a generalized version of our lattice model of low-field gel electrophoresis that allows us to treat the case of macromolecules such as short linear or circular oligomers and semi-flexible rods. We show that free-solution electrophoresis problems can be seen as random walks in the conformational space of the analyte. For sufficiently small molecules, our mathematical approach provides exact mobilities. In a quenched gel-like environment, however, both conformational and positional degrees of freedom must be used, but exact solutions can also be obtained. As an example, we then investigate several two-dimensional model gels, as well as a simple channel system where we see evidence of entropic effects that cannot be captured by the traditional Ogston concept of free volume.     

An exactly solvable Ogston model of gel electrophoresis: X. Application to high-field separation techniques

Recently, we generalized our lattice model of gel electrophoresis to study the net velocity of particles being pulled by a high-intensity electric field through an arbitrary distribution of immobile obstacles (Gauthier, M. G., Slater, G. W., J. Chem. Phys. 2002, 117, 6745-6756). In this article, we show how the high-field version of our model can be used to compare the velocity of particles with different electric charges and/or physical sizes. We then investigate specific two-dimensional distributions of obstacles that can be used to separate particles, e.g., in a microfluidic device. More precisely, we compare the velocity of differently charged or sized analytes in sieving, trapping and deflecting systems to model various electrophoretic separation techniques. In particular, we study the nonlinear effects present in ratchet systems and how they can be combined with time-asymmetric pulsed fields to provide new modes of separation.     

An exactly solvable model for a ternary solution with three-body interactions and orientationally dependent bonding

A model is considered in which the bonds of a honeycomb lattice are covered by rodlike molecules of types AA, BB, and AB. Neighboring molecular ends have three-body and orientation-dependent interactions. The model is shown to be equivalent to a spin-1/2 Ising model on the same lattice with a field, but with only pairwise interactions. Symmetric and asymmetric coexistence surfaces for the separation into an AA-rich and a BB-rich phase are calculated exactly.     

An exactly solvable Ogston model of gel electrophoresis. V. Attractive gel-analyte interactions and their effects on the Ferguson plot

We examine the effect of attractive analyte-gel interactions within the framework of our recently developed lattice model of gel electrophoresis. We show that it is possible to take into account such interactions and still calculate exact mobilities for various analytes and gel structures. Our study then focuses on two main issues: (i) the effect of these interactions on the separation efficiency of the Ogston regime; and (ii) the presence of inflection points (changes of curvature) in Ferguson plots. We establish some general principles, and we describe the results for selected two- and three-dimensional model systems. Numerous practical problems, such as chiral separations and affinity electrophoresis, can be treated using this approach.     

A class of exactly solvable matrix models

Some exactly solvable matrix models are discussed. Possible applications to problems in physical chemistry are pointed out, in particular the Hückel problem, the problem of torsional vibrations of polyatomic molecules, and of vibrations of finite polymer chains. This revised version was published online in July 2006 with corrections to the Cover Date.     

Does really Born–Oppenheimer approximation break down in charge transfer processes? An exactly solvable model

Effects of deviation from the Born–Oppenheimer approximation (BOA) on the non-adiabatic transition probability for the transfer of a quantum particle in condensed media are studied within an exactly solvable model. The particle and the medium are modeled by a set of harmonic oscillators. The dynamic interaction of the particle with a single local mode is treated explicitly without the use of BOA. Two particular situations (symmetric and non-symmetric systems) are considered. It is shown that the difference between the exact solution and the true BOA is negligibly small at realistic parameters of the model. However, the exact results differ considerably from those of the crude Condon approximation (CCA) which is usually considered in the literature as a reference point for BOA (Marcus–Hush–Dogonadze formula). It is shown that the exact rate constant can be smaller (symmetric system) or larger (non-symmetric one) than that obtained in CCA. The non-Condon effects are also studied.     

An exactly solvable Ogston model of gel electrophoresis: VIII. Nonconducting gel fibers, curved field lines, and the Nernst-Einstein relation.

In this article, we examine the low-field electrophoretic migration of infinitely small analytes in dilute sieving media made of nonconducting gel fibers. Using an Ogston obstruction model, we show that the electrophoretic mobility is not affected by the presence of curved field lines. In other words, the Nernst-Einstein relation between the mobility and the diffusion coefficient is valid regardless of the electrical properties of the gel fibers. Although this finding may greatly simplify the development of obstruction models of electrophoretic sieving, it also represents a critical test for any analytical or computational approach.     

On the applicability of the Born-Oppenheimer approximation in a theory of nonadiabatic charge transfer reactions: An exactly solvable model

The question concerning the applicability of the Born-Oppenheimer approximation (BOA) for the calculation of the transition probability for a nonadiabatic process of charge transfer in a polar environment with allowance made for temperature effects is investigated theoretically. Considered is the transfer of a quantum particle (proton) that interacts with a local vibration mode in a model of bound harmonic oscillators. The model admits an exact solution for wavefunctions of the initial and final states. A calculation shows that BOA is applicable even for very large distances of the proton transfer. At the same time, the exact result and the BOA result severely differ from a probability calculated in a crude Condon approximation. It is demonstrated that the non-Condon effects are in a general case temperature-dependent and may substantially influence the calculated values of the transition probability.     

An exactly soluble base problem for atomic systems

An exactly soluble base problem for atomic systems is presented which approximates the atomic Hamiltonian as a sum of identical one‐electron operators. The eigenfunctions of the one‐electron operator consist of a radial function multiplied by a spherical harmonic. Energies for many‐electron atoms are found by summing the one‐electron energy eigenvalues according to the Pauli principle. These energies are rigorous lower bounds to the exact energies. This revised version was published online in July 2006 with corrections to the Cover Date.     

Length dependence of the optical transition energies of the exactly solvable Hubbard model for linear polyenes

Some accurate results on the length dependence of the excitation energies from the ground state to ionic excited states in the Hubbard model of linear polyenes are obtained based on the method of Lieb and Wu. To this end, it is first shown that singly ionic excited states with “plus” alternancy symmetry in the Hubbard model are described by the wave functions in which the two electron operator [∑(?)ⁿCC] is acted on (N ? 2)-electron covalent eigenstates. Then by solving the Lieb-Wu equations the exact excitation energies of the lowest ionic state, which corresponds to the E state in this model, are calculated for systems with up to 50 electrons. The result, together with a correction for the end effect, indicates that the excitation energies do not decrease as 1/N but converge to the limiting value more rapidly when the number of electrons N becomes large.     

The dynamics of equally spaced multilevel model systems

Coherent dynamics of multilevel systems is of interest for the elucidation of energy disposal and multiphoton chemistry of polyatomic molecules. We compare the analytical solutions for coherent multiphoton excitation of two sparse multilevel model systems: (1) The spin-J system consisting of N = (2J + 1) equally spaced levels with the radiative coupling being determined by the xth component of the angular momentum. (2) The equally coupled N level system. The transition amplitudes for system (1) are always periodic, which is in contrast to the nonperiodic behavior exhibited by system (2) for N > 3.     

Markovian approximation in the relaxation of open quantum systems

In this paper, we examine the validity of the Markovian approximation and the slippage scheme used to incorporate short time transient memory effects in the Markovian master equations (Redfield equations). We argue that for a bath described by a spectral function, J(omega), that is dense and smoothly spread out over the range omega(d), a time scale of tau(b) approximately 1/omega(d) exists; for times of t > tau(b), the Markovian approximation is applicable. In addition, if J(omega) decays to zero reasonably fast in both the omega --> 0 and omega --> infinity limits, then the bath relaxation time, tau(b), is determined by the width of the spectral function and is weakly dependent on the temperature of the bath. On the basis of this criterion of tau(b), a scheme to incorporate transient memory effects in the Markovian master equation is suggested. Instead of using slipped initial conditions, we propose a concatenation scheme that uses the second-order perturbation theory for short time dynamics and the Markovian master equation at long times. Application of this concatenation scheme to the spin-boson model shows that it reproduces the reduced dynamics obtained from the non-Markovian master equation for all parameters studied, while the simple slippage scheme breaks down at high temperatures.     

Nuclear magnetic resonance restricted diffusion between parallel planes in a cosine magnetic field: an exactly solvable model

We propose a theoretical and numerical analysis of restricted diffusion between parallel planes in a cosine magnetic field. The specific choice of this spatial profile as proportional to an eigenfunction of the Laplace operator in this confining geometry considerably simplifies the underlying mathematics. In particular, exact and explicit relations for several moments of the total phase accumulated by diffusing spins are derived. These relations are shown to provide good approximations for the typical case of a linear magnetic field gradient, for which the theoretical analysis was in general limited to the second moment. We study the structure and the properties of the higher order moments which are responsible for the breakdown of the "Gaussian phase approximation" (GPA) at intense magnetic fields. The limits of applicability of the GPA for nonlinear magnetic fields and the transition to the localization regime are discussed. In particular, a diagram of different restricted diffusion regimes is presented.     

Vibrational relaxation in hydrogen bonded systems: I. Theory

A theory is proposed which expresses the influence of intermolecular interaction in H-bonded solutions on the solute IR band shape in terms of parameters characterising the modulation and strength of this interaction. It is found that appreciable vibrational relaxation can occur by interaction between oscillators with different frequencies, provided that the energy gap is closed by the hydrogen bond spectral density. Isotopic substitution of the solvent provides a means for experimental verification, which is the subject of a subsequent paper.     

The microscopic structure of an exactly solvable model binary solution that exhibits two closed loops in the phase diagram

An exactly solvable lattice model describing a binary solution is considered where rodlike molecules of types AA and BB cover the links of a honeycomb lattice, the neighboring molecular ends having three-body and orientation-dependent bonding interactions. At phase coexistence of AA-rich and BB-rich phases, the average fraction of each type of triangle of neighboring molecular ends is calculated exactly. The fractions of the different types of triangles are then used to deduce the local microscopic structure of the coexisting phases for a case of the model that contains two closed loops in the phase diagram.     

Ground-state energy of an s-state model of the inhomogeneous electron liquid in relation to an exactly solvable model of He with additional radial correlation

Following Temkin’s proposal of the so-called s-wave model, Amovilli, Howard and March (AHM) displayed an exactly solvable model Hamiltonian in which an additional radial correlation is added. The ground-state wave function for this Hamiltonian, for modelling He-like atomic ions with nuclear charge Ze, is used here to compare and contrast with the Temkin Model. The differences only appear near the critical charges Z_ce at which one electron ionises, Z_c being precisely unity for the AHM model and having the value 0.948768 in Serra’s variational study of the s-wave model. These two models are then compared with He-like ions having the correct e²/r₁₂ interaction.