Accurate ab initio potential energy surfaces for the (3)A '' and (3)A ' electronic states of the O(P-3) plus HBr system

In this work, we report the construction of potential energy surfaces for the (3)A '' and (3)A' states of the system O(P-3) + HBr. These surfaces are based on extensive ab initio calculations employing the MRCI+Q/CBS+SO level of theory. The complete basis set energies were estimated from extrapolation of MRCI+Q/aug-cc-VnZ(-PP) (n = Q, 5) results and corrections due to spin-orbit effects obtained at the CASSCF/aug-cc-pVTZ(-PP) level of theory. These energies, calculated over a region of the configuration space relevant to the study of the reaction O(P-3) + HBr -> OH + Br, were used to generate functions based on the many-body expansion. The three-body potentials were interpolated using the reproducing kernel Hilbert space method. The resulting surface for the (3)A '' electronic state contains van der Waals minima on the entrance and exit channels and a transition state 6.55 kcal/mol higher than the reactants. This barrier height was then scaled to reproduce the value of 5.01 kcal/mol, which was estimated from coupled cluster benchmark calculations performed to include high-order and core-valence correlation, as well as scalar relativistic effects. The (3)A' surface was also scaled, based on the fact that in the collinear saddle point geometry these two electronic states are degenerate. The vibrationally adiabatic barrier heights are 3.44 kcal/mol for the (3)A '' and 4.16 kcal/mol for the (3)A' state. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4705428]

Internal oscillating current-sustained RF plasmas : Parameters, stability, and potential for surface engineering

A new source of low-frequency (0.46 MHz) inductively coupled plasmas sustained by the internal planar "unidirectional" RF current driven through a specially designed internal antenna configuration has been developed. The experimental results of the investigation of the optical and global argon plasma parameters by the optical and Langmuir probes are presented. It is shown that the spatial profiles of the electron density, the effective electron temperature and plasma potential feature a great deal of the radial and axial uniformity compared with conventional sources of inductively coupled plasmas with external at coil configurations. The measurements also reveal a weak azimuthal dependence of the global plasma parameters at low values of the input RF power, which was earlier predicted theoretically. The azimuthal dependence of the global plasma parameters vanishes at high input RF powers. Moreover, under certain conditions, the plasma becomes unstable due to spontaneous transitions between low-density (electrostatic, E) and high-density (electromagnetic, H) operating modes. Excellent uniformity of high-density plasmas makes the plasma reactor promising for various plasma processing applications and surface engineering.

Potential-dependent surface-enhanced resonance raman spectroscopy at nanostructured TiO2 : A case study on cytochrome b5

Abstract: Nanostructured titanium dioxide (TiO2) electrodes, prepared by anodization of titanium, are employed to probe the electron-transfer process of cytochrome b5 (cyt b5) by surface-enhanced resonance Raman (SERR) spectroscopy. Concomitant with the increased nanoscopic surface roughness of TiO2, achieved by raising the anodization voltage from 10 to 20 V, the enhancement factor increases from 2.4 to 8.6, which is rationalized by calculations of the electric field enhancement. Cyt b 5 is immobilized on TiO2 under preservation of its native structure but it displays a non-ideal redox behavior due to the limited conductivity of the electrode material. The electron-transfer efficiency which depends on the crystalline phase of TiO2 has to be improved by appropriate doping for applications in bioelectrochemistry. Nanostructured TiO2 electrodes are employed to probe the electron-transfer process of cytochrome b5 by surface-enhanced resonance Raman spectroscopy. Concomitant with the increased nanoscopic surface roughness of TiO2, the enhancement factor increases, which can be attributed to the electric field enhancement. Copyright © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

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.

Potential energy model for turbulent entrainment across density interface in a linearly stratified fluid

A potential energy model is developed for turbulent entrainment in the absence of mean shear in a linearly stratified fluid. The relation between the entrainment distance D and the time t and the relation between dimensionless entrainment rate E and the local Richardson number are obtained. An experiment is made for examination. The experimental results are in good agreement with the model, in which the dimensionless entrainment distance D is given by DBAR = A(i)(SBAR)-1/4(fBAR)1/2(tBAR)1/8, where A(i) is the proportional coefficient, S is the dimensionless stroke, fBAR is the dimensionless frequency of the grid oscillation, tBAR the dimensionless time.

Viscosity-induced mode splitting and potential energy criterion for mode stability

This paper points out that viscosity can induce mode splitting in a uniform infinite cylinder of an incompressible fluid with self-gravitation, and that the potential energy criterion cannot be appropriate to all normal modes obtained, i.e., there will be stable modes with negative potential energy (<0). Therefore the condition >0 is not necessary, although sufficient, for the stability of a mode in an incompressible static fluid or magnetohydrodynamics (MHD) system, which is a correction of both Hare's [Philos. Mag. 8, 1305 (1959)] and Chandrasekhar's [Hydrodynamic and Hydromagnetic Stability (Oxford U.P., Oxford, 1961), p. 604] stability criterion for a mode. These results can also be extended to compressible systems with a polytropic exponent.

Feeding on protists and particulates by the leptocephali of the worm eels Myrophis spp. (Teleostei: Anguilliformes: Ophichthidae), and the potential energy contribution of large aloricate protozoa

The food sources of the leptocephali of the teleostean superorder Elopomorpha have been controversial, yet observations on the leptocephali of the worm eels, Myrophis spp. (family Ophichthidae) collected in the northern Gulf of Mexico indicate active, not passive, feeding. Leptocephali had protists in their alimentary canals. Estimates of the physiological energetics of worm eels indicate that large aloricate protozoa including ciliates could provide substantial energy to these leptocephali toward the end of the premetamorphic and metamorphic stages, given the low energy requirements of metamorphosing leptocephali. Global ocean warming will likely force a shift in oceanic food webs; a shift away from large protozoa toward smaller protists is possible. Such a disruption of the oceanic food webs could further compromise the survival of leptocephali.

Gaussian Approximation Potential: an interatomic potential derived from first principles Quantum Mechanics

Simulation of materials at the atomistic level is an important tool in studying microscopic structure and processes. The atomic interactions necessary for the simulation are correctly described by Quantum Mechanics. However, the computational resources required to solve the quantum mechanical equations limits the use of Quantum Mechanics at most to a few hundreds of atoms and only to a small fraction of the available configurational space. This thesis presents the results of my research on the development of a new interatomic potential generation scheme, which we refer to as Gaussian Approximation Potentials. In our framework, the quantum mechanical potential energy surface is interpolated between a set of predetermined values at different points in atomic configurational space by a non-linear, non-parametric regression method, the Gaussian Process. To perform the fitting, we represent the atomic environments by the bispectrum, which is invariant to permutations of the atoms in the neighbourhood and to global rotations. The result is a general scheme, that allows one to generate interatomic potentials based on arbitrary quantum mechanical data. We built a series of Gaussian Approximation Potentials using data obtained from Density Functional Theory and tested the capabilities of the method. We showed that our models reproduce the quantum mechanical potential energy surface remarkably well for the group IV semiconductors, iron and gallium nitride. Our potentials, while maintaining quantum mechanical accuracy, are several orders of magnitude faster than Quantum Mechanical methods.

Free energy surface reconstruction from umbrella samples using Gaussian process regression

We demonstrate how the Gaussian process regression approach can be used to efficiently reconstruct free energy surfaces from umbrella sampling simulations. By making a prior assumption of smoothness and taking account of the sampling noise in a consistent fashion, we achieve a significant improvement in accuracy over the state of the art in two or more dimensions or, equivalently, a significant cost reduction to obtain the free energy surface within a prescribed tolerance in both regimes of spatially sparse data and short sampling trajectories. Stemming from its Bayesian interpretation the method provides meaningful error bars without significant additional computation. A software implementation is made available on www.libatoms.org.