Correlation molecular system simulations

The approach to the evaluation of vibrational spectra described above is based on classical simulations for which quantum corrections are possible. The incorporation of quantum effects directly in simulations of large molecular systems is one of the most challenging areas in theoretical chemistry today. The development of quantum simulation methods is particularly important in the area of molecular spectroscopy for which quantum effects can be important and where the goal is to use simulations to help understand the structural and dynamical origins of changes in spectral lineshapes with environmental variables such as the temperature. The direct evaluation of quantum time- correlation functions for anharmonic systems is extremely difficult. Our initial approach to the evaluation of finite temperature anharmonic effects on vibrational lineshapes is derived from the fact that the moments of the vibrational lineshape spectrum can be expressed as functions of expectation values of positional and momentum operators. These expectation values can be evaluated using extremely efficient quantum Monte-Carlo techniques. The main points are summarized below. 

The work described in this paper is an illustration of the potential to be derived from the availability of supercomputers for research in chemistry. The domain of application is the area of new materials which are expected to play a critical role in the future development of molecular electronic and optical devices for information storage and communication. Theoretical simulations of the type presented here lead to detailed understanding of the electronic structure and properties of these systems, information which at times is hard to extract from experimental data or from more approximate theoretical methods. It is clear that the methods of quantum chemistry have reached a point where they constitute tools of semi-quantitative accuracy and have predictive value. Further developments for quantitative accuracy are needed. They involve the application of methods describing electron correlation effects to large molecular systems. The need for supercomputer power to achieve this goal is even more acute. 

The strong point of molecular dynamic simulations is that, for the particular model, the results are (nearly) exact. In particular, the simulations take all necessary excluded-volume correlations into account. However, still it is not advisable to have blind confidence in the predictions of MD. The simulations typically treat the system classically, many parameters that together define the force field are subject to fine-tuning, and one always should be cautious about the statistical certainty. In passing, we will touch upon some more limitations when we discuss more details of MD simulation of lipid systems. We will not go into all the details here, because the use of MD simulation to study the lipid bilayer has recently been reviewed by other authors already [31,32]. Our idea is to present sufficient information to allow a critical evaluation of the method, and to set the stage for comparison with alternative approaches. 

Equilibrium properties are surprisingly accurately predicted by molecular-level SCF calculations. MC simulations help us to understand why the SCF theory works so well for these densely packed layers. In effect, the high density screens the correlations for chain packing and chain conformation effects to such a large extent that the properties of a single chain in an external field are rather accurate. Cooperative fluctuations, such as undulations, are not included in the SCF approach. Also, undulations cannot easily develop in an MD box. To see undulations, one needs to perform molecularly realistic simulations on very large membrane systems, which are extremely expensive in terms of computation time. 

Semiclassical techniques like the instanton approach [211] can be applied to tunneling splittings. Finally, one can exploit the close correspondence between the classical and the quantum treatment of a harmonic oscillator and treat the nuclear dynamics classically. From the classical trajectories, correlation functions can be extracted and transformed into spectra. The particular charm of this method rests in the option to carry out the dynamics on the fly, using Born Oppenheimer or fictitious Car Parrinello dynamics [212]. Furthermore, multiple minima on the hypersurface can be treated together as they are accessed by thermal excitation. This makes these methods particularly useful for liquid state or other thermally excited system simulations. Nevertheless, molecular dynamics and Monte Carlo simulations can also provide insights into cold gas-phase cluster formation [213], if a reliable force field is available [189]. 

Fig. 4.1 a Typical time evolution of a given correlation function in a glass-forming system for different temperatures (T >T2>...>T ), b Molecular dynamics simulation results [105] for the time decay of different correlation functions in polyisoprene at 363 K normalized dynamic structure factor at the first static structure factor maximum solid thick line)y intermediate incoherent scattering function of the hydrogens solid thin line), dipole-dipole correlation function dashed line) and second order orientational correlation function of three different C-H bonds measurable by NMR dashed-dotted lines)... 

Abstract. The physical nature of nonadditivity in many-particle systems and the methods of calculations of many-body forces are discussed. The special attention is devoted to the electron correlation contributions to many-body forces and their role in the Be r and Li r cluster formation. The procedure is described for founding a model potential for metal clusters with parameters fitted to ab initio energetic surfaces. The proposed potential comprises two-body, three-body, and four body interation energies each one consisting of exchange and dispersion terms. Such kind of ab initio model potentials can be used in the molecular dynamics simulation studies and in the cinalysis of binding in small metal clusters. 

The main advantage of the MFA is that it permits one to dramatically reduce the computational requisites associated with the study of solvent effects. This allows one to focus attention on the solute description, and it consequently becomes possible to use calculation levels similar to those usually employed in the study of systems and processes in the gas phase. Furthermore, in the case of ASEP/MD this high level description of the solute is combined with a detailed description of the solvent structure obtained from molecular dynamics simulations. Thanks to these features ASEP/MD [8] enables the study of systems and processes where it is necessary to have simultaneously a good description of the electron correlation of the solute and the explicit consideration of specific solute-solvent interactions, such as for VIS-UV spectra [9] or chemical reactivity [10]. 

One of the features that makes Equation (1) such a good starting point for our work is that it can be, in principal, exact. It is possible to show, without ever explicitly evaluating T and 37, that these crucial functions really do exist and are well defined (31). These formal definitions are rarely, if ever, useful in practical numerical calculations, but one can also work backwards from the exact dynamics x(t) (e.g., from a molecular dynamics simulation) to derive what the friction in particular must look like (32). The analysis tells us, moreover, that the exact T and are actually related to one another (28,31). The requirement that the relaxed system must be in equilibrium at some temperature T can be shown to set the magnitude and correlations of the fluctuating force ... 

A system of fundamental theoretical importance in many-body theory is the uniform-density electron gas. After decades of effort, exchange-correlation effects in this special though certainly not trivial system are by now well understood. In particular, sophisticated Monte Carlo simulations have provided very useful information (5) and have been conveniently parametrized by several authors (6). If the exchange-correlation hole function at a given reference point r in an atomic or molecular system is approximated by the hole function of a uniform electron gas with spin-densities given by the local values of p (r) and Pp(C obtain an... 

When partially hydrated samples are cooled down to 77 K, no crystallization peak is detected by differential thermal analysis. The x-ray and neutrons show that an amorphous form is obtained and its structure is different from those of low-and high-density amorphous ices already known [5]. Samples with lower levels of hydration corresponding to one monolayer coverage of water molecules are under investigation. This phenomenon looks similar in both hydrophilic and hydrophobic model systems. However, in order to characterize more precisely the nature of the amorphous phase, the site-site partial correlation functions need to be experimentally obtained and compared with those deduced from molecular dynamic simulations. 

---