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.

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.

Temperature dependent free energy surface of polymer folding from equilibrium and quench studies

Langevin dynamics simulation studies have been employed to calculate the temperature dependent free energy surface and folding characteristics of a 500 monomer long linear alkane (polyethylene) chain with a realistic interaction potential. Both equilibrium and temperature quench simulation studies have been carried out. Using the shape anisotropy parameter (S) of the folded molecule as the order parameter, we find a weakly first order phase transition between the high-temperature molten globule and low-temperature rodlike crystalline states separated by a small barrier of the order of k(B)T. Near the melting temperature (580 K), we observe an intriguing intermittent fluctuation with pronounced ``1/f noise characteristics'' between these two states with large difference in shape and structure. We have also studied the possibilities of different pathways of folding to states much below the melting point. At 300 K starting from the all-trans linear configuration, the chain folds stepwise into a very regular fourfold crystallite with very high shape anisotropy. Whereas, when quenched from a high temperature (900 K) random coil regime, we identify a two step transition from the random coiled state to a molten globulelike state and, further, to a anisotropic rodlike state. The trajectory reveals an interesting coupling between the two order parameters, namely, radius of gyration (R-g) and the shape anisotropy parameter (S). The rodlike final state of the quench trajectory is characterized by lower shape anisotropy parameter and significantly larger number of gauche defects as compared to the final state obtained through equilibrium simulation starting from all-trans linear chain. The quench study shows indication of a nucleationlike pathway from the molten globule to the rodlike state involving an underlying rugged energy landscape. (C) 2010 American Institute of Physics. doi:10.1063/1.3509398]

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.

Density Functional Theoretical Modeling, Electrostatic Surface Potential and Surface Enhanced Raman Spectroscopic Studies on Biosynthesized Silver Nanoparticles: Observation of 400 pM Sensitivity to Explosives

Interaction of adsorbate on charged surfaces, orientation of the analyte on the surface, and surface enhancement aspects have been studied. These aspects have been explored in details to explain the surface-enhanced Raman spectroscopic (SERS) spectra of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (HNIW or CL-20), a well-known explosive, and 2,4,6-trinitrotoluene (TNT) using one-pot synthesis of silver nanoparticles via biosynthetic route using natural precursor extracts of clove and pepper. The biosynthesized silver nanoparticles (bio Ag Nps) have been characterized using UV-vis spectroscopy, scanning electron microscopy and atomic force microscopy. SERS studies conducted using bio Ag Nps on different water insoluble analytes, such as CL-20 and TNT, lead to SERS signals at concentration levels of 400 pM. The experimental findings have been corroborated with density functional computational results, electrostatic surface potential calculations, Fukui functions and potential measurements.

Multidimensional free energy surface of unfolding of HP-36: Microscopic origin of ruggedness

The protein folding funnel paradigm suggests that folding and unfolding proceed as directed diffusion in a multidimensional free energy surface where a multitude of pathways can be traversed during the protein's sojourn from initial to final state. However, finding even a single pathway, with the detail chronicling of intermediates, is an arduous task. In this work we explore the free energy surface of unfolding pathway through umbrella sampling, for a small globular a-helical protein chicken-villin headpiece (HP-36) when the melting of secondary structures is induced by adding DMSO in aqueous solution. We find that the unfolding proceeds through the initial separation or melting of aggregated hydrophobic core that comprises of three phenylalanine residues (Phe7, Phe11, and Phe18). This separation is accompanied by simultaneous melting of the second helix. Unfolding is found to be a multistage process involving crossing of three consecutive minima and two barriers at the initial stage. At a molecular level, Phe18 is observed to reorient itself towards other hydrophobic grooves to stabilize the intermediate states. We identify the configuration of the intermediates and correlate the intermediates with those obtained in our previous works. We also give an estimate of the barriers for different transition states and observe the softening of the barriers with increasing DMSO concentration. We show that higher concentration of DMSO tunes the unfolding pathway by destabilizing the third minimum and stabilizing the second one, indicating the development of a solvent modified, less rugged pathway. The prime outcome of this work is the demonstration that mixed solvents can profoundly transform the nature of the energy landscape and induce unfolding via a modified route. A successful application of Kramer's rate equation correlating the free energy simulation results shows faster rate of unfolding with increasing DMSO concentration. This work perhaps presents the first systematic theoretical study of the effect of a chemical denaturant on the microscopic free energy surface and rates of unfolding of HP-36. (C) 2014 AIP Publishing LLC.

Towards an exact factorization of the molecular wave function

An exact single-product factorisation of the molecular wave function for the timedependent Schrodinger equation is investigated by using an ansatz involving a phasefactor. By using the Frenkel variational method, we obtain the Schrodinger equations for the electronic and nuclear wave functions. The concept of a potential energy surface (PES) is retained by introducing a modified Hamiltonian as suggested earlier by Cederbaum. The parameter in the phase factor is chosen such that the equations of motion retain the physically appealing Born- Oppenheimer-like form, and is therefore unique.

Relationship between entropy and diffusion: A statistical mechanical derivation of Rosenfeld expression for a rugged energy landscape

Diffusion-a measure of dynamics, and entropy-a measure of disorder in the system are found to be intimately correlated in many systems, and the correlation is often strongly non-linear. We explore the origin of this complex dependence by studying diffusion of a point Brownian particle on a model potential energy surface characterized by ruggedness. If we assume that the ruggedness has a Gaussian distribution, then for this model, one can obtain the excess entropy exactly for any dimension. By using the expression for the mean first passage time, we present a statistical mechanical derivation of the well-known and well-tested scaling relation proposed by Rosenfeld between diffusion and excess entropy. In anticipation that Rosenfeld diffusion-entropy scaling (RDES) relation may continue to be valid in higher dimensions (where the mean first passage time approach is not available), we carry out an effective medium approximation (EMA) based analysis of the effective transition rate and hence of the effective diffusion coefficient. We show that the EMA expression can be used to derive the RDES scaling relation for any dimension higher than unity. However, RDES is shown to break down in the presence of spatial correlation among the energy landscape values. (C) 2015 AIP Publishing LLC.

Dynamics of Barrierless and Activated Chemical Reactions in a Dispersive Medium within the Fractional Diffusion Equation Approach

Barrierless chemical reactions have often been modeled as a Brownian motion on a one-dimensional harmonic potential energy surface with a position-dependent reaction sink or window located near the minimum of the surface. This simple (but highly successful) description leads to a nonexponential survival probability only at small to intermediate times but exponential decay in the long-time limit. However, in several reactive events involving proteins and glasses, the reactions are found to exhibit a strongly nonexponential (power law) decay kinetics even in the long time. In order to address such reactions, here, we introduce a model of barrierless chemical reaction where the motion along the reaction coordinate sustains dispersive diffusion. A complete analytical solution of the model can be obtained only in the frequency domain, but an asymptotic solution is obtained in the limit of long time. In this case, the asymptotic long-time decay of the survival probability is a power law of the Mittag−Leffler functional form. When the barrier height is increased, the decay of the survival probability still remains nonexponential, in contrast to the ordinary Brownian motion case where the rate is given by the Smoluchowski limit of the well-known Kramers' expression. Interestingly, the reaction under dispersive diffusion is shown to exhibit strong dependence on the initial state of the system, thus predicting a strong dependence on the excitation wavelength for photoisomerization reactions in a dispersive medium. The theory also predicts a fractional viscosity dependence of the rate, which is often observed in the reactions occurring in complex environments.

Structure and Bonding in Cyclic Isomers of B2AlHnm (n = 3−6, m = −2 to +1): A Comparative Study with B3Hnm, BAl2Hnm and Al3Hnm†

The structure, bonding and energetics of B2AlHnm (n = 3−6, m = −2 to +1) are compared with corresponding homocyclic boron, aluminum analogues and BAl2Hnm using density functional theory (DFT). Divalent to hexacoordinated boron and aluminum atoms are found in these species. The geometrical and bonding pattern in B2AlH4− is similar to that for B2SiH4. Species with lone pairs on the divalent boron and aluminum atoms are found to be minima on the potential energy surface of B2AlH32−. A dramatic structural diversity is observed in going from B3Hnm to B2AlHnm, BAl2Hnm and Al3Hnm and this is attributable to the preference of lower coordination on aluminum, higher coordination on boron and the higher multicenter bonding capability of boron. The most stable structures of B3H6+, B2AlH5 and BAl2H4− and the trihydrogen bridged structure of Al3H32− show an isostructural relationship, indicating the isolobal analogy between trivalent boron and divalent aluminum anion.

Structure and bonding of MCB5H7 and its sandwiched dimer CB5H6M–MCB5H6 (M = Si, Ge, Sn): Isomer stability and preference for slip distorted structure

We find sandwiched metal dimers CB5H6M–MCB5H6 (M = Si, Ge, Sn) which are minima in the potential energy surface with a characteristic M–M single bond. The NBO analysis and the M–M distances (Å) (2.3, 2.44 and 2.81 for M = Si, Ge, Sn) indicate substantial M–M bonding. Formal generation of CB5H6M–MCB5H6 has been studied theoretically. Consecutive substitution of two boron atoms in B7H−27 by M (Si, Ge, Sn) and carbon, respectively followed by dehydrogenation may lead to our desired CB5H6M–MCB5H6. We find that the slip distorted geometry is preferred for MCB5H7 and its dehydrogenated dimer CB5H6M–MCB5H6. The slip-distortion of M–M bond in CB5H6M–MCB5H6 is more than the slip distortion of M–H bond in MCB5H7. Molecular orbital analysis has been done to understand the slip distortion. Larger M–M bending (CB5H6M–MCB5H6) in comparison with M–H bending (MCB5H7) is suspected to be encouraged by stabilization of one of the M–M π bonding MO’s. Preference of M to occupy the apex of pentagonal skeleton of MCB5H7 over its icosahedral analogue MCB10H11 has been observed.

Ab initio molecular orbital calculations on ion-molecule and ion pair-molecule complexes of the water-lithium cyanide system

Ab initio molecular orbital (MO) calculations with the 3-21G and 6-31G basis sets were performed on a series of ion-molecule and ion pair-molecule complexes for the H2O + LiCN system. Stabilisation energies (with counter-poise corrections), geometrical parameters, internal force constants and harmonic vibrational frequencies were evaluated for 16 structures of interest. Although the interaction energies are smaller, the geometries and relative stabilities of the monohydrated contact ion pair are reminiscent of those computed for the complexes of the individual ions. Thus, interaction of the oxygen lone pair with lithium leads to a highly stabilised C2v structure, while the coordination of water to the cyanide ion involves a slightly non-linear hydrogen bond. Symmetrical bifurcated structures are computed to be saddle points on the potential energy surface, and to have an imaginary frequency for the rocking mode of the water molecule. On optimisation the geometries of the solvent shared ion pair structures (e.g. Li+cdots, three dots, centered OH2cdots, three dots, centered CN−) revealed a proton transfer from the water molecule leading to hydrogen bonded forms such as Li-O-Hcdots, three dots, centered HCN. The variation in the force constants and harmonic frequencies in the various structures considered are discussed in terms of ion-molecular and ion pair-molecule interactions.

Effect of Sorbate-Zeolite Interaction on Cluster Lifetime and Size of Sorbates: Xe in NaY

Molecular dynamics calculations are reported for Xe in sodium Y zeolite with varying strengths of sorbate-zeolite dispersion interaction. In the absence of any dispersion interaction between the sorbate and the zeolite, the presence of the zeolite has a purely geometrical role. Increase in the strength of the sorbate-zeolite interaction increases the monomer population and decreases the population of dimers and higher sized clusters. The lifetime of the monomers as well as dimers increases with the strength of the dispersion interaction. The observed variations in the lifetime and the population of the different sized clusters is explained in terms of the changes in the potential energy surface caused by the increase in the strength of the dispersion interaction.