Some of the new unique features of the MOLCAS quantum chemistry package version 7 are presented in this report. In particular, the Cholesky decomposition method applied to some quantum chemical methods is described. This approach is used both in the context of a straight forward approximation of the two-electron integrals and in the generation of so-called auxiliary basis sets. The article describes how the method is implemented for most known wave functions models: self-consistent field, density functional theory, 2nd order perturbation theory, complete-active space self-consistent field multiconfigurational reference 2nd order perturbation theory, and coupled-cluster methods. The report further elaborates on the implementation of a restricted-active space self-consistent field reference function in conjunction with 2nd order perturbation theory. The average atomic natural orbital basis for relativistic calculations, covering the whole periodic table, are described and associated unique properties are demonstrated. Furthermore, the use of the arbitrary order Douglas-Kroll-Hess transformation for one-component relativistic calculations and its implementation are discussed. This section especially focuses on the implementation of the so-called picture-change-free atomic orbital property integrals. Moreover, the ElectroStatic Potential Fitted scheme, a version of a quantum mechanics/molecular mechanics hybrid method implemented in MOLCAS, is described and discussed. Finally, the report discusses the use of the MOLCAS package for advanced studies of photo chemical phenomena and the usefulness of the algorithms for constrained geometry optimization in MOLCAS in association with such studies.

to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added and two surface hopping algorithms are included to enable nonadiabatic calculations. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.

The primary event that initiates vision is the photoinduced isomerization of retinal in the visual pigment rhodopsin (Rh). Here, we use a scaled quantum mechanics/molecular mechanics potential that reproduces the isomerization path determined with multiconfigurational perturbation theory to follow the excited-state evolution of bovine Rh. The analysis of a 140-fs trajectory provides a description of the electronic and geometrical changes that prepare the system for decay to the ground state. The data uncover a complex change of the retinal backbone that, at Ϸ60-fs delay, initiates a space saving ''asynchronous bicycle-pedal or crankshaft'' motion, leading to a conical intersection on a 110-fs time scale. It is shown that the twisted structure achieved at decay features a momentum that provides a natural route toward the photoRh structure recently resolved by using femtosecond-stimulated Raman spectroscopy.photoisomerization ͉ rhodopsin ͉ vision T he visual pigment rhodopsin (Rh) (1, 2) is a G protein-coupled receptor containing a 11-cis retinal chromophore (PSB11) bounded to a lysine residue (Lys-296) via a protonated Schiff base linkage (see Fig. 1). While the biological activity of Rh is triggered by the light-induced 11-cis all-trans isomerization of PSB11, this reaction owes its efficiency (e.g., short time scale and high quantum yields) to the protein cavity (1). Recently, the mechanism of the isomerization of retinal in Rh has been investigated by using femtosecond-stimulated Raman spectroscopy (FSRS) (3). Kukura et al. (3) have reported on experimentally derived structures of photoRh and bathoRh, namely the first and second ground-state intermediates of the Rh photocycle.While such progress has provided information on the structural changes achieved 200 fs after light absorption, the faster structural changes prompting the excited-state decay of PSB11 (i.e., the central event of the isomerization mechanism) remain to be established. Indeed, it has been suggested that such decay may occur on a 60-fs time scale through fast hydrogen out-ofplane (HOOP) motion (3), whereas the traditional view points to a slower Ϸ150-fs decay driven by cis-trans isomerization motion (4). In principle, molecular dynamics simulations featuring a quantum chemical description of the chromophore can be used to address such issues. This fact was shown by Warshel (5) using semiempirical quantum chemistry to describe PSB11 and geometrical constraints to account for the protein environment. Later, Birge and Hubbard (6) reported a different semiempirical study of an explicit chromophore-counterion pair evolving along a single coordinate. While the first simulation of the retinal photoisomerization using a full atomic-level protein model (7) was reported for the related receptor bacterio-Rh (bR), attempts to simulate the PSB11 excited-state motion in a complete Rh model are more recent (8-10). On the other hand, a quantitative evaluation of the isomerization coordinate and time scale requires, as a prerequisite, an accurate excited-st...