Buckled nano rod - a two state system: quantum effects on its dynamics

We consider a suspended elastic rod under longitudinal compression. The compression can be used to adjust potential energy for transverse displacements from the harmonic to the double well regime. The two minima in potential energy curve describe two possible buckled states. Using transition state theory (TST) we have calculated the rate of conversion from one state to other. If the strain epsilon = 4 epsilon c the simple TST rate diverges. We suggest a method to correct this divergence for quantum calculations. We also find that zero point energy contributions can be quite large so that single mode calculations can lead to large errors in the rate.

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.

Influence of Cooling Rate on the Enthalpy Relaxation and Fragility of a Metallic Glass

Structural relaxation behavior of a rapidly quenched (RQ) and a slowly cooled Pd40Cu30Ni10P20 metallic glass was investigated and compared. Differential scanning calorimetry was employed to monitor the relaxation enthalpies at the glass transition temperature, T-g , and the Kolrausch-Williams-Watts (KWW) stretched exponential function was used to describe its variation with annealing time. It was found that the rate of enthalpy recovery is higher in the ribbon, implying that the bulk is more resistant to relaxation at low temperatures of annealing. This was attributed to the possibility of cooling rate affecting the locations where the glasses get trapped within the potential energy landscape. The RQ process traps a larger amount of free volume, resulting in higher fragility, and in turn relaxes at the slightest thermal excitation (annealing). The slowly cooled bulk metallic glass (BMG), on the other hand, entraps lower free volume and has more short-range ordering, hence requiring a large amount of perturbation to access lower energy basins.

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.

Vibrational dynamics and boson peak in a supercooled polydisperse liquid

Vibrational density of states (VDOS) in a supercooled polydisperse liquid is computed by diagonalizing the Hessian matrix evaluated at the potential energy minima for systems with different values of polydispersity. An increase in polydispersity leads to an increase in the relative population of localized high-frequency modes. At low frequencies, the density of states shows an excess compared to the Debye squared-frequency law, which has been identified with the boson peak. The height of the boson peak increases with polydispersity and shows a rather narrow sensitivity to changes in temperature. While the modes comprising the boson peak appear to be largely delocalized, there is a sharp drop in the participation ratio of the modes that exist just below the boson peak indicative of the quasilocalized nature of the low-frequency vibrations. Study of the difference spectrum at two different polydispersity reveals that the increase in the height of boson peak is due to a population shift from modes with frequencies above the maximum in the VDOS to that below the maximum, indicating an increase in the fraction of the unstable modes in the system. The latter is further supported by the facilitation of the observed dynamics by polydispersity. Since the strength of the liquid increases with polydispersity, the present result provides an evidence that the intensity of boson peak correlates positively with the strength of the liquid, as observed earlier in many experimental systems.

Fugacity coefficients. Enthalpy, and second virial coefficients of boron trifluoride

Fugacity coefficients and isothermal changes of enthalpy have been calculated and reported. The calculations cover a temperature range of 0° to 75°C. up to gas densities of 1.0 gram per cc. The generalized Benedict-Webb-Rubin constants evaluated from generalized PVT relations is found to predict the experimental data with an over-all absolute deviation of 3.1%. Second virial coefficients and potential energy parameters for Lennard-Jones (12-6) potential energy function are reported also.

Protein sequence design on the basis of topology optimization techniques -Using continuous modeling of discrete amino acid types

The notion of optimization is inherent in protein design. A long linear chain of twenty types of amino acid residues are known to fold to a 3-D conformation that minimizes the combined inter-residue energy interactions. There are two distinct protein design problems, viz. predicting the folded structure from a given sequence of amino acid monomers (folding problem) and determining a sequence for a given folded structure (inverse folding problem). These two problems have much similarity to engineering structural analysis and structural optimization problems respectively. In the folding problem, a protein chain with a given sequence folds to a conformation, called a native state, which has a unique global minimum energy value when compared to all other unfolded conformations. This involves a search in the conformation space. This is somewhat akin to the principle of minimum potential energy that determines the deformed static equilibrium configuration of an elastic structure of given topology, shape, and size that is subjected to certain boundary conditions. In the inverse-folding problem, one has to design a sequence with some objectives (having a specific feature of the folded structure, docking with another protein, etc.) and constraints (sequence being fixed in some portion, a particular composition of amino acid types, etc.) while obtaining a sequence that would fold to the desired conformation satisfying the criteria of folding. This requires a search in the sequence space. This is similar to structural optimization in the design-variable space wherein a certain feature of structural response is optimized subject to some constraints while satisfying the governing static or dynamic equilibrium equations. Based on this similarity, in this work we apply the topology optimization methods to protein design, discuss modeling issues and present some initial results.

Infrared Spectra of Dimethylquinolines in the Gas Phase: Experiment and Theory

Infrared spectra of atmospherically important dimethylquinolines (DMQs), namely 2,4-DMQ, 2,6-DMQ, 2,7-DMQ, and 2,8-DMQ in the gas phase at 80 degrees C were recorded using a long variable path-length cell. DFT calculations were carried out to assign the bands in the experimentally observed spectra at the B3LYP/6-31G* level of theory. The spectral assignments particularly for the C-H stretching modes could not be made unambiguously using calculated anharmonic or scaled harmonic frequencies. To resolve this problem, a scaled force field method of assignment was used. Assignment of fundamental modes was confirmed by potential energy distributions (PEDs) of the normal modes derived by the scaled force fields using a modified version of the UMAT program in the QCPE package. We demonstrate that for large molecules such as the DMQs, the scaling of the force field is more effective in arriving at the correct assignment of the fundamentals for a quantitative vibrational analysis. An error analysis of the mean deviation of the calculated harmonic, anharmonic, and force field fitted frequencies from the observed frequency provides strong evidence for the correctness of the assignment.

A rectangular laminated anisotropic shallow thin shell finite element

This report contains the details of the development of the stiffness matrix for a rectangular laminated anisotropic shallow thin shell finite element. The derivation is done under linear thin shell assumptions. Expressing the assumed displacement state over the middle surface of the shell as products of one-dimensional first-order Hermite interpolation polynomials, it is possible to insure that the displacement state for the assembled set of such elements, to be geometrically admissible. Monotonic convergence of the total potential energy is therefore possible as the modelling is successively refined. The element is systematically evaluated for its performance considering various examples for which analytical or other solutions are available

Investigation of the effect of nonlocal scale on ultrasonic wave dispersion characteristics of a monolayer graphene

In this article, an ultrasonic wave propagation in graphene sheet is studied using nonlocal elasticity theory incorporating small scale effects. The graphene sheet is modeled as an isotropic plate of one-atom thick. For this model, the nonlocal governing differential equations of motion are derived from the minimization of the total potential energy of the entire system. An ultrasonic type of wave propagation model is also derived for the graphene sheet. The nonlocal scale parameter introduces certain band gap region in in-plane and flexural wave modes where no wave propagation occurs. This is manifested in the wavenumber plots as the region where the wavenumber tends to infinite or wave speed tends to zero. The frequency at which this phenomenon occurs is called the escape frequency. The explicit expressions for cutoff frequencies and escape frequencies are derived. The escape frequencies are mainly introduced because of the nonlocal elasticity. Obviously these frequencies are function of nonlocal scaling parameter. It has also been obtained that these frequencies are independent of y-directional wavenumber. It means that for any type of nanostructure, the escape frequencies are purely a function of nonlocal scaling parameter only. It is also independent of the geometry of the structure. It has been found that the cutoff frequencies are function of nonlocal scaling parameter (e(0)a) and the y-directional wavenumber (k(y)). For a given nanostructure, nonlocal small scale coefficient can be obtained by matching the results from molecular dynamics (MD) simulations and the nonlocal elasticity calculations. At that value of the nonlocal scale coefficient, the waves will propagate in the nanostructure at that cut-off frequency. In the present paper, different values of e(o)a are used. One can get the exact e(0)a for a given graphene sheet by matching the MD simulation results of graphene with the results presented in this paper. (C) 2010 Elsevier B.V. All rights reserved.

Theoretical studies on alpha-helix--DNA interactions

The interaction energies between (Ala)10 and alpha-helix fragment and different nucleotide sequences in right-handed B-form have been optimized using semi-empirical potential energy functions. The energies are calculated for two different orientations of the alpha-helix, viz., when the alpha-helix axis taken in the N----C direction is (i) parallel and (ii) antiparallel to the 5'-3' ascending strand of DNA, proximal to it. When both the DNA molecule as well as the alpha-helix are treated as rigid molecules it is found that a polyalanine alpha-helix has slightly more favourable contacts when it is in the proximity of a four nucleotide sequence of 5'-(N-A-T-N)-3' type, where N is either a purine or a pyrimidine. However, when the two interacting molecules are allowed to undergo local structural variations then the interaction energy appears to be independent of the base sequence confirming the non-specific nature of these interactions.

A molecular dynamics study of cage-to-cage migration in sodium Y zeolite: role of surface-mediated diffusion

The details of cage-to-cage migration have been obtained from an analysis of the molecular dynamics trajectory of a probe adsorbate. It is observed that particles utilize the region within a radius of 2 angstrom from the window center but with diffusion taking place predominantly at 1.6 angstrom from the window center and a potential energy of nearly -12 kJ/mol. A barrier of about 0.5 kJ/mol is observed for surface-mediated diffusion. Surprisingly, for diffusion without surface mediation for a particle going from one cage center to another, there is an attractive well near the window instead of a barrier. At low adsorbate concentrations and room temperature, the predominant mode for cage-to-cage migration is surface-mediated diffusion. The analysis suggests that particles slide along the surface of the inner walls of the alpha-cages during migration from one cage to another.

Vibrational analysis of 1,3-thiazolidine-2-thione and its selenium analogue

The i.r. spectra of 1,3-thiazolidine-2-thione, ?2-selone and their N-deuteriated derivatives have been investigated in the region between 4000 and 20 cm?1. A complete assignment of the fundamental vibrational frequencies has been made based upon the normal coordinate analysis carried out using a simple Urey�Bradley force field supplemented by the valence force constants for the out-of-plane vibrations of the planar molecular skeleton. The proposed assignments are discussed in relation to the group frequencies in structurally related molecules and in terms of the computed potential energy distributions among the symmetry coordinates.

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.