Vibrational Relaxation and Fragmentation in Icosahedral (Ar_{2}^{+})Ar_{12} Clusters

Bulletin of the Korean Chemical Society.
2014.
Sep,
35(9):
2774-2780

- Received : March 06, 2014
- Accepted : May 28, 2014
- Published : September 20, 2014

Download

PDF

e-PUB

PubReader

PPT

Export by style

Article

Metrics

Cited by

TagCloud

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.
^{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
ion embedded in an icosahedral argon cluster, (
)Ar
_{12}
, and study the energy flow from the ion to host atoms and the fragmentation process (
)Ar
_{12}
→(
)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
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
.
Cluster Geometry.
The host atoms in the icosahedral cluster (
)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
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 (
)Ar
_{6}
of (
)Ar
_{n}
, where
n
= 1-12. The two atoms of
are identified with a subscript as Ar
_{ion(1)}
and Ar
_{ion(2)}
to distinguish them from host atoms Ar(
n
). A vibrationally excited
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 (
)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.
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
where
r
is the instantaneous bond distance of
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
and
n
th host atom,
and
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
, 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
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 (
)Ar
_{12}
is a sum of the intramolecular potential of
, host-host interactions, host oscillations and host-ion interactions including the polarization contribution,
The first term is expressed in the Morse form
where
r
_{e}
is the equilibrium bond distance 2.4221 Å
^{20}
and the exponential range parameter calculated from
, with
D
_{e}
= 1.3903 eV
^{7,20}
and
=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
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
. 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
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
. 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:
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 (
)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
)
^{1/2}
= 103 cm
^{−1}
, where
n
= 5 and
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 (
)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.
Equations of Motion.
The equations of motion needed to solve for the time evolution of atom-atom distances and angles in the cluster (
)Ar
_{12}
are
where
is the reduced mass of the ion and
. We take a sample size of 5000 trajectories and integrate the equations of motion for 5 nanoseconds. Initial vibrational energies of
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 (
)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}
Where
E
is the energy of collision between the ion and host atoms.
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
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.
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
-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.
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
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.
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
, 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.
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
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 (
)Ar
_{12}
clusters with
(
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
(
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 (
)Ar
_{12}
with
(
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 (
)Ar
_{12}
→(
)Ar
_{n}
+ (12-
n
)Ar are clusters with
n
between 7 and 10, the distribution peaking at
n
= 9 for
v
= 1. For
(
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 (
)Ar
_{n}
clusters divided by the total number of (
)Ar
_{12}
sampled. The distribution curve for (
v
= 1) at 300 K is peaked at the size 9 or (
)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 (
)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 (
)Ar
_{7}
and (
)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 (
)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 (
)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
.
)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
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 (
)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 (
)Ar
_{9}
for
v
= 1 and is nearly independent of temperature. The peak shifts to (
)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.

Introduction

The ions of inert gas clusters have been the subject of much research during the past decades both experimentally and theoretically.
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

Model and Methods

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

Numbering of host atoms and coordinate systems in icosahedral ()Ar_{12} 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 q _{5}.

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

- Θ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.

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

(a) Rotational potential energy surface for ()Ar_{12}. The peak height is 0.0189 eV and lowest valley depth is −0.5478 eV. (b) Rotational potential energy surface for ()Ar_{1}. The peak height is 0.6458 eV and lowest valley depth is −0.0840 eV.

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

Results and Discussion

PPT Slide

Lager Image

PPT Slide

Lager Image

Power spectra of ()Ar_{12} 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
PPT Slide

Lager Image

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

PPT Slide

Lager Image

PPT Slide

Lager Image

Time dependence of the vibrational energy of in ()Ar_{12} 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
PPT Slide

Lager Image

PPT Slide

Lager Image

(a) Temperature dependence of short- and long-time relaxation rate constants (k _{1}, k _{3}) for the core ion (v = 1) embedded in ()Ar_{12}. (b) Dependence of short- and long-time relaxation rate constants on initial vibrational excitation of the core ion (v ) at 300 K.

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

Dependence of the fragmentation probability on cluster size for ()Ar_{n}. 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

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

Acknowledgements

J. R. gratefully acknowledges the financial support from Chonnam National University.

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**

Citing 'Vibrational Relaxation and Fragmentation in Icosahedral (Ar_{2}^{+})Ar_{12} Clusters
'

@article{ JCGMCS_2014_v35n9_2774}
,title={Vibrational Relaxation and Fragmentation in Icosahedral (Ar_{2}^{+})Ar_{12} Clusters}
,volume={9}
, url={http://dx.doi.org/10.5012/bkcs.2014.35.9.2774}, DOI={10.5012/bkcs.2014.35.9.2774}
, number= {9}
, journal={Bulletin of the Korean Chemical Society}
, publisher={Korean Chemical Society}
, author={Ree, Jongbaik
and
Kim, Yoo Hang
and
Shin, Hyung Kyu}
, year={2014}
, month={Sep}