Advanced
Vibrational Relaxation and Fragmentation in Icosahedral (Ar<sub>2</sub><sup>+</sup>)Ar<sub>12</sub> Clusters
Vibrational Relaxation and Fragmentation in Icosahedral (Ar2+)Ar12 Clusters
Bulletin of the Korean Chemical Society. 2014. Sep, 35(9): 2774-2780
Copyright © 2014, Korea Chemical Society
  • Received : March 06, 2014
  • Accepted : May 28, 2014
  • Published : September 20, 2014
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Jongbaik Ree
Yoo Hang Kim
Department of Chemistry, Inha University, Incheon 402-751, Korea
Hyung Kyu Shin
Department of Chemistry, University of Nevada, Reno, Nevada 89557, USA

Abstract
A dynamics study of relaxation and fragmentation of icosahedral argon cluster with a vibrationally excited ( v ) is presented. Local translation is shown to be responsible for inducing energy flow from the embedded ion to host atoms and fragmentation of the cluster consisting of various low frequency modes. The total potential energy of ( )Ar 12 is formulated using a building-up procedure of host-guest and host-host interactions. The time dependence of ion-to-host energy transfer is found to be tri-exponential, with the short-time process of ~100 ps contributing most to the overall relaxation process. Relaxation timescales are weakly dependent on both temperature (50-300 K) and initial vibrational excitation ( v = 1-4). Nearly 27% of host atoms in the cluster with ( v = 1) fragment immediately after energy flow, the extent increasing to ~43% for v = 4. The distribution of fragmentation products of ( )Ar 12 → ( )Ar n + (12- n )Ar are peaked around ( )Ar 8 . The distribution of dissociation times reveals fragmentation from one hemisphere dominates that from the other. This effect is attributed to the initial fragmentation causing a sequential perturbation of adjacent atoms on the same icosahedral five-atom layer.
Keywords
Introduction
The ions of inert gas clusters have been the subject of much research during the past decades both experimentally and theoretically. 1-19 These studies have focused on fragmentation, 1,4,14-16 thermochemistry, 5 ground-state geometries, 3,6,11,13,17 interaction potential energies 7-9 and energy transfer dynamics 10,19 among others. Small to medium-size clusters occupy an intermediate position between gas and condensed phases. Such clusters with embedded guest (or impurity) molecules provide an attractive model system for studying the physical and chemical processes, such as solvent effects, evaporation, melting or vibrational relaxation. 2,10,19 Many-body host-host and guest-host interactions induced by nonbonding interactions can lead to a potential energy surface with periodic local minima, which can influence the motions of guest molecules and the fragmentation of host atoms. Argon clusters provide a large enough space where a diatomic or small polyatomic molecule can be embedded to play its role on intra-cluster processes such as energy flow and subsequent fragmentation.
In this paper, we model a vibrationally excited
PPT Slide
Lager Image
ion embedded in an icosahedral argon cluster, (
PPT Slide
Lager Image
)Ar 12 , and study the energy flow from the ion to host atoms and the fragmentation process (
PPT Slide
Lager Image
)Ar 12 →(
PPT Slide
Lager Image
)Ar n + (12− n )Ar. The contents are organized as follows. The Model and Methods section describes the environment of the ion embedded in the first solvation shell, formulation of
PPT Slide
Lager Image
intramolecular and Ar-Ar intramolecular potential functions and method of solving the equations of motion. The extents of energy flow and fragmentation are presented in the Results and Discussion section followed by the Concluding Remarks section. Dynamics studies are carried out by integrating the equations of motion for 5000 clusters, which contain 60 000 host atoms and 10 000 atoms of the guest ion, over the temperature range of 50-300 K and initial vibrational excitation v = 1-4 of
PPT Slide
Lager Image
.
Model and Methods
Cluster Geometry. The host atoms in the icosahedral cluster (
PPT Slide
Lager Image
)Ar 12 are numbered as in Figure 1 , where the 6th and 12th atoms are at the axial positions on the Z axis. The diatomic ion
PPT Slide
Lager Image
is placed in parentheses. The n = 1-5 atoms are in the upper 5-atom layer at θ = 58° and ϕ = 0, 72, …, 288°, respectively, whereas n = 7-11 are in the lower 5- atom layer at θ = 122° and ϕ = 36, 108, …, 324°, respectively. Individual host atoms are identified as Ar(1), Ar(2), …, Ar(12), respectively. Subscript numbers are reserved for cluster sizes in discussing fragmentation products; e.g ., subscript 6 in (
PPT Slide
Lager Image
)Ar 6 of (
PPT Slide
Lager Image
)Ar n , where n = 1-12. The two atoms of
PPT Slide
Lager Image
are identified with a subscript as Ar ion(1) and Ar ion(2) to distinguish them from host atoms Ar( n ). A vibrationally excited
PPT Slide
Lager Image
embedded in the cluster interacts with 12 host atoms, while undergoing local translation ( δ ). In Figure 1 , we take Ar(5) as an example and denote its displacement from its equilibrium position by q 5 ; see the shaded Ar′(5). The instantaneous distance between Ar(5) and a nearest neighbor Ar( n ′) is ( d q 5 + q n′ ) or in general, ( d q n + q n′ ) between n th and n ′th atoms, where d is the equlibrium distance between the first-neighbors. The distance (− q n + q n′ ) represents the distortion of the (
PPT Slide
Lager Image
)Ar 12 con-figuration from the regular icosahedron geometry, the distortion which can become particularly severe in the presence of an embedded guest. In addition to 24 ion-host ( ion-h ) atom interactions, the cluster consists of 30 firstneighbor, 30 second-neighbor and six third-neighbor host-host ( h-h ) interactions.
PPT Slide
Lager Image
Numbering of host atoms and coordinate systems in icosahedral ()Ar12 around the cage center (c.c.) are shown. Each host atom is at one of 12 vertices. Interatomic distances and angles are defined. The displacement of the center-of-mass (c.m.) of is displaced from the c.c. by δ. A shaded circle is the displaced Ar(5) from its equilibrium position by q5.
Interaction Potential Energy. The distances from Ar′(5) to Ar ion(1) and Ar ion(2) are denoted by R Ar(5)-ion(1) and R Ar(5)-ion(2) , respectively. These ion-h distances are
PPT Slide
Lager Image
where r is the instantaneous bond distance of
PPT Slide
Lager Image
and c = −1. The same expression with c = 1 is for R Ar(n)-ion(2) . Here L n is the distance between the center of mass ( c.m .) of
PPT Slide
Lager Image
and n th host atom,
PPT Slide
Lager Image
and
PPT Slide
Lager Image
and 𝓵=3.892 Å represents the distance between the cluster center or cage center ( cc ) and the initial position of Ar(5) or the radius of the present cluster. 19 This distance forms a cage space large enough to embed
PPT Slide
Lager Image
, whose equilibrium bond distance is 2.4221 Å. 20 The explicit forms of the angle Θ n , defining the orientation of Ar n in the cage are
  • Θn=arccos{0.5257cosθ+ 0.8506sinθcos[2(n− 1)π/5 −ϕ]} forn= 1−5
  • Θn=arccos{0.5257cosθ− 0.8506sinθcos[λπ/5 −ϕ]} forn= 7−11 andλ= 1, 3, 5, 7, 9.
For n = 6 and 12, Θ n = θ . The angle factors cos ρ n and cosχ n in L n are replaced by the orientation averaged values 0 and ½, respectively. 21
The total interaction energy needed to describe the dynamic processes taking place in the cluster (
PPT Slide
Lager Image
)Ar 12 is a sum of the intramolecular potential of
PPT Slide
Lager Image
, host-host interactions, host oscillations and host-ion interactions including the polarization contribution,
PPT Slide
Lager Image
The first term is expressed in the Morse form
PPT Slide
Lager Image
where r e is the equilibrium bond distance 2.4221 Å 20 and the exponential range parameter calculated from
PPT Slide
Lager Image
, with D e = 1.3903 eV 7,20 and
PPT Slide
Lager Image
=307 cm −1 22 is 0.317 Å. When the latter potential parameters are used, Eq. (3) is in close agreement with the observed and computationally determined energies. 7,20,23
The host-host interaction potential is
PPT Slide
Lager Image
where the summation is over all pair-wise interactions and the Lennard-Jones ( LJ ) constants are ε = 143.224 K (or 0.01234 eV) and d = 3.757 Å. 24,25 For n = 1, the firstneighbors are Ar(2), Ar(5), Ar(6), Ar(7), Ar(11), and the second-neighbors are Ar(3), Ar(4), Ar(8), Ar(10), Ar(12), and the third-neighbor is Ar(9). It is important to establish the range of distances, which best describe the present cluster environment embedded with
PPT Slide
Lager Image
. For a regular icosahedron with 𝓵 = 3.892 Å, the nearest distance between the apex atoms is 4.092 Å. In this expanded cluster, host atoms oscillate around their nearest equilibrium positions separated from their first-neighbors by 3.757 Å. If we take a 4% displacement for each atom, the maximum and minimum distances between the atoms of the first-neighbor (or a dimer) are 4.057 and 3.456 Å, respectively. Thus the apex distance 4.092 Å in the expanded cluster creates a large enough shell area for the neighbors to move without suffering significant repulsive interactions. In the expanded cluster, the distance between the second-neighbors is 6.621 Å, whereas that between the third-neighbors is 2𝓵 = 7.784 Å.
We describe the oscillatory motion of each host atom around its equilibrium position in the form
PPT Slide
Lager Image
where D Ar measures the energy exerted by each neighbor, which we set equal to ε . The exponential range parameter is a Ar = ( D Ar /2 m ) 1/2 / ω Ar = 0.298 Å with the dimer frequency 30.68 cm −1 for ω Ar 26,27 and the reduced mass of m .
The dimer ion represented by Eq. (3) is quite stable ( D e = 1.39 eV), but its interaction with host atoms is comparatively much weaker. The main contribution to the ion-h interaction will be the polarization energy of Ar in the field of
PPT Slide
Lager Image
. Thus the interaction of the ion with 12 host atom can be expressed as a sum of the neutral Ar ion(1) -Ar and Ar ion(2) -Ar interactions at the ion-h distances defined in Eq. (1) and the polarization energy at the c.m .-to-host distance:
PPT Slide
Lager Image
where α Ar = 1.64 Å 3.28 As the distances in Eq. (6) are dependent on the coordinates set ( r , θ , ϕ , δ ,{ q }), the ion-h energy determines the coupling of all motions in the cluster, thus playing a central role in the dynamics processes.
Potential Energy Surfaces. Figure 2(a) shows the rotational potential energy surface (PES) for (
PPT Slide
Lager Image
)Ar 12 . The five peaks (0.0189 eV) appearing at θ = 58° are for the repulsive interactions of Ar (ion1) with the five upper-layer atoms located at the ϕ angle 72°-intervals. The maximum (0.0477 eV) at θ = 0° is for the interaction of Ar ion(1) with the axial atom Ar(6), whereas that at θ = 180° is for the interaction of Ar ion(2) with Ar(12). The five peaks fit nearly perfectly the 5-fold symmetry expression V h E h (1 + cos5 ϕ ) 1/2 , indicating that the rotational motion in the ϕ -direction is hindered strongly. Here E h = 0.363 eV is the peak height measured from the valley. The frequency of the hindered motion is ν h = ( n /2 π )( E h /2
PPT Slide
Lager Image
) 1/2 = 103 cm −1 , where n = 5 and
PPT Slide
Lager Image
is the moment of inertia of the ion. As θ approaches 90°, the interaction energy tends to the valley floor (−0.5478 eV). The perspective view between θ = 90° and 180° illustrates the interaction of Ar ion(2) with Ar(7), Ar(8), …, Ar(11) in the lower hemisphere. For comparison, we show the rotational PES for the one host atom case (
PPT Slide
Lager Image
)Ar, where only two repulsive peaks (0.646 eV at θ = 57°, ϕ = 72° and at θ = 122°, ϕ = 252°) are present; see Figure 2(b) . Here we consider the host atom Ar(2), but any other host atom results in the same topology except that the two peaks are oriented differently on the θ ϕ plane.
PPT Slide
Lager Image
(a) Rotational potential energy surface for ()Ar12. The peak height is 0.0189 eV and lowest valley depth is −0.5478 eV. (b) Rotational potential energy surface for ()Ar1. The peak height is 0.6458 eV and lowest valley depth is −0.0840 eV.
Equations of Motion. The equations of motion needed to solve for the time evolution of atom-atom distances and angles in the cluster (
PPT Slide
Lager Image
)Ar 12 are
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
where
PPT Slide
Lager Image
is the reduced mass of the ion and
PPT Slide
Lager Image
. We take a sample size of 5000 trajectories and integrate the equations of motion for 5 nanoseconds. Initial vibrational energies of
PPT Slide
Lager Image
correspond to the quantum values of 3/2 ħω , 5/2 ħω , 7/2 ħω and 9/2 ħω , which will be referred to v = 1, 2, 3 and 4, respectively. Numerical procedures using the standard routines based on the Adams-Moulton method 29 and randomsampling of initial conditions for r , θ , ϕ and { q } are given in ref. 21 .
Local translation mediates energy flow from vibration and rotation to host atoms, the process which can lead to fragmentation of (
PPT Slide
Lager Image
)Ar 12 . The rotation PES shown in Figure 2(a) indicates that the ion-h attraction dominates in the valleys between steeply rising peaks. The ion tends to travel along the attractive-interaction region undergoing an oscillatory motion. The motion describing the oscillatory behavior can be determined from the orientation-average of Eq. (1), which can be converted to a quadratic form to represent the oscillatory motion: < V ( δ )> = < V (0)> +½ |<∂ 2 V ( δ )/∂ δ 2 >| o δ 2 , where the second derivative is evaluated at r = r e , { q } = 0 an δ = 0. The solution of the equation of local translation described by the quadratic function gives 30
PPT Slide
Lager Image
PPT Slide
Lager Image
Where E is the energy of collision between the ion and host atoms.
Results and Discussion
Energy Relaxation. The solution of the equations of motion given above for the total interaction energy describes how the local translation, rotation and vibration of the ion, ion⋅⋅⋅host and host⋅⋅host distances evolve with time, and, in particular, their asymptotic behaviors at t →+∞. Before discussing the time evolution, we examine the dominant vibrations in the cluster; see Figure 3 , where the Fourier transforms of the temporal development of q n , δ and ion bond distance are shown. The host vibration is significantly structured between 15 and 60 cm −1 , the main peak appearing at 33 cm −1 . We note the presence of many other frequencies for the vibrations of host atoms; the peak shown here is a representative of them. The well-defined frequency of local translation obtained from the Fourier transform of the time dependence of δ is 121 cm −1 , which is close to
PPT Slide
Lager Image
obtained from Eq. (8-1). Also shown in Figure 3 is the 303 cm −1 peak for the embedded ion, which is red-shifted slightly from the fundamental 307 cm −1 of the free ion. 22 Although not shown in Figure 3 , we note the presence of yet another low frequency mode, the hindered rotational motion of a frequency 103 cm −1 . Thus the cluster is filled with low-frequency modes, the collection of which plays an essential role in promoting the vibrational relaxation of embedded ion and in turn, fragmentation of host atoms.
PPT Slide
Lager Image
Power spectra of ()Ar12 for Ar vibration (highest peak 33 cm−1), local translation (121 cm−1) and embedded (highest peak 303 cm−1).
We now show the time evolution of intra-cluster motions. Figure 4(a) displays several trajectories representing shorttime events occurring at 300 K for v = 1. Host atoms dissociating during the very early period are presented by trajectory A. Some host atoms dissociate after a short period of large amplitude motion (see C and D). Also shown is a representative of nondissociating events for comparison: see trajectory E. The time evolution of trajectory B is particularly interesting. For a clearer presentation, we reproduce it on an expanded scale; see Figure 4(b) . After some irregular motion during a very early period, the atom jumps to another site for a brief residency of ~70 ps, during which it performs a large amplitude oscillation before dissociating at t ≈ 100 ps. The
PPT Slide
Lager Image
-Ar( n ) distance of this transient state is close to twice the first-neighbor Ar-Ar distance, suggesting the formation of a second shell site for a brief period. We find that a significant number of trajectories undergo such site-hopping before fragmentation, while others lead to direct fragmentation as shown by trajectories A, C and D.
PPT Slide
Lager Image
(a) Time evolution of ion-h distances for trajectories representing dissociative (A-D) and nondissociative (E) events for (v = 1) at 300 K. (b) Plot of trajectory B on an expanded scale.
The ensemble-averaged energy transfer to host atoms closely fits the triple-exponential decay expression
PPT Slide
Lager Image
where for example, for the ion in v = 1 at 300 K shown in Figure 5 , c o = 0.0106 eV, c 1 = 0.0266 eV, c 2 = 0.0172 eV, c 3 = 0.0067 eV, k 1 = 0.101 ps −1 , k 2 = 0.014 ps −1 and k 3 = 1.05 ns −1 . A tri-exponential form is needed to describe a rapid decrease in the energy during the early period of interaction and a slow variation toward the limiting value at a long time (>1 ns). The relaxation rate constants k 1 and k 3 represent these short and long timescales of energy transfer process; i.e., k 1 −1 ≈ 10 ps, whereas k 3 −1 ≈ 1 ns, respectively. The intermediate timescale k 2 −1 ≈ 70 ps represents the period immediately following the early drop of vibrational energy, so it can be considered as part of the lower end of the short timescale. The energy c o = 0.0106 eV is the amount left in the ion in the limit t →∞ and is about half the zero-point energy 0.0190 eV of
PPT Slide
Lager Image
in the present quasi-classical calculation. For v = 4 also at 300 K, the vibrational energy tends to the limiting value c o = 0.0135 eV, more rapidly than the v = 1 case. Energy flow is shown to be still very active beyond 1 ns. In fact, local translation continues to pump vibrational energy even after time as long as 3 ns, when energy flow is no longer active in the low vibrational excitation case of v = 1.
PPT Slide
Lager Image
Time dependence of the vibrational energy of in ()Ar12 at 300 K for v = 1 and 4. The tri-exponential fitting curve for each case is shown by a dashed line, but the plot is not clearly seen as it overlaps nearly perfectly with the solid curve.
Figure 6(a) shows the temperature dependence of the relaxation rate constants. The value of k 1 increases only slightly from 7.83 × 10 10 at 50 K to 1.01 × 10 11 s −1 at 300 K. Similarly, the long-time constant k 3 increases from 4.64 × 10 8 s −1 at 50 K to 1.05 × 10 9 s −1 at 300 K. The temperature dependence is weak in each case, but the scales of short- and long-time energy relaxation processes differ by two orders of magnitudes. The effects of initial vibrational excitation on energy transfer at 300 K shown in Figure 6(b) are not significant. At 300 K, the rate constants vary weakly with the vibrational excitation of
PPT Slide
Lager Image
, also showing the longtime process is slower than the short-time event by two orders of magnitude. A negative dependence of k 3 on the initial vibrational energy results from the difficulty of transferring the last trace of vibrational energy from the highly excited ion to local translation.
PPT Slide
Lager Image
(a) Temperature dependence of short- and long-time relaxation rate constants (k1, k3) for the core ion (v = 1) embedded in ()Ar12. (b) Dependence of short- and long-time relaxation rate constants on initial vibrational excitation of the core ion (v) at 300 K.
Fragmentation. In the cage environment, local translation causes the embedded ion to collide continuously with host atoms which oscillate around their equilibrium positions. The energy released by the excited
PPT Slide
Lager Image
and thermal energy of local translation lead to the fragmentation of host atoms, which oscillate around their equilibrium positions. The fragmentation process can then be analyzed in terms of the time evolution of the distance between the host atom and the c.m. of the embedded ion; see Figure 4(a) . The extent of host atoms dissociated from 5000 (
PPT Slide
Lager Image
)Ar 12 clusters with
PPT Slide
Lager Image
( v = 1) increases from 25% at 50 K only slightly to 29% at 300 K. For v = 4, where a greater amount of vibrational energy is released to host atoms, the extents at 50 and 300 K are 39 and 43%, respectively. The axial atoms Ar(6) and Ar(12) positioned on the Z axis dissociate much more readily than others. At 300 K, 77 and 96% of Ar(6) and Ar(12) of clusters with
PPT Slide
Lager Image
( v = 1) dissociate, respectively, whereas fragmentation is significantly less for layer atoms. In the lower layer, the extent of fragmentation sequentially decreases from 85% for Ar(11), 48% for Ar(10), 18% for Ar(9), 5% for Ar(8) and 1% for Ar(7). On the other hand, the extent of fragmentation of host atoms in the upper layer is very low, the percentages being 13% for Ar(5) and only 2% of Ar(4). The host atoms Ar(3), Ar(2) and Ar(1) do not dissociate at all.
The above numbers indicate the extent of fragmentation from the lower hemisphere is nearly three times that from the upper hemisphere in the model shown in Figure 1 . The extent of fragmentation is nearly independent of temperature and vibrational excitation. The most frequently dissociating atoms are the two axial atoms and Ar(10) and Ar(11). Further, these atoms dissociate on a short-time scale characterized by the k 1 −1 timescale discussed above. For example, 65% of dissociative events of Ar(12) in 5000 (
PPT Slide
Lager Image
)Ar 12 with
PPT Slide
Lager Image
( v = 1) at 300 K occur within 100 ps from the start of ion-h interaction. Only 6% dissociate in 2 to 5 ns. As the initial burst of fragmentation consumes most of the available energy, the remaining bound atoms have a difficulty in gaining energy from various low frequency motions, thus taking a longer time to dissociate or do not dissociate at all. The most important mode of fragmentation is due to the energy flow from the embedded ion sequentially to host atoms Ar(11), Ar(10), Ar(9), Ar(8) and Ar(7) via the apex atom Ar(12), where collective motions of all lower-layer atoms facilitate energy flow and fragmentation induced by the apex atom. The apex atom Ar(6) also initiates a fragmentation string in the upper layer but its effect does not propagate long. This result leads to the formation of intermediate-sized clusters mainly composed of (undissociated) layer atoms in a short-time period. We note that when the icosahedron sketched in Figure 1 is rotated by 180° in the θ direction, Ar(1), Ar(2), …, Ar(6) now form the lower hemisphere, but the same fragmentation pattern is found as renaming the axes or renumbering host atoms does not lead to a different configuration. In fact, any host atom and its five nearest neighbors in an icosahedron form a hemisphere. The important point to be made is that fragmentation occurs efficiently from one hemisphere over the other via a chain of consecutive events from axial to layer atoms. For a fixed cluster structure, there are thus preferred sites and pathways of fragmentation facilitated by collective motions of the host atoms in one of the hemispheres at and below room temperature.
Figure 7 shows that the principal fragmentation products of (
PPT Slide
Lager Image
)Ar 12 →(
PPT Slide
Lager Image
)Ar n + (12- n )Ar are clusters with n between 7 and 10, the distribution peaking at n = 9 for v = 1. For
PPT Slide
Lager Image
( v = 1) at 300 K, 29% of host atoms dissociate as mentioned above, leaving 71% of them intact at various sites in the cluster. The majority of these bound atoms is in the n = 10, 9 and 8 clusters. The distribution of cluster sizes shown in Figure 7 is the number of (
PPT Slide
Lager Image
)Ar n clusters divided by the total number of (
PPT Slide
Lager Image
)Ar 12 sampled. The distribution curve for ( v = 1) at 300 K is peaked at the size 9 or (
PPT Slide
Lager Image
)Ar 9 . The distribution for the same vibrational excitation at 50 K is very similar to that at 300 K. Below n = 8, the probability decreases rapidly with decreasing cluster size. As noted above, after several host atoms dissociate, the available energy remaining in the cluster is low, so further dissociation forming smaller clusters becomes increasingly difficult. In fact, no clusters with n = 1 and 2 are found in all cases shown in Figure 7 . The figure also shows that the number of undissociated clusters (
PPT Slide
Lager Image
)Ar 12 is negligible, the probabilities being only 0.004 and 0.021 at 300 and 50 K, respectively. When the embedded ion is excited to v = 4, the amount of energy transfer to host atoms is significantly larger than the low excitation case (see Figure 5 ) and the probability peak shifts downward with the dominant products being (
PPT Slide
Lager Image
)Ar 7 and (
PPT Slide
Lager Image
)Ar 8 ; see Figure 7 for the v = 4 curves at 300 and 50 K. The size distribution shows that n = 12 and 11 clusters are now absent, while small but noticeable probabilities are found for n = 1 and 2, which reflects host atoms have difficulty of gaining the last trace of energy for dissociation, as noted in Figure 6(b) . For example, the probability of producing (
PPT Slide
Lager Image
)Ar 10 for v = 4 at 300 K is only 0.015, far less than the v = 1 probability 0.20, whereas the probabilities for small clusters such as (
PPT Slide
Lager Image
)Ar 2 are 0.013 and 0. Although their curves are not shown, we note that the cluster size distribution for v = 2 and 3 takes intermediate positions between the v = 1 and 4 curves shown in Figure 7 .
PPT Slide
Lager Image
Dependence of the fragmentation probability on cluster size for ()Arn. For initial excitation v = 1; filled circles with solid lines are at 300 K and filled squares with solid lines are at 50 K. For initial excitation v = 4: open circles with dashed lines are at 300 K and open squares with dashed lines are at 50 K.
Conclusions
We report the results of dynamics of intra-cluster energy flow and fragmentation of icosahedral (
PPT Slide
Lager Image
)Ar 12 with a vibrationally excited core ion over the temperature range of 50-300 K. The time evolution of vibrational energy transfer calculated in classical mechanics procedures using an analytic representation of the many-body interaction potential can be fit closely to a tri-exponential expression, which shows energy relaxation taking place over a wide range of time. The main portion of the vibrational energy of
PPT Slide
Lager Image
flows to host atoms via local translation on a picosecond scale and leads to fragmentation of the cluster. The dependence of relaxation rate constants on temperature and initial vibrational excitation of the diatomic ion is not strong.
About one-quarter of (
PPT Slide
Lager Image
)Ar 12 with the ion in v = 1 are found to fragment into various sizes in the temperature range 50-300 K. The extent of fragmentation nearly doubles for the ion in v = 4. The size distribution of fragmentation products is peaked for (
PPT Slide
Lager Image
)Ar 9 for v = 1 and is nearly independent of temperature. The peak shifts to (
PPT Slide
Lager Image
)Ar 8 for v = 4. The overall feature of size distribution varies significantly as vibrational excitation increases. Nearly all axial atoms Ar(6) and Ar(12) dissociate within a short-time period of ~100 ps. Further fragmentation is dominated sequentially by the layer atoms in one hemisphere of the icosahedron over another, which is a sign of regioselectivity showing the preference of energy flow in one region over all others. Thus efficient pathways of fragmentation in the vibrationally excited ionembedded cluster appear to be collective motions of nearest host atoms surrounding the apex site induced by the local translation of the core ion.
Acknowledgements
J. R. gratefully acknowledges the financial support from Chonnam National University.
References
Woodward C. A. , Upham J. E. , Stace A. J. , Murrell J. N. 1989 J. Chem. Phys. 91 7612 -    DOI : 10.1063/1.457283
Tanner G. M. , Bhattacharya A. , Nayak S. K. , Mahanti S. D. 1997 Phys. Rev. E 55 322 -    DOI : 10.1103/PhysRevE.55.322
Böhmer H.-U. , Peyerimhoff S. D. Z. 1989 Phys. D-Atoms, Molecules and Clusters 11 239 -    DOI : 10.1007/BF01438010
Snodgrass J. T. , Roehl C. M. , Bowers M. T. 1989 Chem. Phys. Lett. 159 10 -    DOI : 10.1016/S0009-2614(89)87444-5
Hiraoka K. , Mori T. 1989 J. Chem. Phys. 90 7143 -    DOI : 10.1063/1.456245
Gadea F. X. , Savrda J. , Paidarova I. 1994 Chem. Phys. Lett. 223 369 -    DOI : 10.1016/0009-2614(94)00453-6
Gadea F. X. , Paidarova I. 1996 Chem. Phys. 209 281 -    DOI : 10.1016/0301-0104(96)00107-3
Doltsinis N. L. , Knowles P. J. 1998 Mol. Phys. 94 981 -    DOI : 10.1080/002689798167557
Doltsinis N. L. , Knowles P. J. , Naumkin F. Y. 1999 Mol. Phys. 96 749 -
Ming L. , Markovic N. , Svanberg M. , Pettersson J. B. C. 1997 J. Phys. Chem. A 101 4011 -    DOI : 10.1021/jp964060w
Froudakis G. E. , Fanourgakis G. S. , Farantos S. C. , Xantheas S. S. 1998 Chem. Phys. Lett, 294 109 -    DOI : 10.1016/S0009-2614(98)00844-6
Greenblatt B. J. , Zanni M. T. , Neumark D. M. 1999 J. Chem. Phys. 111 10566 -    DOI : 10.1063/1.480433
Hrivnak D. , Kalus R. 2001 Chem. Phys. 264 319 -    DOI : 10.1016/S0301-0104(00)00363-3
Barat M. , Brenot J. C. , Fayeton J. A. , Picard Y. J. 2002 J. Chem. Phys. 117 1497 -    DOI : 10.1063/1.1485067
Ismail I. M. , Lepere V. , Barat M. , Fayeton J. A. , Picard Y. J. , Wohrer K. 2006 J. Chem. Phys. 124 164305 -    DOI : 10.1063/1.2191039
Lepere V. , Picard Y. J. , Barat M. , Fayeton J. A. , Lucas B. , Beroff K. 2009 J. Chem. Phys. 130 194301 -    DOI : 10.1063/1.3130974
Ritschel T. , Zuhrt C. , Zülicke L. , Kuntz P. J. 2007 Eur. Phys. J. D 41 127 -    DOI : 10.1140/epjd/e2006-00191-2
da Silva F. F. , Bartle P. , Denifl S. , Echt O. , Mark T. D. , Scheier P. 2009 Phys. Chem. Chem. Phys. 11 9791 -    DOI : 10.1039/b913175b
Shin H. K. 2011 J. Chem. Phys. 134 124301 -    DOI : 10.1063/1.3567958
Wüest A. , Merkt F. 2004 J. Chem. Phys. 120 638 -    DOI : 10.1063/1.1621618
Shin H. K. 2010 J. Chem. Phys. 132 104302 -    DOI : 10.1063/1.3339385
Signorell R. , Wüest A. , Merkt F. 1997 J. Chem. Phys. 107 10819 -    DOI : 10.1063/1.474199
Moseley J. T. , Saxon R. P. , Huber B. A. , Cosby P. C. , Abouaf R. , Tadjeddine M. 1977 J. Chem. Phys. 67 1659 -    DOI : 10.1063/1.434998
Aziz R. A. , Chen H. H. 1977 J. Chem. Phys. 67 5719 -    DOI : 10.1063/1.434827
Aziz R. A. 1993 J. Chem. Phys. 99 4518 -    DOI : 10.1063/1.466051
Dehmer P. M. , Dehmer J. L. 1978 J. Chem. Phys. 69 125 -    DOI : 10.1063/1.436396
Herman P. R. , LaRocque P. E. , Stoicheff B. P. 1988 J. Chem. Phys. 89 4535 -    DOI : 10.1063/1.454794
Bell R. J. , Kingston A. E. 1966 Proc. Phys. Soc. 88 901 -    DOI : 10.1088/0370-1328/88/4/311
1991 Fortran Subroutines for Mathematical Applications, Version 2.0 IMSL Houston 755 - 771
Shin H. K. 1974 Chem. Phys. Lett. 26 450 -    DOI : 10.1016/0009-2614(74)89070-6