• Keine Ergebnisse gefunden

(1)Notizen 505 Hydrogen Stretching Band of Hydrogen Bonded Complex: Molecular Dynamics Approach D

N/A
N/A
Protected

Academic year: 2022

Aktie "(1)Notizen 505 Hydrogen Stretching Band of Hydrogen Bonded Complex: Molecular Dynamics Approach D"

Copied!
2
0
0

Wird geladen.... (Jetzt Volltext ansehen)

Volltext

(1)

Notizen 505

Hydrogen Stretching Band of Hydrogen Bonded Complex: Molecular Dynamics Approach D. Pumpernik, B. Borštnik, and A. Ažman Chemical Institute Boris Kidrič, Ljubljana, Yugoslavia Z. Naturforsch. 3 8 a , 5 0 5 - 5 0 6 (1978);

received February 23, 1978

The molecular dynamics is used to evaluate the hydrogen stretching band in the hydrogen bonded system. The vibra- tional relaxation is the main mechanism influencing the band profile.

Rahman et al. [1] have recently reported a molecular dynamics study of water. This seems to be the first calculation avoiding the assumption of rigidness of the molecules by explicitly introducing the vibrational degree of freedom.

The purpose of this note is to show a new kind of information available through the molecular dy- namics simulation. As a model compound we studied the ( C H3)2S O - H F complex where we have been interesting in the band shape of the stretching vibration v(HF).

The potential of the system w as composed from pair-wise additive atom-atom potentials. Three kinds of potentials were used. Within (CH3)2SO we used harmonic potentials derived from infrared data [2]. The Lippincott-Schroeder potential [3]

was used for atoms which are directly involved in hydrogen bonds (0, H, F). The parameters of this potential were obtained by fitting Lippincott- Schroeder potential to ab-initio calculated quan- tities such as equilibrium distances, dissociation energies and vibrational frequencies of the hydrogen bonded dimers ( C H3)2S O - H F and H F - H F . The potential between (CH3)2SO molecules has been composed from Lennard-Jones potentials [4]. W e have studied 54 molecules of both kinds ((CH3)2SO, HF), at the concentration ratio 1:1, density

1 g/cm3 and room temperature. More details about the potentials and the strategy used in the molecular dynamics simulation may be found elsewhere [5].

The intensity I (w) of an infrared absorption is given by

oo

I(w) ~ $G(t)eiwtdt (1)

3600 3800 £000 £200 {[cm"]

Fig. 1. a) Band profile determined by Equation (1). b) Band profile calculated with the assumption of uncorrelated rotational and vibrational motions.

with G(t) = <(r(0) — ro) (r(t) — ro)) as an auto- correlation function of the vector distance r in the molecule HF, where ro is the equilibrium distance.

Figure l a gives the stretching bands v(HF). The band at higher frequency corresponds to unbonded HF molecules while the band at lower frequency belongs to HF molecules bonded into complexes.

If r is taken as a scalar quantity, the autocorrelation function (called vibrational correlation function Gx) does not include the rotational degree of freedom.

The stretching band v(HF) calculated with Gv does not differ significantly from the band in Figure l a . This seems to be a reasonable result since one expects that in hydrogen bonded systems the vibrational relaxation is the predominant mecha- nism in determining the band profile.

Very often it is claimed that rotational and vibrational motions are uncorrelated. With this assumption the intensity is

oo

I'(w) ~ jG^Gy^e^dt, (2)

Reprint requests to Dr. A. Azman, Chemical Institute Boris Kidric, Hajdrihova 19, 61001 Ljubljana, Jugoslawien.

where G\ (t) is the autocorrelation function

<u(0) « ( 0 )

(2)

506 Notizen

(u is the unit vector) and Gy is the vibrational correlation function. Figure 1 b gives the hydrogen stretching band deduced with the mentioned assumption. From the results it is obvious that the supposition of uncorrelated motions can not be valid. The motion of the hydrogen bonded complex is complicated to such an extent that any trivial

separation between the vibrational and the rota- tional degrees of freedom is impossible. Such a separation is also avoided in more advanced theories [6J.

A cknowledge merits

This work was suported by RC of Slovenia,

[1] A. Rahman, F. H. Stillinger, and H. L. Lemberg, J.

Chem. Phys. 63, 5223 (1975).

[2] M. F. Forel and M. Tranquille, Spectrochim. Acta 26 A, 1023 (1970).

[3] E. R, Lippincott and R. Schroeder, J. Chem. Phys. 23, 1099 (1955).

[4] J. O. Hirschfelder, Ch. F. Curtiss, and R. B. Bird, Molecular Theory of Gases and Liquids (John Wiley

1954). H. A. Scheraga, Advances Physical Organic Chemistry, Vol. 6., Ed. V. Gold (Academic Press, 1968), p . 103.

[5] D. Pumpernik, B. Borstnik, D. Janefcic, and A. Azman, to be published.

[6] S. Bratos, J. Chem. Phys. 63, 3499 (1975). S. Bratos and J. P. Chestier, Phys. Rev. A 9, 2136 (1974).

Nachdruck — auch auszugsweise — nur mit schriftlicher Genehmigung des Verlags Verantwortlich für den Inhalt: A. Klemm

Satz und Druck: Konrad Triltsch, Würzburg

Referenzen

ÄHNLICHE DOKUMENTE

The mechanism of the formation of [C,Hs]' ions by hydrogen rearrangement in the molecular ions of 1-phenylpropane and 1,3-diphenylpropane has been investigated by

Further, to examine the roles of K234, K235 and R237 in minimizing DNA leakage, several force-probe simulations were performed using a modified connector similar to the

We have presented how molecular dynamics (MD) simulations can benefit from machine learning (ML) potentials and provided some background for the imple-.. mentation of this

With respect to internal dynamics, useful information is obtained when splittings of the rotational transitions, due to tunneling motions connecting equivalent minima, are observed

The results of our calculation are in agreement with the above experimental data, indicating that the deuterium exchange found proceeds with about the same rate

The main objective of the present study is to use molecular dynamics computer simulations for quan- titative characterization of the local structural envi- ronments of ammonium ions

The smaller substituents in general lead to simpler packing patterns such as layer structures; more complex substituents can lead to three-dimensional patterns of great complexity.

As a result of the differ- ent hydrogen bonding ability of the counterions, in 1 and 2 ten (from a total of twelve) NH groups of each [Co(DMU) 6 ] 2+ building block participate