We performed first-principles calculations of the optical response of the green fluorescent protein (GFP) within a combined quantum-mechanical molecular-mechanics and time-dependent density-functional theory approach. The computed spectra are in excellent agreement with experiments assuming the presence of two, protonated and deprotonated, forms of the photoreceptor in a approximately 4:1 ratio, which supports the conformation model of photodynamics in GFP. Furthermore, we discuss charge transfer, isomerization, and environment effects. The present approach allows for systematic studies of excited-state electron-ion dynamics in biological systems.

Knowledge of the pK(a) of phosphoranes is important for the interpretation of phosphate ester hydrolysis. Calculated pK(a)'s of the model phosphorane, ethylene phosphorane, are reported. The method of calculation is based on the use of dimethyl phosphate as a reference state for evaluating relative pK(a) values, and on the optimization of the oxygen and acidic hydrogen van der Waals radii to give reasonable pK(1)(a), pK(2)(a), and pK(3)(a) for phosphoric acid in solution. Density functional theory is employed to calculate the gas-phase protonation energies, and continuum dielectric methods are used to determine the solvation corrections. The calculated pK(1)(a) and p(2)(a) for the model phosphorane are 7.9 and 14.3, respectively. These values are within the range of proposed experimental values, 6.5-11.0 for pK(1)(a), and 11.3-15.0 for pK(2)(a). The mechanistic implications of the calculated pK(a)'s are discussed.

An explicit formulation of the Piris cumulant λΔ,Π matrix is described herein, and used to reconstruct the two-particle reduced density matrix (2-RDM). Then, we have derived a natural orbital functional, the Piris Natural Orbital Functional 5, PNOF5, constrained to fulfill the D, Q, and G positivity necessary conditions of the N-representable 2-RDM. This functional yields a remarkable accurate description of systems bearing substantial (near)degeneracy of one-particle states. The theory is applied to the homolitic dissociation of selected diatomic molecules and to the rotation barrier of ethylene, both paradigmatic cases of near-degeneracy effects. It is found that the method describes correctly the dissociation limit yielding an integer number of electrons on the dissociated atoms. PNOF5 predicts a barrier of 65.6 kcal/mol for the ethylene torsion in an outstanding agreement with Complete Active Space Second-order Perturbation Theory (CASPT2). The obtained occupation numbers and pseudo one-particle energies at the ethylene transition state account for fully degenerate π orbitals. The calculated equilibrium distances, dipole moments, and binding energies of the considered molecules are presented. The values obtained are accurate comparing those obtained by the complete active space self-consistent field method and the experimental data.