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.

Effect of quantum nuclear motion on hydrogen bonding

This work considers how the properties of hydrogen bonded complexes, X-H center dot center dot center dot Y, are modified by the quantum motion of the shared proton. Using a simple two-diabatic state model Hamiltonian, the analysis of the symmetric case, where the donor (X) and acceptor (Y) have the same proton affinity, is carried out. For quantitative comparisons, a parametrization specific to the O-H center dot center dot center dot O complexes is used. The vibrational energy levels of the one-dimensional ground state adiabatic potential of the model are used to make quantitative comparisons with a vast body of condensed phase data, spanning a donor-acceptor separation (R) range of about 2.4-3.0 angstrom, i.e., from strong to weak hydrogen bonds. The position of the proton (which determines the X-H bond length) and its longitudinal vibrational frequency, along with the isotope effects in both are described quantitatively. An analysis of the secondary geometric isotope effect, using a simple extension of the two-state model, yields an improved agreement of the predicted variation with R of frequency isotope effects. The role of bending modes is also considered: their quantum effects compete with those of the stretching mode for weak to moderate H-bond strengths. In spite of the economy in the parametrization of the model used, it offers key insights into the defining features of H-bonds, and semi-quantitatively captures several trends. (C) 2014 AIP Publishing LLC.

Quantum simulation using fidelity-profile optimization

Experimental quantum simulation of a Hamiltonian H requires unitary operator decomposition (UOD) of its evolution unitary U = exp(-iHt) in terms of native unitary operators of the experimental system. Here, using a genetic algorithm, we numerically evaluate the most generic UOD (valid over a continuous range of Hamiltonian parameters) of the unitary operator U, termed fidelity-profile optimization. The optimization is obtained by systematically evaluating the functional dependence of experimental unitary operators (such as single-qubit rotations and time-evolution unitaries of the system interactions) to the Hamiltonian (H) parameters. Using this technique, we have solved the experimental unitary decomposition of a controlled-phase gate (for any phase value), the evolution unitary of the Heisenberg XY interaction, and simulation of the Dzyaloshinskii-Moriya (DM) interaction in the presence of the Heisenberg XY interaction. Using these decompositions, we studied the entanglement dynamics of a Bell state in the DM interaction and experimentally verified the entanglement preservation procedure of Hou et al. Ann. Phys. (N.Y.) 327, 292 (2012)] in a nuclear magnetic resonance quantum information processor.

Majorana edge modes in the Kitaev model

We study the Majorana modes, both equilibrium and Floquet, which can appear at the edges of the Kitaev model on the honeycomb lattice. We first present the analytical solutions known for the equilibrium Majorana edge modes for both zigzag and armchair edges of a semi-infinite Kitaev model and chart the parameter regimes in which they appear. We then examine how edge modes can be generated if the Kitaev coupling on the bonds perpendicular to the edge is varied periodically in time as periodic delta-function kicks. We derive a general condition for the appearance and disappearance of the Floquet edge modes as a function of the drive frequency for a generic d-dimensional integrable system. We confirm this general condition for the Kitaev model with a finite width by mapping it to a one-dimensional model. Our numerical and analytical study of this problem shows that Floquet Majorana modes can appear on some edges in the kicked system even when the corresponding equilibrium Hamiltonian has no Majorana mode solutions on those edges. We support our analytical studies by numerics for a finite sized system which show that periodic kicks can generate modes at the edges and the corners of the lattice.

Topological blocking in quantum quench dynamics

We study the nonequilibrium dynamics of quenching through a quantum critical point in topological systems, focusing on one of their defining features: ground-state degeneracies and associated topological sectors. We present the notion of ``topological blocking,'' experienced by the dynamics due to a mismatch in degeneracies between two phases, and we argue that the dynamic evolution of the quench depends strongly on the topological sector being probed. We demonstrate this interplay between quench and topology in models stemming from two extensively studied systems, the transverse Ising chain and the Kitaev honeycomb model. Through nonlocal maps of each of these systems, we effectively study spinless fermionic p-wave paired topological superconductors. Confining the systems to ring and toroidal geometries, respectively, enables us to cleanly address degeneracies, subtle issues of fermion occupation and parity, and mismatches between topological sectors. We show that various features of the quench, which are related to Kibble-Zurek physics, are sensitive to the topological sector being probed, in particular, the overlap between the time-evolved initial ground state and an appropriate low-energy state of the final Hamiltonian. While most of our study is confined to translationally invariant systems, where momentum is a convenient quantum number, we briefly consider the effect of disorder and illustrate how this can influence the quench in a qualitatively different way depending on the topological sector considered.

Studies on a frustrated Heisenberg spin chain with alternating ferromagnetic and antiferromagnetic exchanges

We study Heisenberg spin-1/2 and spin-1 chains with alternating ferromagnetic (J(1)(F)) and antiferromagnetic (J(1)(A)) nearest-neighbor interactions and a ferromagnetic next-nearest-neighbor interaction (J(2)(F)). In this model frustration is present due to the non-zero J(2)(F). The model with site spin s behaves like a Haldane spin chain, with site spin 2s in the limit of vanishing J(2)(F) and large J(1)(F)/J(1)(A). We show that the exact ground state of the model can be found along a line in the parameter space. For fixed J(1)(F), the phase diagram in the space of J(1)(A)-J(2)(F) is determined using numerical techniques complemented by analytical calculations. A number of quantities, including the structure factor, energy gap, entanglement entropy and zero temperature magnetization, are studied to understand the complete phase diagram. An interesting and potentially important feature of this model is that it can exhibit a macroscopic magnetization jump in the presence of a magnetic field; we study this using an effective Hamiltonian.

A FUNCTIONAL MODEL FOR PURE Gamma-CONTRACTIONS

A pair of commuting operators (S,P) defined on a Hilbert space H for which the closed symmetrized bidisc Gamma = {(z(1) + z(2), z(1)z(2)) : vertical bar z(1)vertical bar <= 1, vertical bar z(2)vertical bar <= 1} subset of C-2 is a spectral set is called a Gamma-contraction in the literature. A Gamma-contraction (S, P) is said to be pure if P is a pure contraction, i.e., P*(n) -> 0 strongly as n -> infinity Here we construct a functional model and produce a set of unitary invariants for a pure Gamma-contraction. The key ingredient in these constructions is an operator, which is the unique solution of the operator equation S - S*P = DpXDp, where X is an element of B(D-p), and is called the fundamental operator of the Gamma-contraction (S, P). We also discuss some important properties of the fundamental operator.

Edge states of a three-dimensional topological insulator

We use the bulk Hamiltonian for a three-dimensional topological insulator such as Bi-2 Se-3 to study the states which appear on its various surfaces and along the edge between two surfaces. We use both analytical methods based on the surface Hamiltonians (which are derived from the bulk Hamiltonian) and numerical methods based on a lattice discretization of the bulk Hamiltonian. We find that the application of a potential barrier along an edge can give rise to states localized at that edge. These states have an unusual energy-momentum dispersion which can be controlled by applying a potential along the edge; in particular, the velocity of these states can be tuned to zero. The scattering and conductance across the edge is studied as a function of the edge potential. We show that a magnetic field in a particular direction can also give rise to zero energy states on certain edges. We point out possible experimental ways of looking for the various edge states.

Improved quasiparticle wave functions and mean field for G(0)W(0) calculations: Initialization with the COHSEX operator

The GW approximation to the electron self-energy has become a standard method for ab initio calculation of excited-state properties of condensed-matter systems. In many calculations, the G W self-energy operator, E, is taken to be diagonal in the density functional theory (DFT) Kohn-Sham basis within the G0 W0 scheme. However, there are known situations in which this diagonal Go Wo approximation starting from DFT is inadequate. We present two schemes to resolve such problems. The first, which we called sc-COHSEX-PG W, involves construction of an improved mean field using the static limit of GW, known as COHSEX (Coulomb hole and screened exchange), which is significantly simpler to treat than GW W. In this scheme, frequency-dependent self energy E(N), is constructed and taken to be diagonal in the COHSEX orbitals after the system is solved self-consistently within this formalism. The second method is called off diagonal-COHSEX G W (od-COHSEX-PG W). In this method, one does not self-consistently change the mean-field starting point but diagonalizes the COHSEX Hamiltonian within the Kohn-Sham basis to obtain quasiparticle wave functions and uses the resulting orbitals to construct the G W E in the diagonal form. We apply both methods to a molecular system, silane, and to two bulk systems, Si and Ge under pressure. For silane, both methods give good quasiparticle wave functions and energies. Both methods give good band gaps for bulk silicon and maintain good agreement with experiment. Further, the sc-COHSEX-PGW method solves the qualitatively incorrect DFT mean-field starting point (having a band overlap) in bulk Ge under pressure.

Vertex Cover Gets Faster and Harder on Low Degree Graphs

The problem of finding an optimal vertex cover in a graph is a classic NP-complete problem, and is a special case of the hitting set question. On the other hand, the hitting set problem, when asked in the context of induced geometric objects, often turns out to be exactly the vertex cover problem on restricted classes of graphs. In this work we explore a particular instance of such a phenomenon. We consider the problem of hitting all axis-parallel slabs induced by a point set P, and show that it is equivalent to the problem of finding a vertex cover on a graph whose edge set is the union of two Hamiltonian Paths. We show the latter problem to be NP-complete, and also give an algorithm to find a vertex cover of size at most k, on graphs of maximum degree four, whose running time is 1.2637(k) n(O(1)).

Relative entropy in higher spin holography

We examine relative entropy in the context of the higher spin/CFT duality. We consider 3D bulk configurations in higher spin gravity which are dual to the vacuum and a high temperature state of a CFT with W-algebra symmetries in the presence of a chemical potential for a higher spin current. The relative entropy between these states is then evaluated using the Wilson line functional for holographic entanglement entropy. In the limit of small entangling intervals, the relative entropy should vanish for a generic quantum system. We confirm this behavior by showing that the difference in the expectation values of the modular Hamiltonian between the states matches with the difference in the entanglement entropy in the short-distance regime. Additionally, we compute the relative entropy of states corresponding to smooth solutions in the SL(2, Z) family with respect to the vacuum.

The Tetrablock as a Spectral Set

The tetrablock, roughly speaking, is the set of all linear fractional maps that map the open unit disc to itself. A formal definition of this inhomogeneous domain is given below. This paper considers triples of commuting bounded operators (A,B,P) that have the tetrablock as a spectral set. Such a triple is named a tetrablock contraction. The motivation comes from the success of model theory in another inhomogeneous domain, namely, the symmetrized bidisc F. A pair of commuting bounded operators (S,P) with Gamma as a spectral set is called a Gamma-contraction, and always has a dilation. The two domains are related intricately as the Lemma 3.2 below shows. Given a triple (A, B, P) as above, we associate with it a pair (F-1, F-2), called its fundamental operators. We show that (A,B,P) dilates if the fundamental operators F-1 and F-2 satisfy certain commutativity conditions. Moreover, the dilation space is no bigger than the minimal isometric dilation space of the contraction P. Whether these commutativity conditions are necessary, too, is not known. what we have shown is that if there is a tetrablock isometric dilation on the minimal isometric dilation space of P. then those commutativity conditions necessarily get imposed on the fundamental operators. En route, we decipher the structure of a tetrablock unitary (this is the candidate as the dilation triple) and a tertrablock isometry (the restriction of a tetrablock unitary to a joint invariant sub-space). We derive new results about r-contractions and apply them to tetrablock contractions. The methods applied are motivated by 11]. Although the calculations are lengthy and more complicated, they beautifully reveal that the dilation depends on the mutual relationship of the two fundamental operators, so that certain conditions need to be satisfied. The question of whether all tetrablock contractions dilate or not is unresolved.

Synthesis, crystal structure and spectroscopic and electrochemical properties of bridged trisbenzoato copper-zinc heterobinuclear complex of 2,2 `-bipyridine

The synthesis of the heterobinuclear copper-zinc complex CuZn(bz)(3)(bpy)(2)]ClO4 (bz = benzoate) from benzoic acid and bipyridine is described. Single crystal X-ray diffraction studies of the heterobinuclear complex reveals the geometry of the benzoato bridged Cu(II)-Zn(II) centre. The copper or zinc atom is pentacoordinate, with two oxygen atoms from bridging benzoato groups and two nitrogen atoms from one bipyridine forming an approximate plane and a bridging oxygen atom from a monodentate benzoate group. The Cu-Zn distance is 3.345 angstrom. The complex is normal paramagnetic having mu(eff) value equal to 1.75 BM, ruling out the possibility of Cu-Cu interaction in the structural unit. The ESR spectrum of the complex in CH3CN at RT exhibit an isotropic four line spectrum centred at g = 2.142 and hyperfine coupling constants A(av) = 63 x 10(-4) cm(-1), characteristic of a mononuclear square-pyramidal copper(II) complexes. At LNT, the complex shows an isotropic spectrum with g(parallel to) = 2.254 and g(perpendicular to) =2.071 and A(parallel to) = 160 x 10(-4) cm(-1). The Hamiltonian parameters are characteristic of distorted square pyramidal geometry. Cyclic voltammetric studies of the complex have indicated quasi-reversible behaviour in acetonitrile solution. (C) 2014 Elsevier B.V. All rights reserved.

A Family of Domains Associated with mu-Synthesis

We introduce a family of domains-which we call the -quotients-associated with an aspect of -synthesis. We show that the natural association that the symmetrized polydisc has with the corresponding spectral unit ball is also exhibited by the -quotient and its associated unit `` -ball''. Here, is the structured singular value for the case Specifically: we show that, for such an E, the Nevanlinna-Pick interpolation problem with matricial data in a unit `` -ball'', and in general position in a precise sense, is equivalent to a Nevanlinna-Pick interpolation problem for the associated -quotient. Along the way, we present some characterizations for the -quotients.

Heterochromatic paths in edge colored graphs without small cycles and heterochromatic-triangle-free graphs

Conditions for the existence of heterochromatic Hamiltonian paths and cycles in edge colored graphs are well investigated in literature. A related problem in this domain is to obtain good lower bounds for the length of a maximum heterochromatic path in an edge colored graph G. This problem is also well explored by now and the lower bounds are often specified as functions of the minimum color degree of G - the minimum number of distinct colors occurring at edges incident to any vertex of G - denoted by v(G). Initially, it was conjectured that the lower bound for the length of a maximum heterochromatic path for an edge colored graph G would be 2v(G)/3]. Chen and Li (2005) showed that the length of a maximum heterochromatic path in an edge colored graph G is at least v(G) - 1, if 1 <= v(G) <= 7, and at least 3v(G)/5] + 1 if v(G) >= 8. They conjectured that the tight lower bound would be v(G) - 1 and demonstrated some examples which achieve this bound. An unpublished manuscript from the same authors (Chen, Li) reported to show that if v(G) >= 8, then G contains a heterochromatic path of length at least 120 + 1. In this paper, we give lower bounds for the length of a maximum heterochromatic path in edge colored graphs without small cycles. We show that if G has no four cycles, then it contains a heterochromatic path of length at least v(G) - o(v(G)) and if the girth of G is at least 4 log(2)(v(G)) + 2, then it contains a heterochromatic path of length at least v(G) - 2, which is only one less than the bound conjectured by Chen and Li (2005). Other special cases considered include lower bounds for the length of a maximum heterochromatic path in edge colored bipartite graphs and triangle-free graphs: for triangle-free graphs we obtain a lower bound of 5v(G)/6] and for bipartite graphs we obtain a lower bound of 6v(G)-3/7]. In this paper, it is also shown that if the coloring is such that G has no heterochromatic triangles, then G contains a heterochromatic path of length at least 13v(G)/17)]. This improves the previously known 3v(G)/4] bound obtained by Chen and Li (2011). We also give a relatively shorter and simpler proof showing that any edge colored graph G contains a heterochromatic path of length at least (C) 2015 Elsevier Ltd. All rights reserved.