Fractal and complex network analyses of protein molecular dynamics

Based on protein molecular dynamics, we investigate the fractal properties of energy, pressure and volume time series using the multifractal detrended fluctuation analysis (MF-DFA) and the topological and fractal properties of their converted horizontal visibility graphs (HVGs). The energy parameters of protein dynamics we considered are bonded potential, angle potential, dihedral potential, improper potential, kinetic energy, Van der Waals potential, electrostatic potential, total energy and potential energy. The shape of the h(q)h(q) curves from MF-DFA indicates that these time series are multifractal. The numerical values of the exponent h(2)h(2) of MF-DFA show that the series of total energy and potential energy are non-stationary and anti-persistent; the other time series are stationary and persistent apart from series of pressure (with H≈0.5H≈0.5 indicating the absence of long-range correlation). The degree distributions of their converted HVGs show that these networks are exponential. The results of fractal analysis show that fractality exists in these converted HVGs. For each energy, pressure or volume parameter, it is found that the values of h(2)h(2) of MF-DFA on the time series, exponent λλ of the exponential degree distribution and fractal dimension dBdB of their converted HVGs do not change much for different proteins (indicating some universality). We also found that after taking average over all proteins, there is a linear relationship between 〈h(2)〉〈h(2)〉 (from MF-DFA on time series) and 〈dB〉〈dB〉 of the converted HVGs for different energy, pressure and volume.

The Broad-Line Type Ic Supernoava SN 2007ru: Adding To The Diversity Of Type Ic Supernovae

Photometric and spectral evolution of the Type Ic supernova SN 2007ru until around 210 days after maximum are presented. The spectra show broad spectral features due to very high expansion velocity, normally seen in hypernovae. The photospheric velocity is higher than other normal Type Ic supernovae (SNe Ic). It is lower than SN 1998bw at similar to 8 days after the explosion, but is comparable at later epochs. The light curve (LC) evolution of SN 2007ru indicates a fast rise time of 8 +/- 3 days to B-band maximum and postmaximum decline more rapid than other broad-line SNe Ic. With an absolute V magnitude of -19.06, SN 2007ru is comparable in brightness with SN 1998bw and lies at the brighter end of the observed SNe Ic. The ejected mass of Ni-56 is estimated to be similar to 0.4 M-circle dot. The fast rise and decline of the LC and the high expansion velocity suggest that SN 2007ru is an explosion with a high kinetic energy/ejecta mass ratio (E-K/M-ej). This adds to the diversity of SNe Ic. Although the early phase spectra are most similar to those of broad-line SN 2003jd, the [O I] line profile in the nebular spectrum of SN 2007ru shows the singly peaked profile, in contrast to the doubly peaked profile in SN 2003jd. The singly peaked profile, together with the high luminosity and the high expansion velocity, may suggest that SN 2007ru could be an aspherical explosion viewed from the polar direction. Estimated oxygen abundance 12 + log(O/H) of similar to 8.8 indicates that SN 2007ru occurred in a region with nearly solar metallicity.

Power system stability implications from electromechanical wave propagation

Electromechanical wave propagation characterizes the first-swing dynamic response in a spatially delayed manner. This paper investigates the characteristics of this phenomenon in two-dimensional and one-dimensional power systems. In 2-D systems, the wave front expands as a ripple in a pond. In 1-D systems, the wave front is more concentrated, retains most of its magnitude, and travels like a pulse on a string. This large wave front is more impactful upon any weak link and easily causes transient instability in 1-D systems. The initial disturbance injects both high and low frequency components, but the lumped nature of realistic systems only permits the lower frequency components to propagate through. The kinetic energy split at a junction is equal to the generator inertia ratio in each branch in an idealized continuum system. This prediction is approximately valid in a realistic power system. These insights can enhance understanding and control of the traveling waves.

An investigation of bond formation in the weakly bound first excited 1Σ and lowest 3Σ states of HeH+

The role of the electronic kinetic energy and its Cartesian components is examined during the formation of the first excited 1�£ and the lowest 3�£ states of HeH+ employing wavefunctions of multi-configuration type with basis orbitals in elliptic coordinates. Results show that the bond formation in these states is preceded primarily by a charge transfer from H to He+ rather than by polarisation of the H-orbital by He+

Ab initio molecular dynamics simulations of excited hydrogen halides and methyl halides

Using excited-state ab initio molecular dynamics simulations employing the complete-active-space self-consistent-field approach, we study the mechanism of photodissociation in terms of time evolution of structure, kinetic energy, charges and potential energy for the first excited state of hydrogen halides and methyl halides. Although the hydrogen halides and methyl halides are similar in the photodissociation mechanism, their dynamics are slightly different. The presence of the methyl group causes delay in photodissociation as compared to hydrogen halides.

Lifting the singular nature of a model for peeling of an adhesive tape

We investigate the dynamics of peeling of an adhesive tape subjected to a constant pull speed. Due to the constraint between the pull force, peel angle and the peel force, the equations of motion derived earlier fall into the category of differential-algebraic equations (DAE) requiring an appropriate algorithm for its numerical solution. By including the kinetic energy arising from the stretched part of the tape in the Lagrangian, we derive equations of motion that support stick-slip jumps as a natural consequence of the inherent dynamics itself, thus circumventing the need to use any special algorithm. In the low mass limit, these equations reproduce solutions obtained using a differential-algebraic algorithm introduced for the earlier singular equations. We find that mass has a strong influence on the dynamics of the model rendering periodic solutions to chaotic and vice versa. Apart from the rich dynamics, the model reproduces several qualitative features of the different waveforms of the peel force function as also the decreasing nature of force drop magnitudes.

Structural studies of decaying fluid turbulence: Effect of initial conditions

We present results from a systematic numerical study of structural properties of an unforced, incompressible, homogeneous, and isotropic three-dimensional turbulent fluid with an initial energy spectrum that develops a cascade of kinetic energy to large wave numbers. The results are compared with those from a recently studied set of power-law initial energy spectra [C. Kalelkar and R. Pandit, Phys. Rev. E 69, 046304 (2004)] which do not exhibit such a cascade. Differences are exhibited in plots of vorticity isosurfaces, the temporal evolution of the kinetic energy-dissipation rate, and the rates of production of the mean enstrophy along the principal axes of the strain-rate tensor. A crossover between "non-cascade-type" and "cascade-type" behavior is shown numerically for a specific set of initial energy spectra.

Type I seesaw mechanism for quasi-degenerate neutrinos

We discuss symmetries and scenarios leading to quasi-degenerate neutrinos in type I seesaw models. The existence of degeneracy in the present approach is not linked to any specific structure for the Dirac neutrino Yukawa coupling matrix y(D) and holds in general. Basic input is the application of the minimal flavour violation principle to the leptonic sector. Generalizing this principle, we assume that the structure of the right-handed neutrino mass matrix is determined by y(D) and the charged lepton Yukawa coupling matrix y(l) in an effective theory invariant under specific groups G(F) contained in the full symmetry group of the kinetic energy terms. G(F) invariance also leads to specific structure for the departure from degeneracy. The neutrino mass matrix (with degenerate mass m(0)) resulting after seesaw mechanism has a simple form Mv approximate to m(0)(I - py(l)y(l)(T)) in one particular scenario based on supersymmetry. This form is shown tolead to correct description of neutrino masses and mixing angles. The thermal leptogenesis after inclusion of flavour effects can account for the observed baryon asymmetry of the universe within the present scenario. Rates for lepton flavour violating processes can occur at observable levels in the supersymmetric version of the scenario. (c) 2010 Elsevier B.V. All rights reserved.

On the alleged equivalence of certain field theories

Abstract. We critically examine some recent claims that certain field theories with and without boson kinetic energy terms are equivalent. We point out that the crucial element in these claims is the finiteness or otherwise of the boson wavefunction renormalisation constant. We show that when this constant is finite, the equivalence proof offered in the literature fails in a direct way. When the constant is divergent, the claimed equivalence is only a consequence of improper use of divergent quantities.

Thermal gap conductance at low contact pressures (<1 MPa): Effect of gold plating and plating thickness

Thermal contact conductance (TCC) measurements are made on bare and gold plated (<= 0.5 mu m) oxygen free high conductivity (OFHC) Cu and brass contacts in vacuum, nitrogen, and argon environments. It is observed that the TCC in gaseous environment is significantly higher than that in vacuum due to the enhanced thermal gap conductance. It is found that for a given contact load and gas pressure, the thermal gap conductance for bare OFHC Cu contacts is higher than that for gold plated contacts. It is due to the difference in the molecular weights of copper and gold, which influences the exchange of kinetic energy between the gas molecules and contact surfaces. Furthermore, the gap conductance is found to increase with increasing thickness of gold plating. Topography measurements and scanning electron microscopy (SEM) analysis of contact surfaces revealed that surfaces become smoother with increasing gold plating thickness, thus resulting in smaller gaps and consequently higher gap conductance. (C) 2010 Elsevier Ltd. All rights reserved.

Systematic Organic UV Photoelectron Spectroscopy

Photoelectron spectroscopy (PES) provides valuable information on the ionization energies of atoms and molecules. The ionization energy (IE) is given by the relation.hv = IE + T where hv is t h e energy of the radiation and T i s the kinetic energy of the electron. The IEs are directly related to the orbital energies (Koopmans' theorem). By employing UV radiation (HeI. 21.2 eV. or HeII. 40.8 eV). extensive data on the ionization of valence electrons in organic molecules have been obtained in recent years. These studies of UV photoelectron spectroscopy. originated by Turner, have provided a direct probe into the energy levels of organic molecules. Molecular orbital calculations of various degrees of sophistication are generally employed to make assignments of the PES bands. Analysis of the vibrational structure of PES bands has not only provided structural information on the molecular ions, but has also been of value in band assignments. Dewar and co-workers [1, 2) presented summaries of available PES data on organic molecules in 1969 and 1970. Turner et al. [3] published a handbook of Hel spectra of organic molecules in 1970. Since then, a few books [4-7] discussing the principles and applications of UV photoelectron spectroscopy have appeared of which special mention should be made of the recent article by Heilbronner and Maier [7]. There has, however, been no comprehensive review of the vast amount of data on the UV-PES of organic molecules published in the literature since 1970.

Metal-Insulator Transition in the Hubbard Model on a Triangular Lattice

We discuss the results of an extensive mean-field investigation of the half-filled Hubbard model on a triangular lattice at zero temperature. At intermediate U we find a first-order metal-insulator transition from an incommensurate spiral magnetic metal to a semiconducting state with a commensurate linear spin density wave ordering stabilized by the competition between the kinetic energy and the frustrated nature of the magnetic interaction. At large U the ground state is that of a classical triangular antiferromagnet within our approximation. In the incommensurate spiral metallic phase the Fermi surface has parts in which the wave function renormalization Z is extremely small. The evolution of the Fermi surface and the broadening of the quasi-particle band along with the variation of the plasma frequency and a charge stiffness constant with U/t are discussed.

Paradoxical reflection in quantum mechanics

This article concerns a phenomenon of elementary quantum mechanics that is quite counter-intuitive, very non-classical, and apparently not widely known: a quantum particle can get reflected at a downward potential step. In contrast, classical particles get reflected only at upward steps. The conditions for this effect are that the wave length is much greater than the width of the potential step and the kinetic energy of the particle is much smaller than the depth of the potential step. This phenomenon is suggested by non-normalizable solutions to the time-independent Schroedinger equation, and we present evidence, numerical and mathematical, that it is also indeed predicted by the time-dependent Schroedinger equation. Furthermore, this paradoxical reflection effect suggests, and we confirm mathematically, that a quantum particle can be trapped for a long time (though not forever) in a region surrounded by downward potential steps, that is, on a plateau.

Inherent structures of phase-separating binary mixtures: Nucleation, spinodal decomposition, and pattern formation

An energy landscape view of phase separation and nonideality in binary mixtures is developed by exploring their potential energy landscape (PEL) as functions of temperature and composition. We employ molecular dynamics simulations to study a model that promotes structure breaking in the solute-solvent parent binary liquid, at low temperatures. The PEL of the system captures the potential energy distribution of the inherent structures (IS) of the system and is obtained by removing the kinetic energy (including that of intermolecular vibrations). The broader distribution of the inherent structure energy for structure breaking liquid than that of the structure making liquid demonstrates the larger role of entropy in stabilizing the parent liquid of the structure breaking type of binary mixtures. At high temperature, although the parent structure of the structure breaking binary mixture is homogenous, the corresponding inherent structure is found to be always phase separated, with a density pattern that exhibits marked correlation with the energy of its inherent structure. Over a broad range of intermediate inherent structure energy, bicontinuous phase separation prevails with interpenetrating stripes as signatures of spinodal decomposition. At low inherent structure energy, the structure is largely phase separated with one interface where as at high inherent structure energy we find nucleation type growth. Interestingly, at low temperature, the average inherent structure energy (< EIS >) exhibits a drop with temperature which signals the onset of crystallization in one of the phases while the other remains in the liquid state. The nonideal composition dependence of viscosity is anticorrelated with average inherent structure energy.

High swirl-inducing piston bowls in small diesel engines for emission reduction

Detailed three-dimensional CFD simulations involving flow and combustion chemistry are used to study the effect of swirl induced by re-entrant piston bowl geometries on pollutant emissions from a single-cylinder diesel engine. The baseline engine configuration consists of a hemispherical piston bowl and an injector with finite sac volume. The first iteration involved using a torroidal, slightly re-entrant bowl geometry, and a sac-less injector. Pollutant emission measurements indicated a reduction in emissions with this modification. Simulations on both configurations were then conducted to understand the effect of the changes. The simulation results indicate that the selected piston bowl geometry could actually be reducing the in-cylinder swirl and turbulence and the emission reduction may be entirely due to the introduction of the sac-less injector. In-cylinder air motion was then studied in a number of combustion chamber geometries, and a geometry which produced the highest in-cylinder swirl and Turbulence Kinetic Energy (TKE) around the compression top dead centre (TDC) was identified. The optimal nature of this re-entrant piston bowl geometry is confirmed by detailed combustion simulations and emission predictions. (C) 2010 Elsevier Ltd. All rights reserved.