首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 0 毫秒
1.
    
Forward-backward semiclassical dynamics (FBSD) provides a rigorous and powerful methodology for calculating time correlation functions in condensed phase systems characterized by substantial quantum mechanical effects associated with zero-point motion, quantum dispersion, or identical particle exchange symmetries. The efficiency of these simulations arises from the use of classical trajectories to capture all dynamical information. However, full quantization of the density operator makes these calculations rather expensive compared to fully classical molecular dynamics simulations. This article discusses the convergence properties of various correlation functions and introduces an optimal Monte Carlo sampling scheme that leads to a significant reduction of statistical error. A simple and efficient procedure for normalizing the FBSD results is also discussed. Illustrative examples on model systems are presented.  相似文献   

2.
The nonradiative decay of a π-stacked pair of adenine molecules,one of which was excited by an ultrafast laser pulse,is studied by semiclassical dynamics simulations.This simulation investigation is focused on the effect of the formation of bonded excimer in stacked adenines on the mechanism of ultrafast decay.The simulation finds that the formation of the bonded excimer significantly lowers the energy gap between the LUMO and HOMO and consequently facilitates the deactivation of the electronically excited molecule.On the other hand,the formation of the chemical bond between two stacked adenines restricts the deformation vibration of the pyrimidine of the excited molecule due to the steric effect.This slows down the formation of the coupling between the HOMO and LUMO energy levels and therefore delays the deactivation process of the excited adenine molecule to the electronic ground state.  相似文献   

3.
We have simulated the photoisomerization dynamics of azobenzene, taking into account internal conversion and geometrical relaxation processes, by means of a semiclassical surface hopping approach. Both n-->pi* and pi-->pi* excitations and both cis-->trans and trans-->cis conversions have been considered. We show that in all cases the torsion around the N==N double bond is the preferred mechanism. The quantum yields measured are correctly reproduced and the observed differences are explained as a result of the competition between the inertia of the torsional motion and the premature deactivation of the excited state. Recent time-resolved spectroscopic experiments are interpreted in the light of the simulated dynamics.  相似文献   

4.
A computationally inexpensive energy correction is suggested for radicals described by the equation-of-motion coupled cluster method for ionized states in the singles and doubles approximation (EOMIP-CCSD). The approach is primarily intended for doublet states that are qualitatively described by Koopmans' approximation. Following a strategy similar to those used in multireference coupled cluster theory, the proposed correction accounts for all correlation effects through third order in perturbation theory and also includes selected contributions to higher-order energies. As an initial test of the numerical performance of the method, total energies and energy splittings are calculated for some small prototype radicals.  相似文献   

5.
A very efficient, strict nonstationary multi-term approach has been developed as a generalization of the conventional and the strict nonstationary two-term approximations used for solving the nonstationary , electron Boltzmann equation. As a first application the temporal relaxation of electrons in a model plasma acted upon by a de electric field has been investigated. The results are compared with corresponding ones obtained by the conventional and the strict nonstationary two-term approach as well as with very accurate Monte Carlo simulations. Perfect agreement between nonstationary eight-term Boltzmann and Monte Carlo calculations is found.  相似文献   

6.
Tilinin  I. S.  Werner  W. S. M. 《Mikrochimica acta》1994,114(1):485-503
The study of fast electron interaction with solids in the energy range from 100 eV to several tens of keV is prompted by quickly developing microbeam analysis techniques such as electron probe microanalysis, scanning electron microscopy, electron energy loss spectroscopy and so on. It turned out that for random solids the electron transport problem might be solved on the basis of the generalized radiative field similarity principle. The latter states that the exact differential elastic cross section in the kinetic equation may be replaced by an approximate one provided the conditions of radiative field similarity are fulfilled. Application of the generalized similarity principle to electron scattering in solids has revealed many interesting features of electron transport. Easy to use and effective formulae have been obtained for the angular and energy distribution of electrons leaving a target, total yields of characteristic photons and slow electrons escaping from a sample bombarded by fast primaries, escape probability of Auger electrons as a function of depth etc. The analytical results have been compared with Monte Carlo calculations and experiments in a broad range of electron energies and scattering properties of solids and good agreement has been observed.  相似文献   

7.
    
Ab initio electronic structure methods have reached a satisfactory accuracy for the calculation of static properties, but remain too expensive for quantum dynamical calculations. Recently, an efficient semiclassical method was proposed to evaluate the accuracy of quantum dynamics on an approximate potential without having to perform the expensive quantum dynamics on the accurate potential. Here, this method is applied for the first time to evaluate the accuracy of quantum dynamics on an approximate analytical or interpolated potential in comparison to the quantum dynamics on an accurate potential obtained by an ab initio electronic structure method. Specifically, the vibrational dynamics of H2 on a Morse potential is compared with that on the full CI potential, and the photodissociation dynamics of CO2 on a LEPS potential with that on the excited 1Π surface computed at the EOM‐CCSD/aug‐cc‐pVDZ level of theory. Finally, the effect of discretization of a potential energy surface on the quantum dynamics is evaluated. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110:2426–2435, 2010  相似文献   

8.
    
The UV–vis spectrum of Tyrosine and its response to different backbone protonation states have been studied by applying the Perturbed Matrix Method (PMM) in conjunction with molecular dynamics (MD) simulations. Herein, we theoretically reproduce the UV–vis absorption spectrum of aqueous solution of Tyrosine in its zwitterionic, anionic and cationic forms, as well as of aqua‐p‐Cresol (i.e., the moiety that constitutes the side chain portion of Tyrosine). To achieve a better accuracy in the MD sampling, the Tyrosine Force Field (FF) parameters were derived de novo via quantum mechanical calculations. The UV–vis absorption spectra are computed considering the occurring electronic transitions in the vertical approximation for each of the chromophore configurations sampled by the classical MD simulations, thus including the effects of the chromophore semiclassical structural fluctuations. Finally, the explicit treatment of the perturbing effect of the embedding environment permits to fully model the inhomogeneous bandwidth of the electronic spectra. Comparison between our theoretical–computational results and experimental data shows that the used model captures the essential features of the spectroscopic process, thus allowing to perform further analysis on the strict relationship between the quantum properties of the chromophore and the different embedding environments. © 2018 Wiley Periodicals, Inc.  相似文献   

9.
    
Binding free energies were calculated for the inhibitors lopinavir, ritonavir, saquinavir, indinavir, amprenavir, and nelfinavir bound to HIV-1 protease. An MMPB/SA-type analysis was applied to conformational samples from 3 ns explicit solvent molecular dynamics simulations of the enzyme-inhibitor complexes. Binding affinities and the sampled conformations of the inhibitor and enzyme were compared between different HIV-1 protease protonation states to find the most likely protonation state of the enzyme in the complex with each of the inhibitors. The resulting set of protonation states leads to good agreement between calculated and experimental binding affinities. Results from the MMPB/SA analysis are compared with an explicit/implicit hybrid scheme and with MMGB/SA methods. It is found that the inclusion of explicit water molecules may offer a slight advantage in reproducing absolute binding free energies while the use of the Generalized Born approximation significantly affects the accuracy of the calculated binding affinities.  相似文献   

10.
本文结合玻尔兹曼输运方程和电声散射速率计算研究锐钛矿和金红石二氧化钛中光生载流子的超快动力学过程. 其中,动力学模拟所需的结构参数均通过第一性原理计算获得. 结果表明,由于存在强Fröhlich型电声耦合,纵光学声子模对两个晶相的能量弛豫过程均有十分显著的影响,但是两个晶相的弛豫机理却表现出明显的差异. 对于单条导带内的弛豫过程,锐钛矿和金红石的能量弛豫时间分别为24.0 fs和11.8 fs,前者约为后者的二倍. 这一差异来源于两个晶相中不同的电子扩散分布以及不同的声学模散射贡献. 对于涉及多条导带的弛豫过程,预测的锐钛矿和金红石的总体弛豫时间分别为47 fs和57 fs,其相对大小与单条导带的情况相反. 分析表明金红石相弛豫较慢是因为存在多个速率控制步骤. 这些发现为调控电子动力学以及设计高效的二氧化钛器件提供了有价值的信息.  相似文献   

11.
The temporal evolution of the electron velocity-distribution function(EVDF), the concentration, mean energy, and the drift velocity of theelectrons is studied on a kinetic basis in a weakly ionized Ar/F2mixture plasma under conditions when the electron concentration temporallydecreases as a result of the electron attachment to fluoride molecules. Usingan appropriate relaxation model, the time-dependent electron Boltzmannequation was solved in multiterm and two-term approximations of the velocitydistribution function. The multiterm results confirmed predictions on theoccurrence of negative electron mobilities in such a decaying Ar/F2plasma, which were made in a former study using the conventional two-termapproximation. The investigations particularly showed that this approximationgives almost accurate results for the EVDF and related electron swarm parametersexcept for in the very beginning of the relaxation process. It has been furthershown that for a certain range of the reduced electric field strength, thedrift velocity becomes negative in the process of temporal evolution and remainsnegative even when approaching the hydrodynamic stage of the electronswarm. In addition, the role played by the back heating from the gas byelastic collisions on the EVDF formation is studied and various comparisonswith corresponding Monte Carlo results are performed.  相似文献   

12.
    
  相似文献   

13.
    
A fast stable finite difference Poisson-Boltzmann (FDPB) model for implicit solvation in molecular dynamics simulations was developed using the smooth permittivity FDPB method implemented in the OpenEye ZAP libraries. This was interfaced with two widely used molecular dynamics packages, AMBER and CHARMM. Using the CHARMM-ZAP software combination, the implicit solvent model was tested on eight proteins differing in size, structure, and cofactors: calmodulin, horseradish peroxidase (with and without substrate analogue bound), lipid carrier protein, flavodoxin, ubiquitin, cytochrome c, and a de novo designed 3-helix bundle. The stability and accuracy of the implicit solvent simulations was assessed by examining root-mean-squared deviations from crystal structure. This measure was compared with that of a standard explicit water solvent model. In addition we compared experimental and calculated NMR order parameters to obtain a residue level assessment of the accuracy of MD-ZAP for simulating dynamic quantities. Overall, the agreement of the implicit solvent model with experiment was as good as that of explicit water simulations. The implicit solvent method was up to eight times faster than the explicit water simulations, and approximately four times slower than a vacuum simulation (i.e., with no solvent treatment).  相似文献   

14.
15.
    
The Poisson-Boltzmann equation is widely used to describe the electrostatic potential of molecules in an ionic solution that is treated as a continuous dielectric medium. The linearized form of this equation, applicable to many biologic macromolecules, may be solved using the boundary element method. A single-layer formulation of the boundary element method, which yields simpler integral equations than the direct formulations previously discussed in the literature, is given. It is shown that the electrostatic force and torque on a molecule may be calculated using its boundary element representation and also the polarization charge for two rigid molecules may be rapidly calculated using a noniterative scheme. An algorithm based on a fast adaptive multipole method is introduced to further increase the speed of the calculation. This method is particularly suited for Brownian dynamics or molecular dynamics simulations of large molecules, in which the electrostatic forces must be calculated for many different relative positions and orientations of the molecules. It has been implemented as a set of programs in C++, which are used to study the accuracy and speed of this method for two actin monomers.  相似文献   

16.
17.
胰蛋白酶和苯酰氨类抑制剂结合自由能的预测   总被引:1,自引:0,他引:1  
用基于线性响应近似的自由能预测方法计算胰蛋白酶和苯酰氨类抑制剂的结合 自由能。计算结果表明,单参数,双参数和三参数模型具有相似的线性回归系数, 但三参数和双参数模型的交互验证回归系数要明显优于单参数模型。从预测能力来 看,双参数模型和三参数模型都能够很好地预测测试集中抑制剂的结合自由能,其 中双参数模型预测的结果要略优于三参数模型的预测结果。对测试集中的抑制剂, 双参数模型预测得到的预测自由能和实际自由能之间平均绝对误差仅为1.15 kJ/mol。自由能计算模型以及分子动力学轨迹能很好地解释抑制剂结构和活性的 关系,为药物设计提供了重要的结构信息。  相似文献   

18.
    
IBIsCO is a parallel molecular dynamics simulation package developed specially for coarse-grained simulations with numerical potentials derived by the iterative Boltzmann inversion (IBI) method (Reith et al., J Comput Chem 2003, 24, 1624). In addition to common features of molecular dynamics programs, the techniques of dissipative particle dynamics (Groot and Warren, J Chem Phys 1997, 107, 4423) and Lowe-Andersen dynamics (Lowe, Europhys Lett 1999, 47, 145) are implemented, which can be used both as thermostats and as sources of friction to compensate the loss of degrees of freedom by coarse-graining. The reverse nonequilibrium molecular dynamics simulation method (Müller-Plathe, Phys Rev E 1999, 59, 4894) for the calculation of viscosities is also implemented. Details of the algorithms, functionalities, implementation, user interfaces, and file formats are described. The code is parallelized using PE_MPI on PowerPC architecture. The execution time scales satisfactorily with the number of processors.  相似文献   

19.
Kunz AP  van Gunsteren WF 《Chemphyschem》2011,12(14):2609-2614
A method for conformational Boltzmann sampling of loops in proteins in aqueous solution is presented that is based on adiabatic decoupling molecular dynamics (MD) simulation with temperature or force scaling. To illustrate the enhanced sampling, the loop from residues 33 to 43 in the bovine protein ribonuclease A is adiabatically decoupled from the rest of the protein and the solvent with a mass scaling factor s(m) =1000 and the sampling is enhanced with a scaling of the temperature using s(T) =2 or of the force using s(V) =0.667. Over 5 ns of simulation the secondary structure of the protein remains unaltered while a combined dihedral-angle conformational cluster analysis shows an increase of conformations outside the first most populated cluster of loop conformations for adiabatic decoupling MD with temperature scaling using s(T) =2 or force scaling using s(V) =0.667 compared to the standard MD simulation. The atom-positional root-mean-square fluctuations of the C(α) atoms of the loop show an increase in the movement of the loop as well, indicating that adiabatic decoupling MD with upscaling of the temperature or downscaling of the force is a promising method for conformational Boltzmann sampling.  相似文献   

20.
    
Mesoscale molecular dynamics simulations are performed to analyze the curing process of an epoxy resin with polyfunctional amines on a generic surface. The coarse grained potentials were derived from all-atomistic molecular dynamics simulations using iterative Boltzmann inversion. The reactive scheme incorporates cross-linking between an epoxy resin and an amine, as well as amine adsorption on the surface. The structure of the cured network is examined and compared with equilibrium properties of the uncured system. Special attention has been paid on the implications of the surface that is believed to play a crucial role in the performance of epoxy systems.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号