A three-state model for the photophysics of guanine

The nonadiabatic photochemistry of the guanine molecule (2-amino-6-oxopurine) and some of its tautomers has been studied by means of the high-level theoretical ab initio quantum chemistry methods CASSCF and CASPT2. Accurate computations, based by the first time on minimum energy reaction paths, states minima, transition states, reaction barriers, and conical intersections on the potential energy hypersurfaces of the molecules lead to interpret the photochemistry of guanine and derivatives within a three-state model. As in the other purine DNA nucleobase, adenine, the ultrafast subpicosecond fluorescence decay measured in guanine is attributed to the barrierless character of the path leading from the initially populated (1)(pi pi* L-a) spectroscopic state of the molecule toward the low-lying methanamine-like conical intersection (gs/pi pi* L-a)(CI). On the contrary, other tautomers are shown to have a reaction energy barrier along the main relaxation profile. A second, slower decay is attributed to a path involving switches toward two other states, (1)(pi pi* L-b) and, in particular, (1)(n(o)pi*), ultimately leading to conical intersections with the ground state. A common framework for the ultrafast relaxation of the natural nucleobases is obtained in which the predominant role of a pi pi*-type state is confirmed.

The photoisomerization of retinal in bacteriorhodopsin: Experimental evidence for a three-state model

The primary events in the all-trans to 13-cis photoisomerization of retinal in bacteriorhodopsin have been investigated with femtosecond time-resolved absorbance spectroscopy. Spectra measured over a broad range extending from 7000 to 22,400 cm−1 reveal features whose dynamics are inconsistent with a model proposed earlier to account for the highly efficient photoisomerization process. Emerging from this work is a new three-state model. Photoexcitation of retinal with visible light accesses a shallow well on the excited state potential energy surface. This well is bounded by a small barrier, arising from an avoided crossing that separates the Franck–Condon region from the nearby reactive region of the photoisomerization coordinate. At ambient temperatures, the reactive region is accessed with a time constant of ≈500 fs, whereupon the retinal rapidly twists and encounters a second avoided crossing region. The protein mediates the passage into the second avoided crossing region and thereby exerts control over the quantum yield for forming 13-cis retinal. The driving force for photoisomerization resides in the retinal, not in the surrounding protein. This view contrasts with an earlier model where photoexcitation was thought to access directly a reactive region of the excited-state potential and thereby drive the retinal to a twisted conformation within 100–200 fs.

Thermodynamic analysis of three state denaturation of Peanut Agglutinin

Peanut Agglutinin (PNA) is a homotetrameric protein with a very unusual open quaternary structure. During denaturation, it first dissociates into a molten globule like state, which subsequently undergoes complete denaturation. Urea denaturation of PNA at neutral pH has been studied by intrinsic fluorescence spectroscopy and has been fitted to a three state model, A(4) double left right arrow 4I double left right arrow 4U, to get all the relevant thermodynamic parameters. Urea denaturation leads to continuous red shift of wavelength maxima. The molten globule like state is formed in a short range of urea concentration. Refolding of the denatured PNA has been attempted by intrinsic fluorescence study. Refolding by instantaneous dilution shows the occurrence of the formation of an intermediate at a relatively rapid rate, within few seconds. The transition from PNA tetramer to molten globule like state is found to have a Delta G value of similar to 33 kcal/mole while it is similar to 8 kcal/mole for the transition from molten globule like state to a completely denatured state. This in turn indicates that the tetramerization in PNA contributes significantly to the stability of the oligomer.

Aging and stationary properties of non-equilibrium symmetrical three-state models

We consider a non-equilibrium three-state model whose dynamics is Markovian and displays the same symmetry as the three-state Potts model, i.e. the transition rates are invariant under the cyclic permutation of the states. Unlike the Potts model, detailed balance is, in general, not satisfied. The aging and the stationary properties of the model defined on a square lattice are obtained by means of large-scale Monte Carlo simulations. We show that the phase diagram presents a critical line, belonging to the three-state Potts universality class, that ends at a point whose universality class is that of the Voter model. Aging is considered on the critical line, at the Voter point and in the ferromagnetic phase.

Three-body model of the final state interaction in heavy meson decay

We develop a relativistic three-body model for the final state interaction in D(+) -> K(-) pi(+) pi(+) decay based on the ladder approximation of the Bethe-Salpeter equation. The decay amplitude has the standard two-meson resonant transition matrix modulated by a spectator amplitude that carries a fully interacting three-body contribution and a smooth background term from the partonic amplitude. Our calculations are compared with D decay data.

Relativistic three-body model for final state interaction in D+ → K-π+π+ decay

A challenge in mesonic three-body decays of heavy mesons is to quantify the contribution of re-scattering between the final mesons. D decays have the unique feature that make them a key to light meson spectroscopy, in particular to access the Kn S-wave phase-shifts. We built a relativis-tic three-body model for the final state interaction in D+ → K -π+π+ decay based on the ladder approximation of the Bethe-Salpeter equation projected on the light-front. The decay amplitude is separated in a smooth term, given by the direct partonic decay amplitude, and a three-body fully interacting contribution, that is factorized in the standard two-meson resonant amplitude times a reduced complex amplitude that carries the effect of the three-body rescattering mechanism. The off-shell reduced amplitude is a solution of an inhomogeneous Faddeev type three-dimensional integral equation, that includes only isospin 1/2 K -π+ interaction in the S-wave channel. The elastic K-π+ scattering amplitude is parameterized according to the LASS data[1]. The integral equation is solved numerically and preliminary results are presented and compared to the experimental data from the E791 Collaboration[2, 3] and FOCUS Collaboration[4, 5].

Yang-Lee zeros of the two- and three-state Potts model defined on π3 Feynman diagrams

We present both analytical and numerical results on the position of partition function zeros on the complex magnetic field plane of the q=2 state (Ising) and the q=3 state Potts model defined on phi(3) Feynman diagrams (thin random graphs). Our analytic results are based on the ideas of destructive interference of coexisting phases and low temperature expansions. For the case of the Ising model, an argument based on a symmetry of the saddle point equations leads us to a nonperturbative proof that the Yang-Lee zeros are located on the unit circle, although no circle theorem is known in this case of random graphs. For the q=3 state Potts model, our perturbative results indicate that the Yang-Lee zeros lie outside the unit circle. Both analytic results are confirmed by finite lattice numerical calculations.

Asymptotic scattering and duality in the one-dimensional three-state quantum Potts model on a lattice

We determine numerically the single-particle and the two-particle spectrum of the three-state quantum Potts model on a lattice by using the density matrix renormalization group method, and extract information on the asymptotic (small momentum) S-matrix of the quasiparticles. The low energy part of the finite size spectrum can be understood in terms of a simple effective model introduced in a previous work, and is consistent with an asymptotic S-matrix of an exchange form below a momentum scale p*. This scale appears to vanish faster than the Compton scale, mc, as one approaches the critical point, suggesting that a dangerously irrelevant operator may be responsible for the behaviour observed on the lattice.

Critical behavior in the site-diluted three-dimensional three-state Potts model

We have studied numerically the effect of quenched site dilution on a weak first-order phase transition in three dimensions. We have simulated the site diluted three-states Potts model studying in detail the secondorder region of its phase diagram. We have found that the n exponent is compatible with the one of the three-dimensional diluted Ising model, whereas the h exponent is definitely different.

Mathematical model of COREX melter gasifier: Part I. Steady-state model

The COREX melter gasifier is a countercurrent reactor to produce liquid iron. Directly reduced iron (DRI), noncoking coal, and other additives are charged to the melter gasifier at their respective temperatures, and O-2 is blown through the tuyeres. Functionally, a melter gasifier is divided into three zones: a moving bed, fluidized bed, and free board. A model has been developed for the moving bed, where the tuyere region is two-dimensional (2-D) and the rest is one-dimensional (1-D). It is based on multiphase conservation of mass, momentum, and heat. The fluidized bed has been treated as 1-D. Partial equilibrium is calculated for the free board. The calculated temperature of the hot metal, the top gas, and the chemistry of the top gas agree with the reported plant data. The model has been used to study the effects of bed height, injection of impure O-2, coal chemistry, and reactivity on the process performance.

Influence of light in the mixed-layer on the parameters of a three-component model of phytoplankton size class

Phytoplankton size structure is an important indicator of the state of the pelagic ecosystem. Stimulated by the paucity of in situ observations on size structure, and by the sampling advantages of autonomous remote platforms, new efforts are being made to infer the size-structure of the phytoplankton from oceanographic variables that may be measured at high temporal and spatial resolution, such as total chlorophyll concentration. Large-scale analysis of in situ data has revealed coherent relationships between size-fractionated chlorophyll and total chlorophyll that can be quantified using the three-component model of Brewin et al. (2010). However, there are variations surrounding these general relationships. In this paper, we first revise the three-component model using a global dataset of surface phytoplankton pigment measurements. Then, using estimates of the average irradiance in the mixed-layer, we investigate the influence of ambient light on the parameters of the three-component model. We observe significant relationships between model parameters and the average irradiance in the mixed-layer, consistent with ecological knowledge. These relationships are incorporated explicitly into the three-component model to illustrate variations in the relationship between size-structure and total chlorophyll, ensuing from variations in light availability. The new model may be used as a tool to investigate modifications in size-structure in the context of a changing climate.

Estimation of electronic coupling in π -stacked donor-bridge-acceptor systems: correction of the two-state model

Comparison of donor-acceptor electronic couplings calculated within two-state and three-state models suggests that the two-state treatment can provide unreliable estimates of Vda because of neglecting the multistate effects. We show that in most cases accurate values of the electronic coupling in a π stack, where donor and acceptor are separated by a bridging unit, can be obtained as Ṽ da = (E2 - E1) μ12 Rda + (2 E3 - E1 - E2) 2 μ13 μ23 Rda2, where E1, E2, and E3 are adiabatic energies of the ground, charge-transfer, and bridge states, respectively, μij is the transition dipole moments between the states i and j, and Rda is the distance between the planes of donor and acceptor. In this expression based on the generalized Mulliken-Hush approach, the first term corresponds to the coupling derived within a two-state model, whereas the second term is the superexchange correction accounting for the bridge effect. The formula is extended to bridges consisting of several subunits. The influence of the donor-acceptor energy mismatch on the excess charge distribution, adiabatic dipole and transition moments, and electronic couplings is examined. A diagnostic is developed to determine whether the two-state approach can be applied. Based on numerical results, we showed that the superexchange correction considerably improves estimates of the donor-acceptor coupling derived within a two-state approach. In most cases when the two-state scheme fails, the formula gives reliable results which are in good agreement (within 5%) with the data of the three-state generalized Mulliken-Hush model

A three-compartment model of the hemodynamic response and oxygen delivery to brain

We describe a mathematical model linking changes in cerebral blood flow, blood volume and the blood oxygenation state in response to stimulation. The model has three compartments to take into account the fact that the cerebral blood flow and volume as measured concurrently using laser Doppler flowmetry and optical imaging spectroscopy have contributions from the arterial, capillary as well as the venous compartments of the vasculature. It is an extension to previous one-compartment hemodynamic models which assume that the measured blood volume changes are from the venous compartment only. An important assumption of the model is that the tissue oxygen concentration is a time varying state variable of the system and is driven by the changes in metabolic demand resulting from changes in neural activity. The model takes into account the pre-capillary oxygen diffusion by flexibly allowing the saturation of the arterial compartment to be less than unity. Simulations are used to explore the sensitivity of the model and to optimise the parameters for experimental data. We conclude that the three-compartment model was better than the one-compartment model at capturing the hemodynamics of the response to changes in neural activation following stimulation.

Positronium-hydrogen-atom scattering in a five-state model

The scattering of orthopositronium (Ps) by hydrogen atoms has been investigated in a five-state coupled-channel model allowing for Ps(1s)H(2s,2p) and Ps(2s,2p)H(1s) excitations using a recently proposed electron-exchange model potential. The higher (n greater than or equal to 3) excitations and ionization of the Ps atom are calculated using the first Born approximation. Calculations are reported of scattering lengths, phase shifts. elastic, Ps and H excitation, and total cross sections. Remarkable correlations are observed between the S-wave Ps-H binding energy and the singlet scattering length, effective range, and resonance energy obtained in various model calculations. These correlations suggest that if a Ps-H dynamical model yields the correct result for one of these four observables, it is expected to lead to the correct result for the other three. The present model, which is constructed so as to reproduce the Ps-H resonance at 4.01 eV, automatically yields a Ps-H bound state at - 1.05 eV that compares well with the accurate value of - 1.067 eV. The model leads to a singlet scattering length of 3.72a(0) and effective range of 1.67a(0), whereas the correlations suggest the precise values of 3.50a(0) and 1.65a(0) for these observables, respectively. [S1050-2947(99)07703-3].

Ps-H scattering using the three-state positronium close-coupling approximation

Positronium scattering off a hydrogen target has been studied employing a three-state positronium model close-coupling approximation (CCA) with and without electron exchange. Elastic, excitation and quenching cross sections are reported at low and medium energies. The effect of electron exchange is found to be significant at low energies. The ratio of quenching to the total cross section (the conversion ratio) approaches the value of 0.25 with increase of energy, as expected.