Probing potential energy surfaces in confined systems: behavior of mean-square displacement in zeolites

Time evolution of mean-squared displacement based on molecular dynamics for a variety of adsorbate-zeolite systems is reported. Transition from ballistic to diffusive behavior is observed for all the systems. The transition times are found to be system dependent and show different types of dependence on temperature. Model calculations on a one-dimensional system are carried out which show that the characteristic length and transition times are dependent on the distance between the barriers, their heights, and temperature. In light of these findings, it is shown that it is possible to obtain valuable information about the average potential energy surface sampled under specific external conditions.

Tropical deep convection, convective available potential energy and sea surface temperature

The convective available potential energy (CAFE) based on monthly mean sounding has been shown to be relevant to deep convection in the tropics. The variation of CAFE with SST has been found to be similar to the variation of the frequency of deep convection at one station each in the tropical Atlantic and W. Pacific oceans. This suggests a strong link between the frequency of tropical convection and CAFE. It has been shown that CAFE so derived can be interpreted as the work potential of the atmosphere above the boundary layer with ascent in the convective region and subsidence in the surrounding cloud-free region.

Levitation effect and its relationship with the underlying potential energy landscape

Positions of potential energy minima for spherical monatomic sorbates in zeolite NaY have been identified for different sizes of the sorbate. It is found that for small sorbates (sigma less than or equal to 4.96 Angstrom) there are only six adsorption sites per alpha-cage, which are located close to the inner surface of the alpha-cage. For larger sorbates, additional sites of comparable energies appear close to the 12-ring window which forms the bottleneck for intercage diffusion. Minimum energy paths between these sites have been computed. These suggest that the barriers for both intracage and intercage site-to-site migrations are comparable and decrease with increase in sorbate size. Earlier simulation studies on the diffusion of monatomic sorbates in zeolites indicated that there is a dramatic change in the nature of dependence of D on sorbate size around 4.96 Angstrom, for zeolite NaY. Therefore, the present results suggest that the dependence of D on sorbate size and the changes in the potential energy landscape are correlated. The sorbate-zeolite system is characterized by a flatter potential energy landscape when the sorbate size is large. (C) 1999 American Institute of Physics. [S0021-9606(99)51110-0].

Characterization of the dynamics of glass-forming liquids from the properties of the potential energy landscape

We develop a framework for understanding the difference between strong and fragile behavior in the dynamics of glass-forming liquids from the properties of the potential energy landscape. Our approach is based on a master equation description of the activated jump dynamics among the local minima of the potential energy (the so-called inherent structures) that characterize the potential energy landscape of the system. We study the dynamics of a small atomic cluster using this description as well as molecular dynamics simulations and demonstrate the usefulness of our approach for this system. Many of the remarkable features of the complex dynamics of glassy systems emerge from the activated dynamics in the potential energy landscape of the atomic cluster. The dynamics of the system exhibits typical characteristics of a strong supercooled liquid when the system is allowed to explore the full configuration space. This behavior arises because the dynamics is dominated by a few lowest-lying minima of the potential energy and the potential energy barriers between these minima. When the system is constrained to explore only a limited region of the potential energy landscape that excludes the basins of attraction of a few lowest-lying minima, the dynamics is found to exhibit the characteristics of a fragile liquid.

A multi-sheeted three-dimensional potential-energy surface for the H-atom photodissociation of phenol

Extending the previous work of Lan et al. J. Chem. Phys., 122, 224315 (2005)], a multi-state potential model for the H atom photodissociation is presented. All three ``disappearing coordinates'' of the departing H atom have been considered. Ab initio CASSCF computations have been carried out for the linear COH geometry of C-2v symmetry, and for several COH angles with the OH group in the ring plane and also perpendicular to the ring plane. By keeping the C6H5O fragment frozen in a C-2v-constrained geometry throughout, we have been able to apply symmetry-based simplifications in the constructions of a diabatic model. This model is able to capture the overall trends of twelve adiabats at both torsional limits for a wide range of COH bend angles.

A comparative study of the energetics, structures, and mechanisms of the HCN H HNC and LiCN <-> LiNC isomerizations

The potential energy surfaces of the HCN<->HNC and LiCN<->LiNC isomerization processes were determined by ab initio theory using fully optimized triple-zeta double polarization types of basis sets. Both the MP2 corrections and the QCISD level of calculations were performed to correct for the electron correlation. Results show that electron correlation has a considerable influence on the energetics and structures. Analysis of the intramolecular bond rearrangement processes reveals that, in both cases, H (or Li+) migrates in an almost elliptic path in the plane of the molecule. In HCN<->HNC, the migrating hydrogen interacts with the in-plane pi,pi* orbitals of CN, leading to a decrease in the C-N bond order. In LiCN<->LiNC, Li+ does not interact with the corresponding pi,pi* orbitals of CN.

The Mechanism Of The Hydrogen Abstraction Process By Photo-Excited Ketones

Qualitative potential energy surfaces for hydrogen abstraction from alkanes containing primary, secondary and tertiary C-H bonds by a photo-excited ketone have been reported, The results suggest that the activation barriers for these processes decrease in the order primary > secondary > tertiary in agreement with the observed trend in the rate constants. The analysis of the electronic structure of the transition-state reveal that electron-transfer from hydrocarbon to ketone and formation of a new bond are almost synchronous in the hydrogen transfer process. The tunneling of hydrogen is not important in the normal temperature region even though the barriers are small.

Structure and Bonding in Stannadiphospholes and their Dianions SnC(2)P(2)R(2)(m) (R=H, tBu m=0,-2): A Comparative Study with C(5)H(5)(+) and C(5)H(5)(-) Analogues

The potential energy surfaces of both neutral and dianionic SnC(2)P(2)R(2) (R=H, tBu) ring systems have been explored at the B3PW91/LANL2DZ (Sn) and 6-311 + G* (other atoms) level. In the neutral isomers the global minimum is a nido structure in which a 1,2-diphosphocyclobutadiene ring (1,2-DPCB) is capped by the Sn. Interestingly, the structure established by Xray diffraction analysis, for R=tBu, is a 1,3-DPCB ring capped by Sn and it is 2.4 kcal mol(-1) higher in energy than the 1,2-DPCB ring isomer. This is possibly related to the kinetic stability of the 1,3-DPCB ring, which might originate from the synthetic precursor ZrCp(2)tBu(2)C(2)P(2). In the case of the dianionic isomers we observe only a 6 pi-electron aromatic structure as the global minimum, similarly to the cases of our previously reported results with other types of heterodiphospholes.([1,4,19]) The existence of large numbers of cluster-type isomers in neutral and 6 pi-planar structures in the dianions SnC(2)P(2)R(2)(2-) (R=H, tBu) is due to 3D aromaticity in neutral clusters and to 2D pi aromaticity of the dianionic rings. Relative energies of positional isomers mainly depend on: 1) the valency and coordination number of the Sn centre, 2) individual bond strengths, and 3) the steric effect of tBu groups. A comparison of neutral stannadiphospholes with other structurally related C(5)H(5)(+) analogues indicates that Sn might be a better isolobal analogue to P(+) than to BH or CH(+). The variation in global minima in these C(5)H(5)(+) analogues is due to characteristic features such as 1) the different valencies of C, B, P and Sn, 2) the electron deficiency of B, 3) weaker p pi-p pi bonding by P and Sn atoms, and 4) the tendency of electropositive elements to donate electrons to nido clusters. Unlike the C5H5+ systems, all C(5)H(5)(-) analogues have 6 pi-planar aromatic structures as global minima. The differences in the relative ordering of the positional isomers and ligating properties are significant and depend on 1) the nature of the pi orbitals involved, and 2) effective overlap of orbitals.

Theoretical analysis of photoinduced H-atom elimination in thiophenol

The photoinduced hydrogen elimination reaction in thiophenol via the conical intersections of the dissociative (1)pi sigma* excited state with the bound (1)pi pi* excited state and the electronic ground state has been investigated with ab initio electronic-structure calculations and time-dependent quantum wave-packet calculations. A screening of the coupling constants of the symmetry-allowed coupling modes at the (1)pi pi*-(1)pi sigma* and (1)pi sigma*-S-0 conical intersection shows that the SH torsional mode is by far the most important coupling mode at both conical intersections. A model including three intersecting potential-energy surfaces (S-0, (1)pi pi*, (1)pi sigma*) and two nuclear degrees of freedom (SH stretch and SH torsion) has been constructed on the basis of ab initio complete-active-space self-consistent field and multireference second-order perturbation theory calculations. The nonadiabatic quantum wave-packet dynamics initiated by optical excitation of the (1)pi pi* and (1)pi sigma* states has been explored for this three-state two-coordinate model. The photodissociation dynamics is characterized in terms of snapshots of time-dependent wave packets, time-dependent electronic population probabilities, and the branching ratio of the (2)sigma/(2)pi electronic states of the thiophenoxyl radical. The dependence of the timescale of the photodissociation process and the branching ratio on the initial excitation of the SH stretching and SH torsional vibrations has been analyzed. It is shown that the node structure, which is imposed on the nuclear wave packets by the initial vibrational preparation as well as by the transitions through the conical intersections, has a profound effect on the photodissociation dynamics. The effect of additional weak coupling modes of CC twist (nu(16a)) and ring-distortion (nu(16b)) character has been investigated with three-dimensional and four-dimensional time-dependent wave-packet calculations, and has been found to be minor. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4709608]

Thermal Decomposition of Propargyl Alcohol: Single Pulse Shock Tube Experimental and ab Initio Theoretical Study

Thermal decomposition of propargyl alcohol (C3H3OH), a molecule of interest in interstellar chemistry and combustion, was investigated using a single pulse shock tube in the temperature ranging from 953 to 1262 K. The products identified include acetylene, propyne, vinylacetylene, propynal, propenal, and benzene. The experimentally observed overall rate constant for thermal decomposition of propargyl alcohol was found to be k = 10((10.17 +/- 0.36)) exp(-39.70 +/- 1.83)/RT) s(-1) Ab initio theoretical calculations were carried out to understand the potential energy surfaces involved in the primary and secondary steps of propargyl alcohol thermal decomposition. Transition state theory was used to predict the rate constants, which were then used and refined in a kinetic simulation of the product profile. The first step in the decomposition is C-O bond dissociation, leading to the formation of two important radicals in combustion, OH and propargyl. This has been used to study the reverse OH propargyl radical reaction, about which there appears to be no prior work. Depending on the site of attack, this reaction leads to propargyl alcohol or propenal, one of the major products at temperatures below 1200 K. A detailed mechanism has been derived to explain all the observed products.

Study of structures, energies and vibrational frequencies of (O-2)(n)(+) (n=2-5) clusters by GGA and meta-GGA density functional methods

Using Generalized Gradient Approximation (GGA) and meta-GGA density functional methods, structures, binding energies and harmonic vibrational frequencies for the clusters O-4(+), O-6(+), O-8(+) and O-10(+) have been calculated. The stable structures of O-4(+), O-6(+), O-8(+) and O-10(+) have point groups D-2h, D-3h, D-4h, and D-5h optimized on the quartet, sextet, octet and dectet potential energy surfaces, respectively. Rectangular (D-2h) O-4(+) has been found to be more stable compared to trans-planar (C-2h) on the quartet potential energy surface. Cyclic structure (D-3h) of CA cluster ion has been calculated to be more stable than other structures. Binding energy (B.E.) of the cyclic O-6(+) is in good agreement with experimental measurement. The zero-point corrected B.E. of O-8(+) with D4h symmetry on the octet potential energy surface and zero-point corrected B.E. of O-10(+) with D-5h symmetry on the dectet potential energy surface are also in good agreement with experimental values. The B.E. value for O-4(+) is close to the experimental value when single point energy is calculated by Brueckner coupled-cluster method, BD(T). (C) 2014 Elsevier B.V. All rights reserved.

Microwave Spectrum of Hexafluoroisopropanol and Torsional Behavior of Molecules with a CF3-C-CF3 Group

This paper presents the first microwave spectroscopic investigation on hexafluoroisopropanol (HFIP). A pulsed nozzle Fourier transform microwave spectrometer has been used to determine the rotational constants for HFIP as A = 2105.12166(18) MHz, B = 1053.99503(12) MHz, and C = 932.33959(13) MHz. In addition, five isotopologues of HFIP have been observed experimentally to determine the accurate structure of HFIP. The observed spectrum could be assigned to the most stable conformer of HFIP, called antiperiplanar. Available spectroscopic information and ab initio calculations on five prototype molecules helped in exploring the torsional behavior of molecules having a CF3-C-CF3 group. Two-dimensional potential energy surfaces have been analyzed for all molecules, which explained the presence/absence of doubling in the rotational transitions. With the help of natural bond orbital (NBO) analysis, reasons for the conformational preference of HFIP have been explained.

Prediction of electronically nonadiabatic decomposition mechanisms of isolated gas phase nitrogen-rich energetic salt: Guanidium-triazolate

Electronically nonadiabatic decomposition pathways of guanidium triazolate are explored theoretically. Nonadiabatically coupled potential energy surfaces are explored at the complete active space self-consistent field (CASSCF) level of theory. For better estimation of energies complete active space second order perturbation theories (CASPT2 and CASMP2) are also employed. Density functional theory (DFT) with B3LYP functional and MP2 level of theory are used to explore subsequent ground state decomposition pathways. In comparison with all possible stable decomposition products (such as, N-2, NH3, HNC, HCN, NH2CN and CH3NC), only NH3 (with NH2CN) and N-2 are predicted to be energetically most accessible initial decomposition products. Furthermore, different conical intersections between the S-1 and S-0 surfaces, which are computed at the CASSCF(14,10)/6-31G(d) level of theory, are found to play an essential role in the excited state deactivation process of guanidium triazolate. This is the first report on the electronically nonadiabatic decomposition mechanisms of isolated guanidium triazolate salt. (C) 2015 Elsevier B.V. All rights reserved.

Structure and Dynamics of Octamethylcyclotetrasiloxane Confined between Mica Surfaces

Using a molecular model for octamethylcydotetrasiloxane (OMCTS), molecular dynamics simulations are carried out to probe the phase state of OMCTS confined between two mica surfaces in equilibrium With a reservoir. Molecular dynamics simulations are carried out for elevations ranging from 5 to 35 K above the melting point for the OMCTS model used in this study. The Helmholtz free energy is, computed for a specific confinement using the :two-phase thermodynamic (2PT) method. Analysis of the in-plane pair correlation functions did not reveal signatures of freezing even under an extreme confinement of two layers. OMCTS is found to orient with a wide distribution of orientations with respect to the mica surface, with a distinct preference for the surface parallel configuration in the contact layers. The self-intermediate scattering function is found to decay with increasing relaxation times as the surface separation is decreased, and the two-step relaxation in the scattering function, a signature of glassy dynamics, distinctly evolves as the temperature is lowered. However, even at 5 K above the melting point, we did not observe a freezing transition and the self-intermediate scattering functions relax within 200 ps for the seven-layered confined system. The self diffusivity and relaxation times obtained from the Kohlrausch-Williams-Watts stretched exponential fits to the late alpha-relaxation exhibit power law scalings with the packing fraction as predicted by mode coupling theory. A distinct discontinuity in the Helmholtz free energy, potential energy, and a sharp change in the local bond order parameter, Q(4), was observed at 230 K for a five-layered system upon cooling, indicative of a first-order transition. A freezing point depression of about 30 K was observed for this five-layered confined system, and at the lower temperatures, contact layers were found to be disordered with long-range order present only in the inner layers. These dynamical signatures indicate that confined OMCTS undergoes a slowdown akin to a fluid approaching a glass transition upon increasing confinement, and freezing under confinement would require substantial subcooling below the bulk melting point of OMCTS.

Studies involving electrostatic potential of valinomycin

The electrostatic potential of valinomycin in various conformations as obtained by the crystal structures (uncomplexed, complexed) and theoretical considerations have been evaluated and compared. The potential energy profiles along the æ axis of the bracelet-like structures show a systematic variation from the uncomplexed to the complexed structure. This type of conformational change and the potential variation are probably associated with different states of ion transport, like the capture and release of ions by the ionophore. Also, the asymmetry of the molecule due to D-HyIV on one side and L-Lac on the other side is reflected in the potential values along the Z-axis, the magnitude of which, is considerable in the uncomplexed structure. The evaluation of the potential at the ab-initio level on smaller fragments indicate that the order of liganding capacity of oxygen is amide ether ester. Also, the inductive effects due to alkyl substitution is negligible as evidenced by the potential studies on the substituted amides and esters.