Blueshift and intramolecular tunneling of NH3 umbrella mode in 4He n clusters

We present diffusion Monte Carlo calculations of the ground and first excited vibrational states of NH(3) (4)He(n) for n< or =40. We use the potential energy surface developed by one of us [M. P. Hodges and R. J. Wheatley, J. Chem. Phys. 114, 8836 (2001)], which includes the umbrella mode coordinate of NH(3). Using quantum Monte Carlo calculations of excited states, we show that this potential is able to reproduce qualitatively the experimentally observed effects of the helium environment, namely, a blueshift of the umbrella mode frequency and a reduction of the tunneling splittings in ground and first excited vibrational states of the molecule. These basic features are found to result regardless of whether dynamical approximations or exact calculations are employed.     

Electronically excited rubidium atom in helium clusters and films. II. Second excited state and absorption spectrum

Following our work on the study of helium droplets and film doped with one electronically excited rubidium atom Rb(?) ((2)P) [M. Leino, A. Viel, and R. E. Zillich, J. Chem. Phys. 129, 184308 (2008)], we focus in this paper on the second excited state. We present theoretical studies of such droplets and films using quantum Monte Carlo approaches. Diffusion and path integral Monte Carlo algorithms combined with a diatomics-in-molecule scheme to model the nonpair additive potential energy surface are used to investigate the energetics and the structure of Rb(?)He(n) clusters. Helium films as a model for the limit of large clusters are also considered. As in our work on the first electronic excited state, our present calculations find stable Rb(?)He(n) clusters. The structures obtained are however different with a He-Rb(?)-He exciplex core to which more helium atoms are weakly attached, preferentially on one end of the core exciplex. The electronic absorption spectrum is also presented for increasing cluster sizes as well as for the film.     

Quantum Monte Carlo calculations of the potential energy curve of the helium dimer

We report results of two quantum Monte Carlo methods -- variational Monte Carlo and diffusion Monte Carlo -- on the potential energy curve of the helium dimer. In contrast to previous quantum Monte Carlo calculations on this system, we have employed trial wave functions of the Slater-Jastrow form and used the fixed node approximation for the fermion nodal surface. We find both methods to be in excellent agreement with the best theoretical results at short range. In addition, the diffusion Monte Carlo results give very good agreement across the whole potential energy curve, while the Slater-Jastrow wave function fails to bind the dimer at all.     

剩余函数量子Monte Carlo方法

为量子Ｍｏｎｔｅ Ｃａｒｌｏ方法提出一条新途径－剩余函数法,引入了Ｓｃｈｒｏｅｄｉｎｇｅｒ方程剩余函数的概念,利用剩余函数将一种新的有明显物理意义的试探函数应用到量子Ｍｏｎｔｅ Ｃａｒｌｏ过程中,这种试探函数是通过一种迭进式的方式确定的,它不需要在Ｍｏｎｔｅ Ｃａｒｌｏ过程中优化参数。文中我们将给出这种试探函数的具体形式,证明由这种试探函数求出的能量期望值收敛于体系真实的能量值;文中还给出这种试探     

Complex excited dynamics around a plateau on a retinal-like potential surface: chaos, multi-exponential decays and quantum/classical differences

We investigate the classical and quantum dynamics on the plateau of an excited potential energy surface (PES) whose shape mimics the PES driving the photoisomerization of the protonated Schiff base of retinal (PSBR). We adopt a two-dimensional analytical model of the PES, and perform an extended study by varying the potential parameters, revealing a scenario whose interest goes beyond the relevance for the specific case of PSBR. In fact, we document cases with net differences among classical and quantum dynamical predictions, for barrierless PESs. Classical trajectories released on the PES display the signature of chaos and partial trapping on the plateau, whose origin is purely dynamical, since no barrier exists. At variance, on the same barrierless PESs, quantum dynamics does not predict any trapping, always showing a complete depletion of the excited population according to an approximate mono-exponential law. The plateau on the PES promotes complex and unusual dynamical features, and it is sufficient to introduce a very small barrier along the cis–trans torsional mode to give rise to a multi-exponential decay, also at quantum level. Our results are of general interest because plateaux are often found in the excited states of conjugated chromophores.     

Quantum Monte Carlo for electronic excitations of free-base porphyrin

Accurate calculations of allowed and nonallowed transitions in porphyrin are reported. Using the quantum Monte Carlo method in the diffusion Monte Carlo variant, the vertical transition between the ground state singlet and the second excited state singlet as well as the adiabatic transition between the ground state and the lowest triplet state have been computed for this 162-electron system. The present theoretical results are compared to experiment and to results of other theoretical methods. The diffusion Monte Carlo energy differences are found to be in excellent agreement with experiment.     

A rare event sampling method for diffusion Monte Carlo using smart darting

We identify a set of multidimensional potential energy surfaces sufficiently complex to cause both the classical parallel tempering and the guided or unguided diffusion Monte Carlo methods to converge too inefficiently for practical applications. The mathematical model is constructed as a linear combination of decoupled Double Wells [(DDW)(n)]. We show that the set (DDW)(n) provides a serious test for new methods aimed at addressing rare event sampling in stochastic simulations. Unlike the typical numerical tests used in these cases, the thermodynamics and the quantum dynamics for (DDW)(n) can be solved deterministically. We use the potential energy set (DDW)(n) to explore and identify methods that can enhance the diffusion Monte Carlo algorithm. We demonstrate that the smart darting method succeeds at reducing quasiergodicity for n ? 100 using just 1 × 10(6) moves in classical simulations (DDW)(n). Finally, we prove that smart darting, when incorporated into the regular or the guided diffusion Monte Carlo algorithm, drastically improves its convergence. The new method promises to significantly extend the range of systems computationally tractable by the diffusion Monte Carlo algorithm.     

Computational analysis of the proton translocation from Asp96 to schiff base in bacteriorhodopsin

The potential energy change during the M --> N process in bacteriorhodopsin has been evaluated by ab initio quantum chemical and advanced quantum chemical calculations following molecular dynamics (MD) simulations. Many previous experimental studies have suggested that the proton transfer from Asp96 to the Schiff base occurs under the following two conditions: (1) the hydrogen bond between Thr46 and Asp96 breaks and Thr46 is detached from Asp96 and (2) a stable chain of four water molecules spans an area from Asp96 --> Schiff base. In this work, we successfully reproduced the proton-transfer process occurring under these two conditions by molecular dynamics and quantum chemical calculations. The quantum chemical computation revealed that the proton transfer from Asp96 to Shiff base occurs in two-step reactions via an intermediate in which an H(3)O(+) appears around Ala215. The activation energy for the proton transfer in the first reaction was calculated to be 9.7 kcal/mol, which enables fast and efficient proton pump action. Further QM/MM (quantum mechanical/molecular mechanical) and FMO (fragment molecular orbital) calculations revealed that the potential energy change during the proton transfer is tightly regulated by the composition and the geometry of the surrounding amino acid residues of bacteriorhodopsin. Here, we report in detail the Asp96 --> Schiff base proton translocation mechanism of bacteriorhodopsin. Additionally, we discuss the effectiveness of combining quantum chemical calculations with truncated cluster models followed by advanced quantum chemical calculations applied to a whole protein to elucidate its reaction mechanism.     

Self-optimizing Diffusion Quantum Monte Carlo Calculation: the Potential Energy Curve of C_2

11NTRODUCTIONDiffusionquantumMonteCarlo(DMC)isoneofthesimplestofthevariousMonteCarlotechniquesavailabletosolvetheSchrodingerequation,forarecentre-viewofDMC(seeRef[1i).Foravarietyofsmallatomsandmolecules,DMChasbeenshowntobecapableofprovidinganestimateoftheground-state(nonrelativistic)energywithanaccuracycomparabletogoodqualityClcalculations,evenwhenarelativelysimpletrialwavefunctionisemployed.However,todate,thereareseveralobstaclesinaDMCcalculation:(1)BeforetheDMCcomputation,thepar…     

Nuclear quantum effect on the hydrogen-bonded structure of guanine�Ccytosine pair

The structure of Watson?CCrick type guanine?Ccytosine (G?CC) base pair has been studied by classical hybrid Monte Carlo (HMC) and quantum path integral hybrid Monte Carlo (PIHMC) simulations on the semiempirical PM6 potential energy surface. For the three NH?X hydrogen-bonded moieties, the intramolecular NH bonds are found systematically longer while the H?X distance shorter in the PIHMC simulation than in the HMC simulation. We found that the hydrogen bonded length N?X correlates with the H?X distance, but not with the NH distance. A correlation is also between the neighboring hydrogen bonds in the G?CC base pair.     

精确固定节面量子Monte：Carlo差值法

本文提出了精确固定节面量子Monte Carlo差值法，这个新算法能够在精确固定节面量子Monte Carlo方法的基础上直接计算两个体系之间的能量差，且使计算结果的统计误差达到10-2 kJ/mol 数量级，获得电子相关能90%以上。我们把这个新算法应用于分子势能面的研究中，使用一个“刚性移动”模型，利用Jacobi变换使分子两个几何构型的能量计算具有很好的正相关性，因而能得到准确的能量差值，于是精确的分子势能面就可以得到。这个新算法已经被使用到BH分子基态势能曲线和H₃分子势能面的研究。这个算法还可应用于分子光谱、化学反应能量变化值等领域的研究。     

Optimized Jastrow-Slater wave functions for ground and excited states: application to the lowest states of ethene

A quantum Monte Carlo method is presented for determining multideterminantal Jastrow-Slater wave functions for which the energy is stationary with respect to the simultaneous optimization of orbitals and configuration interaction coefficients. The approach is within the framework of the so-called energy fluctuation potential method which minimizes the energy in an iterative fashion based on Monte Carlo sampling and a fitting of the local energy fluctuations. The optimization of the orbitals is combined with the optimization of the configuration interaction coefficients through the use of additional single excitations to a set of external orbitals. A new set of orbitals is then obtained from the natural orbitals of this enlarged configuration interaction expansion. For excited states, the approach is extended to treat the average of several states within the same irreducible representation of the pointgroup of the molecule. The relationship of our optimization method with the stochastic reconfiguration technique by Sorella et al. is examined. Finally, the performance of our approach is illustrated with the lowest states of ethene, in particular with the difficult case of the 1(1)B(1u) state.     

PROTON TRANSFER IN THE PRIMARY PROCESS OF VISION

Abstract— Proton transfer was theoretically examined as a possible primary process of vision. The motion of protons in the adiabatic potential of the Schiff base hydrogen bond was investigated in terms of quantum mechanics. The probability of proton transfer from the Schiff base nitrogen (i.e. the unprotonation of Schiff base) was found to increase as the retinal rotated around 11–12. double bond by 90°. The results also suggested that the proton transfer can take place before or during the transition from the excited to ground state (excited state proton transfer). We proposed that such excited state proton transfer is one of the elementary processes in primary visual photochemistry, and this process leads to the unprotonated visual pigment, hyposorhodopsin, which has been experimentally verified as one of the primary photoproducts of rhodopsin. The probability of this process could be comparable to the conventional process leading to the protonated intermediate, bathorhodopsin. The relation of these results with the recent experimental data is discussed.     

The electronic structure of low-lying excited states of two model systems of retinal

2,4-pentadienal and 2,4,6,8-nonatetraenal were studied as simple model systems of retinal. Four kinds of CI were performed on low-lying excited states of 2,4-pentadienal by using a split valence basis set. The results show that MR SD π CI is not adequate for the π–π* state and the single excitation σπ CI is a good compromise between cost and effectiveness as far as singly excited states are concerned. This CI was applied to the bigger model system. All-trans and 11-cis forms of aldehyde, Schiff base, and protonated Schiff base of the model system were studied. A fairly large energy lowering of about 1 eV of the first allowed excited state (π → π*) occurs when the Schiff base is protonated for both all-trans and 11-cis forms.     

多维微正则蒙特卡罗方法计算异氰酸S0→T1系间窜跃速率

推导了通用的坐标空间多维蒙特卡罗微正则非绝热速率计算公式,并应用于计算异氰酸基5（S0）到第一激发三态（T1）的系间窜跃（ISC）速率．在CCSD水平上构建了一,二,三维势能面用于蒙特卡罗抽样．计算所得S0→T1 ISC速率与实验符合较好,因此可以预期多维蒙特卡罗微正则非绝热速率理论将成为计算多原子分子非绝热反应速率的有效手段之一．     

Gamma distribution model to provide a direct assessment of the overall quality of quantum Monte Carlo-generated electron distributions

Our objective is to assess the accuracy of simulated quantum Monte Carlo electron distributions of atoms and molecules. Our approach is first to model the exact electron distribution by a linear combination of gamma distribution functions, with parameters chosen to exactly reproduce highly accurate literature values for a number of selected moments for the system of interest. In application to the ground-state electron distributions of helium and dihydrogen, a high level of accuracy of the model was confirmed upon comparing its predicted moments, not used in the model's parametrization, to those calculated from high-level theory. Next, we generated electron-electron and electron-nucleus distributions for dihydrogen from electron positions outputted from a variety of quantum Monte Carlo algorithms. Upon juxtaposition of the simulated distributions with the putatively exact one that we derived from the model, we quantified the error in simulated distributions. The most accurate distributions were obtained from no-compromise reptation quantum Monte Carlo, a recently developed algorithm designed to ameliorate the distributions' time-step bias. Marginally less accurate distributions were generated from fixed-node diffusion Monte Carlo with descendant counting and detailed balance.     

Proton transfer reaction of a new orthohydroxy Schiff base in protic solvents at room temperature

Ground and excited state inter- and intramolecular proton transfer reactions of a new o-hydroxy Schiff base, 7-ethylsalicylidenebenzylamine (ESBA) have been investigated by means of absorption, emission and nanosecond spectroscopy in different protic solvents at room temperature and 77 K. The excited state intramolecular proton transfer (ESIPT) is evidenced by a large Stokes shifted emission (approximately 11000 cm(-1)) at a selected excited energy in alcoholic solvents. Spectral characteristics obtained reveal that ESBA exists in more than one structural form in most of the protic solvents, both in the ground and excited states. From the nanosecond measurements and quantum yield of fluorescence we have estimated the decay rate constants, which are mainly represented by nonradiative decay rates. At 77 K the fluorescence spectra are found to be contaminated with phosphorescence spectra in glycerol and ethylene glycol. It is shown that the fluorescence intensity and nature of the species present are dependent upon the excitation energy.     

Coupling quantum Monte Carlo to a nonlinear polarizable continuum model for spherical solutes

Starting from the nonlinear dielectric response model of Sandberg and Edholm, we derive an analytical expression of the polarization contribution to the solvation free energy in terms of the electronic density of the solute and the dielectric properties of the solvent. The solvent inhomogeneity is taken into account with the use of a smooth switching function whose spacial variation is established on the basis of how the solvent is arranged around the solute. An explicit form of a local potential representing the solvent effect on the solute is thus obtained by functional analysis. This effective potential can be combined with density functional or quantum chemical methods for the quantum mechanical treatment of the solute. Here, we use quantum Monte Carlo techniques for the solute and apply the method to the hydration of atomic ions finding very good agreement with experimental data.