Adiabatic eigenfunction-based approach for coherent excitation transfer: An almost analytical treatment of the Fenna-Matthews-Olson complex

We suggest a method of studying coherence in finite-level systems coupled to the environment and use it for the Hamiltonian that has been used to describe the light-harvesting pigment-protein complex. The method works with the adiabatic states and transforms the Hamiltonian to a form in which the terms responsible for decoherence and population relaxation are separated out. Decoherence is then accounted for nonperturbatively and population relaxation using a Markovian master equation. Almost analytical results can be obtained for the seven-level system, and the calculations are very simple for systems with more levels. We apply the treatment to the seven-level system, and the results are in excellent agreement with the exact numerical results of Nalbach et al. Nalbach, Braun, and Thorwart, Phys. Rev. E 84, 041926 (2011)]. Our approach is able to account for decoherence and population relaxation separately. It is found that decoherence causes only damping of oscillations and does not lead to transfer to the reaction center. Population relaxation is necessary for efficient transfer to the reaction center, in agreement with earlier findings. Our results show that the transformation to the adiabatic basis followed by a Redfield type of approach leads to results in good agreement with exact simulation.

Adiabatic Eigenfunction Based Approach to Coherent Transfer: Application to the Fenna-Matthews-Olson (FMO) Complex and the Role of Correlations in the Efficiency of Energy Transfer

We have recently suggested a method (Pallavi Bhattacharyya and K. L. Sebastian, Physical Review E 2013, 87, 062712) for the analysis of coherence in finite-level systems that are coupled to the surroundings and used it to study the process of energy transfer in the Fenna-Matthews-Olson (FMO) complex. The method makes use of adiabatic eigenstates of the Hamiltonian, with a subsequent transformation of the Hamiltonian into a form where the terms responsible for decoherence and population relaxation could be separated out at the lowest order. Thus one can account for decoherence nonperturbatively, and a Markovian type of master equation could be used for evaluating the population relaxation. In this paper, we apply this method to a two-level system as well as to a seven-level system. Comparisons with exact numerical results show that the method works quite well and is in good agreement with numerical calculations. The technique can be applied with ease to systems with larger numbers of levels as well. We also investigate how the presence of correlations among the bath degrees of freedom of the different bacteriochlorophyll a molecules of the FMO Complex affect the rate of energy transfer. Surprisingly, in the cases that we studied, our calculations suggest that the presence of anticorrelations, in contrast to correlations, make the excitation transfer more facile.

Matrix-Assisted Refolding, Purification and Activity Assessment Using a `Form Invariant' Assay for Matrix Metalloproteinase 2 (MMP2)

Matrix metalloproteinases expression is used as biomarker for various cancers and associated malignancies. Since these proteinases can cleave many intracellular proteins, overexpression tends to be toxic; hence, a challenge to purify them. To overcome these limitations, we designed a protocol where full length pro-MMP2 enzyme was overexpressed in E. coli as inclusion bodies and purified using 6xHis affinity chromatography under denaturing conditions. In one step, the enzyme was purified and refolded directly on the affinity matrix under redox conditions to obtain a bioactive protein. The pro-MMP2 protein was characterized by mass spectrometry, CD spectroscopy, zymography and activity analysis using a simple in-house developed `form invariant' assay, which reports the total MMP2 activity independent of its various forms. The methodology yielded higher yields of bioactive protein compared to other strategies reported till date, and we anticipate that using the protocol, other toxic proteins can also be overexpressed and purified from E. coli and subsequently refolded into active form using a one step renaturation protocol.

Lower bounds on Ricci flow invariant curvatures and geometric applications

We consider Ricci flow invariant cones C in the space of curvature operators lying between the cones ``nonnegative Ricci curvature'' and ``nonnegative curvature operator''. Assuming some mild control on the scalar curvature of the Ricci flow, we show that if a solution to the Ricci flow has its curvature operator which satisfies R + epsilon I is an element of C at the initial time, then it satisfies R + epsilon I is an element of C on some time interval depending only on the scalar curvature control. This allows us to link Gromov-Hausdorff convergence and Ricci flow convergence when the limit is smooth and R + I is an element of C along the sequence of initial conditions. Another application is a stability result for manifolds whose curvature operator is almost in C. Finally, we study the case where C is contained in the cone of operators whose sectional curvature is nonnegative. This allows us to weaken the assumptions of the previously mentioned applications. In particular, we construct a Ricci flow for a class of (not too) singular Alexandrov spaces.

NON-COERCIVE RICCI FLOW INVARIANT CURVATURE CONES

This note is a study of nonnegativity conditions on curvature preserved by the Ricci flow. We focus on a specific class of curvature conditions which we call non-coercive: These are the conditions for which nonnegative curvature and vanishing scalar curvature does not imply flatness. We show, in dimensions greater than 4, that if a Ricci flow invariant nonnegativity condition is satisfied by all Einstein curvature operators with nonnegative scalar curvature, then this condition is just the nonnegativity of scalar curvature. As a corollary, we obtain that a Ricci flow invariant curvature condition, which is stronger than a nonnegative scalar curvature, cannot be strictly satisfied by curvature operators (other than multiples of the identity) of compact Einstein symmetric spaces. We also investigate conditions which are satisfied by all conformally flat manifolds with nonnegative scalar curvature.

Exact Phase Retrieval in Principal Shift-Invariant Spaces

We address the problem of phase retrieval from Fourier transform magnitude spectrum for continuous-time signals that lie in a shift-invariant space spanned by integer shifts of a generator kernel. The phase retrieval problem for such signals is formulated as one of reconstructing the combining coefficients in the shift-invariant basis expansion. We develop sufficient conditions on the coefficients and the bases to guarantee exact phase retrieval, by which we mean reconstruction up to a global phase factor. We present a new class of discrete-domain signals that are not necessarily minimum-phase, but allow for exact phase retrieval from their Fourier magnitude spectra. We also establish Hilbert transform relations between log-magnitude and phase spectra for this class of discrete signals. It turns out that the corresponding continuous-domain counterparts need not satisfy a Hilbert transform relation; notwithstanding, the continuous-domain signals can be reconstructed from their Fourier magnitude spectra. We validate the reconstruction guarantees through simulations for some important classes of signals such as bandlimited signals and piecewise-smooth signals. We also present an application of the proposed phase retrieval technique for artifact-free signal reconstruction in frequency-domain optical-coherence tomography (FDOCT).

Statistics of an adiabatic charge pump

We investigate the effect of time-dependent cyclic-adiabatic driving on the charge transport in a quantum junction. We propose a nonequilibrium Green's function formalism to study the statistics of the charge pumped (at zero bias) through the junction. The formulation is used to demonstrate charge pumping in a single electronic level coupled to two (electronic) reservoirs with time-dependent couplings. An analytical expression for the average pumped current for a general cyclic driving is derived. It is found that for zero bias, for a certain class of driving, the Berry phase contributes only to the odd cumulants. In contrast, a quantum master equation formulation does not show a Berry-phase effect at all.

Evolution of multivariant microstructures with anisotropic misfit: A phase field study

We study, in two dimensions, the effect of misfit anisotropy on microstructural evolution during precipitation of an ordered beta phase from a disordered alpha matrix; these phases have, respectively, 2- and 6-fold rotation symmetries. Thus, precipitation produces three orientational variants of beta phase particles, and they have an anisotropic (and crystallographically equivalent) misfit strain with the matrix. The anisotropy in misfit is characterized using a parameter t = epsilon(yy)/epsilon(xx), where epsilon(xx) and epsilon(yy) are the principal components of the misfit strain tensor. Our phase field, simulations show that the morphology of beta phase particles is significantly influenced by 1, the level of misfit anisotropy. Particles are circular in systems with dilatational misfit (t = 1), elongated along the direction of lower principal misfit when 0 < t < 1 and elongated along the invariant direction when - 1 <= t <= 0. In the special case of a pure shear misfit strain (t = - 1), the microstructure exhibits star, wedge and checkerboard patterns; these microstructural features are in agreement with those in Ti-Al-Nb alloys.

Self-organized criticality in dynamics without branching

We demonstrate the phenomenon of self-organized criticality (SOC) in a simple random walk model described by a random walk of a myopic ant, i.e., a walker who can see only nearest neighbors. The ant acts on the underlying lattice aiming at uniform digging, i.e., reduction of the height profile of the surface but is unaffected by the underlying lattice. In one, two, and three dimensions we have explored this model and have obtained power laws in the time intervals between consecutive events of "digging." Being a simple random walk, the power laws in space translate to power laws in time. We also study the finite size scaling of asymptotic scale invariant process as well as dynamic scaling in this system. This model differs qualitatively from the cascade models of SOC.

Mechanical processing and microstructural control in hot working of hot isostatically pressed P/M IN-100 superalloy

The hot deformation behavior of hot isostatically pressed (HIPd) P/M IN-100 superalloy has been studied in the temperature range 1000-1200 degrees C and strain rate range 0.0003-10 s(-1) using hot compression testing. A processing map has been developed on the basis of these data and using the principles of dynamic materials modelling. The map exhibited three domains: one at 1050 degrees C and 0.01 s(-1), with a peak efficiency of power dissipation of approximate to 32%, the second at 1150 degrees C and 10 s(-1), with a peak efficiency of approximate to 36% and the third at 1200 degrees C and 0.1 s(-1), with a similar efficiency. On the basis of optical and electron microscopic observations, the first domain was interpreted to represent dynamic recovery of the gamma phase, the second domain represents dynamic recrystallization (DRX) of gamma in the presence of softer gamma', while the third domain represents DRX of the gamma phase only. The gamma' phase is stable upto 1150 degrees C, gets deformed below this temperature and the chunky gamma' accumulates dislocations, which at larger strains cause cracking of this phase. At temperatures lower than 1080 degrees C and strain rates higher than 0.1 s(-1), the material exhibits flow instability, manifested in the form of adiabatic shear bands. The material may be subjected to mechanical processing without cracking or instabilities at 1200 degrees C and 0.1 s(-1), which are the conditions for DRX of the gamma phase.

Modeling of sonochemical decomposition of CCI4 in aqueous solutions

In the context of removal of organic pollutants from wastewater, sonolysis of CCl4 dissolved in water has been widely investigated. These investigations are either completely experimental or correlate data empirically. In this work, a quantitative model is developed to predict the rate of sonolysis of aqueous CCl4. The model considers the isothermal growth and partially adiabatic collapse of cavitation bubbles containing gas and vapor leading to conditions of high temperatures and pressures in them, attainment of thermodynamic equilibrium at the end of collapse, release of bubble contents into the liquid pool, and reactions in the well-mixed pool. The model successfully predicts the extent of degradation of dissolved CCl4, and the influence of various parameters such as initial concentration of CCl4, temperature, and nature of gas atmosphere above the liquid. in particular, it predicts the results of Hua and Hoffmann (Environ. Sci Technol, 1996, 30, 864-871), who found that degradation is first order with CCl4 and that Argon as well as Ar-O-3 atmospheres give the same results. The framework of the model is capable of quantitatively describing the degradation of many dissolved organics by considering all the involved species.

Gaussian Schell-model beams and general shape invariance

We present two six-parameter families of anisotropic Gaussian Schell-model beams that propagate in a shape-invariant manner, with the intensity distribution continuously twisting about the beam axis. The two families differ in the sense or helicity of this beam twist. The propagation characteristics of these shape-invariant beams are studied, and the restrictions on the beam parameters that arise from the optical uncertainty principle are brought out. Shape invariance is traced to a fundamental dynamical symmetry that underlies these beams. This symmetry is the product of spatial rotation and fractional Fourier transformation.

PocketMatch: A new algorithm to compare binding sites in protein structures

Recognizing similarities and deriving relationships among protein molecules is a fundamental requirement in present-day biology. Similarities can be present at various levels which can be detected through comparison of protein sequences or their structural folds. In some cases similarities obscure at these levels could be present merely in the substructures at their binding sites. Inferring functional similarities between protein molecules by comparing their binding sites is still largely exploratory and not as yet a routine protocol. One of the main reasons for this is the limitation in the choice of appropriate analytical tools that can compare binding sites with high sensitivity. To benefit from the enormous amount of structural data that is being rapidly accumulated, it is essential to have high throughput tools that enable large scale binding site comparison. Results: Here we present a new algorithm PocketMatch for comparison of binding sites in a frame invariant manner. Each binding site is represented by 90 lists of sorted distances capturing shape and chemical nature of the site. The sorted arrays are then aligned using an incremental alignment method and scored to obtain PMScores for pairs of sites. A comprehensive sensitivity analysis and an extensive validation of the algorithm have been carried out. A comparison with other site matching algorithms is also presented. Perturbation studies where the geometry of a given site was retained but the residue types were changed randomly, indicated that chance similarities were virtually non-existent. Our analysis also demonstrates that shape information alone is insufficient to discriminate between diverse binding sites, unless combined with chemical nature of amino acids. Conclusion: A new algorithm has been developed to compare binding sites in accurate, efficient and high-throughput manner. Though the representation used is conceptually simplistic, we demonstrate that along with the new alignment strategy used, it is sufficient to enable binding comparison with high sensitivity. Novel methodology has also been presented for validating the algorithm for accuracy and sensitivity with respect to geometry and chemical nature of the site. The method is also fast and takes about 1/250(th) second for one comparison on a single processor. A parallel version on BlueGene has also been implemented.

Localisation of Plasmodium falciparum uroporphyrinogen III decarboxylase of the heme-biosynthetic pathway in the apicoplast and characterisation of its catalytic properties

Uroporphyrinogen decarboxylase (UROD) is a key enzyme in the heme-biosynthetic pathway and in Plasmodium falciparum it occupies a strategic position in the proposed hybrid pathway for heme biosynthesis involving shuttling of intermediates between different subcellular compartments in the parasite. In the present study, we demonstrate that an N-terminally truncated recombinant P. falciparum UROD (r(Δ)PfUROD) over-expressed and purified from Escherichia coli cells, as well as the native enzyme from the parasite were catalytically less efficient compared with the host enzyme, although they were similar in other enzyme parameters. Molecular modeling of PfUROD based on the known crystal structure of the human enzyme indicated that the protein manifests a distorted triose phosphate isomerase (TIM) barrel fold which is conserved in all the known structures of UROD. The parasite enzyme shares all the conserved or invariant amino acid residues at the active and substrate binding sites, but is rich in lysine residues compared with the host enzyme. Mutation of specific lysine residues corresponding to residues at the dimer interface in human UROD enhanced the catalytic efficiency of the enzyme and dimer stability indicating that the lysine rich nature and weak dimer interface of the wild-type PfUROD could be responsible for its low catalytic efficiency. PfUROD was localised to the apicoplast, indicating the requirement of additional mechanisms for transport of the product coproporphyrinogen to other subcellular sites for its further conversion and ultimate heme formation.

Synergy between tribo-oxidation and strain rate response on governing the dry sliding wear behavior of titanium.

The present work provides an insight into the dry sliding wear behavior of titanium based on synergy between tribo-oxidation and strain rate response. Pin-on-disc tribometer was used to characterize the friction and wear behavior of titanium pin in sliding contact with polycrystalline alumina disk under ambient and vacuum condition. The sliding speed was varied from 0.01 to 1.4 ms(-1), normal load was varied from 15.3 to 76 N and with a sliding distance of 1500 m. It was seen that dry sliding wear behavior of titanium was governed by combination of tribo-oxidation and strain rate response in near surface region of titanium. Strain rate response of titanium was recorded by conducting uni-axial compression tests at constant true strain rate of 100 s(-1) in the temperature range from 298 to 873 K. Coefficient of friction and wear rate were reduced with increased sliding speed from 0.01 to 1.0 ms(-1). This is attributed to the formation of in situ self lubricating oxide film (TiO) and reduction in the intensity of adiabatic shear band cracking in the near surface region. This trend was confirmed by performing series of dry sliding tests under vacuum condition of 2 x 10(-4) Torr. Characterization tools such as optical microscopy, scanning electron microscopy, and X-ray diffractometer provided evidence of such processes. These experimental findings can be applied to enhance the dry sliding wear behavior of titanium with proper choice of operating conditions such as sliding speed, normal load, and environment.