Weak coupling of acoustic-like phonons and magnon dynamics in incommensurate spin ladder compound Sr14Cu24O41

Intriguing lattice dynamics has been predicted for aperiodic crystals that contain incommensurate substructures. Here we report inelastic neutron scattering measurements of phonon and magnon dispersions in Sr14Cu24O41, which contains incommensurate one-dimensional (1D) chain and two-dimensional (2D) ladder substructures. Two distinct acoustic phonon-like modes, corresponding to the sliding motion of one sublattice against the other, are observed for atomic motions polarized along the incommensurate axis. In the long wavelength limit, it is found that the sliding mode shows a remarkably small energy gap of 1.7-1.9 meV, indicating very weak interactions between the two incommensurate sublattices. The measurements also reveal a gapped and steep linear magnon dispersion of the ladder sublattice. The high group velocity of this magnon branch and weak coupling with acoustic phonons can explain the large magnon thermal conductivity in Sr14Cu24O41 crystals. In addition, the magnon specific heat is determined from the measured total specific heat and phonon density of states, and exhibits a Schottky anomaly due to gapped magnon modes of the spin chains. These findings offer new insights into the phonon and magnon dynamics and thermal transport properties of incommensurate magnetic crystals that contain low-dimensional substructures.

On Provable Algorithms for Determination of Continuous Protein Interdomain Motions from Residual Dipolar Couplings

Dynamics of biomolecules over various spatial and time scales are essential for biological functions such as molecular recognition, catalysis and signaling. However, reconstruction of biomolecular dynamics from experimental observables requires the determination of a conformational probability distribution. Unfortunately, these distributions cannot be fully constrained by the limited information from experiments, making the problem an ill-posed one in the terminology of Hadamard. The ill-posed nature of the problem comes from the fact that it has no unique solution. Multiple or even an infinite number of solutions may exist. To avoid the ill-posed nature, the problem needs to be regularized by making assumptions, which inevitably introduce biases into the result.

Here, I present two continuous probability density function approaches to solve an important inverse problem called the RDC trigonometric moment problem. By focusing on interdomain orientations we reduced the problem to determination of a distribution on the 3D rotational space from residual dipolar couplings (RDCs). We derived an analytical equation that relates alignment tensors of adjacent domains, which serves as the foundation of the two methods. In the first approach, the ill-posed nature of the problem was avoided by introducing a continuous distribution model, which enjoys a smoothness assumption. To find the optimal solution for the distribution, we also designed an efficient branch-and-bound algorithm that exploits the mathematical structure of the analytical solutions. The algorithm is guaranteed to find the distribution that best satisfies the analytical relationship. We observed good performance of the method when tested under various levels of experimental noise and when applied to two protein systems. The second approach avoids the use of any model by employing maximum entropy principles. This 'model-free' approach delivers the least biased result which presents our state of knowledge. In this approach, the solution is an exponential function of Lagrange multipliers. To determine the multipliers, a convex objective function is constructed. Consequently, the maximum entropy solution can be found easily by gradient descent methods. Both algorithms can be applied to biomolecular RDC data in general, including data from RNA and DNA molecules.

Sequential anomaly detection in the presence of noise and limited feedback

This paper describes a methodology for detecting anomalies from sequentially observed and potentially noisy data. The proposed approach consists of two main elements: 1) filtering, or assigning a belief or likelihood to each successive measurement based upon our ability to predict it from previous noisy observations and 2) hedging, or flagging potential anomalies by comparing the current belief against a time-varying and data-adaptive threshold. The threshold is adjusted based on the available feedback from an end user. Our algorithms, which combine universal prediction with recent work on online convex programming, do not require computing posterior distributions given all current observations and involve simple primal-dual parameter updates. At the heart of the proposed approach lie exponential-family models which can be used in a wide variety of contexts and applications, and which yield methods that achieve sublinear per-round regret against both static and slowly varying product distributions with marginals drawn from the same exponential family. Moreover, the regret against static distributions coincides with the minimax value of the corresponding online strongly convex game. We also prove bounds on the number of mistakes made during the hedging step relative to the best offline choice of the threshold with access to all estimated beliefs and feedback signals. We validate the theory on synthetic data drawn from a time-varying distribution over binary vectors of high dimensionality, as well as on the Enron email dataset. © 1963-2012 IEEE.

Information relaxations and duality in stochastic dynamic programs

We describe a general technique for determining upper bounds on maximal values (or lower bounds on minimal costs) in stochastic dynamic programs. In this approach, we relax the nonanticipativity constraints that require decisions to depend only on the information available at the time a decision is made and impose a "penalty" that punishes violations of nonanticipativity. In applications, the hope is that this relaxed version of the problem will be simpler to solve than the original dynamic program. The upper bounds provided by this dual approach complement lower bounds on values that may be found by simulating with heuristic policies. We describe the theory underlying this dual approach and establish weak duality, strong duality, and complementary slackness results that are analogous to the duality results of linear programming. We also study properties of good penalties. Finally, we demonstrate the use of this dual approach in an adaptive inventory control problem with an unknown and changing demand distribution and in valuing options with stochastic volatilities and interest rates. These are complex problems of significant practical interest that are quite difficult to solve to optimality. In these examples, our dual approach requires relatively little additional computation and leads to tight bounds on the optimal values. © 2010 INFORMS.

The empirical analysis of cigarette tax avoidance and illicit trade in Vietnam, 1998-2010.

Illicit trade carries the potential to magnify existing tobacco-related health care costs through increased availability of untaxed and inexpensive cigarettes. What is known with respect to the magnitude of illicit trade for Vietnam is produced primarily by the industry, and methodologies are typically opaque. Independent assessment of the illicit cigarette trade in Vietnam is vital to tobacco control policy. This paper measures the magnitude of illicit cigarette trade for Vietnam between 1998 and 2010 using two methods, discrepancies between legitimate domestic cigarette sales and domestic tobacco consumption estimated from surveys, and trade discrepancies as recorded by Vietnam and trade partners. The results indicate that Vietnam likely experienced net smuggling in during the period studied. With the inclusion of adjustments for survey respondent under-reporting, inward illicit trade likely occurred in three of the four years for which surveys were available. Discrepancies in trade records indicate that the value of smuggled cigarettes into Vietnam ranges from $100 million to $300 million between 2000 and 2010 and that these cigarettes primarily originate in Singapore, Hong Kong, Macao, Malaysia, and Australia. Notable differences in trends over time exist between the two methods, but by comparison, the industry estimates consistently place the magnitude of illicit trade at the upper bounds of what this study shows. The unavailability of annual, survey-based estimates of consumption may obscure the true, annual trend over time. Second, as surveys changed over time, estimates relying on them may be inconsistent with one another. Finally, these two methods measure different components of illicit trade, specifically consumption of illicit cigarettes regardless of origin and smuggling of cigarettes into a particular market. However, absent a gold standard, comparisons of different approaches to illicit trade measurement serve efforts to refine and improve measurement approaches and estimates.

The Function and Regulation of Photobodies in Phytochrome Signaling

Light is a critical environmental signal that regulates every phase of the plant life cycle, from germination to floral initiation. Of the many light receptors in the model plant Arabidopsis thaliana, the red- and far-red light-sensing phytochromes (phys) are arguably the best studied, but the earliest events in the phy signaling pathway remain poorly understood. One of the earliest phy signaling events is the translocation of photoactivated phys from the cytoplasm to the nucleus, where they localize to subnuclear foci termed photobodies; in continuous light, photobody localization correlates closely with the light-dependent inhibition of embryonic stem growth. Despite a growing body of evidence supporting the biological significance of photobodies in light signaling, photobodies have also been shown to be dispensable for seedling growth inhibition in continuous light, so their physiological importance remains controversial; additionally, the molecular components that are required for phy localization to photobodies are largely unknown. The overall goal of my dissertation research was to gain insight into the early steps of phy signaling by further defining the role of photobodies in this process and identifying additional intragenic and extragenic requirements for phy localization to photobodies.

Even though the domain structure of phys has been extensively studied, not all of the intramolecular requirements for phy localization to photobodies are known. Previous studies have shown that the entire C-terminus of phys is both necessary and sufficient for their localization to photobodies. However, the importance of the individual subdomains of the C-terminus is still unclear. For example a truncation lacking part of the most C-terminal domain, the histidine kinase-related domain (HKRD), can still localize to small photobodies in the light and behaves like a weak allele. However, a point mutation within the HKRD renders the entire molecule completely inactive. To resolve this discrepancy, I explored the hypothesis that this point mutation might impair the dimerization of the HKRD; dimerization has been shown to occur via the C-terminus of phy and is required for more efficient signaling. I show that this point mutation impairs nuclear localization of phy as well as its subnuclear localization to photobodies. Additionally, yeast-two-hybrid analysis shows that the wild-type HKRD can homodimerize but that the HKRD containing the point mutation fails to dimerize with both itself and with wild-type HKRD. These results demonstrate that dimerization of the HKRD is required for both nuclear and photobody localization of phy.

Studies of seedlings grown in diurnal conditions show that photoactivated phy can persist into darkness to repress seedling growth; a seedling's growth rate is therefore fastest at the end of the night. To test the idea that photobodies could be involved in regulating seedling growth in the dark, I compared the growth of two transgenic Arabidopsis lines, one in which phy can localize to photobodies (PBG), and one in which it cannot (NGB). Despite these differences in photobody morphology, both lines are capable of transducing light signals and inhibiting seedling growth in continuous light. After the transition from red light to darkness, the PBG line was able to repress seedling growth, as well as the accumulation of the growth-promoting, light-labile transcription factor PHYTOCHROME INTERACTING FACTOR 3 (PIF3), for eighteen hours, and this correlated perfectly with the presence of photobodies. Reducing the amount of active phy by either reducing the light intensity or adding a phy-inactivating far-red pulse prior to darkness led to faster accumulation of PIF3 and earlier seedling growth. In contrast, the NGB line accumulated PIF3 even in the light, and seedling growth was only repressed for six hours; this behavior was similar in NGB regardless of the light treatment. These results suggest that photobodies are required for the degradation of PIF3 and for the prolonged stabilization of active phy in darkness. They also support the hypothesis that photobody localization of phys could serve as an instructive cue during the light-to-dark transition, thereby fine-tuning light-dependent responses in darkness.

In addition to determining an intragenic requirement for photobody localization and further exploring the significance of photobodies in phy signaling, I wanted to identify extragenic regulators of photobody localization. A recent study identified one such factor, HEMERA (HMR); hmr mutants do not form large photobodies, and they are tall and albino in the light. To identify other components in the HMR-mediated branch of the phy signaling pathway, I performed a forward genetic screen for suppressors of a weak hmr allele. Surprisingly, the first three mutants isolated from the screen were alleles of the same novel gene, SON OF HEMERA (SOH). The soh mutations rescue all of the phenotypes associated with the weak hmr allele, and they do so in an allele-specific manner, suggesting a direct interaction between SOH and HMR. Null soh alleles, which were isolated in an independent, tall, albino screen, are defective in photobody localization, demonstrating that SOH is an extragenic regulator of phy localization to photobodies that works in the same genetic pathway as HMR.

In this work, I show that dimerization of the HKRD is required for both the nuclear and photobody localization of phy. I also demonstrate a tight correlation between photobody localization and PIF3 degradation, further establishing the significance of photobodies in phy signaling. Finally, I identify a novel gene, SON OF HEMERA, whose product is necessary for phy localization to photobodies in the light, thereby isolating a new extragenic determinant of photobody localization. These results are among the first to focus exclusively on one of the earliest cellular responses to light - photobody localization of phys - and they promise to open up new avenues into the study of a poorly understood facet of the phy signaling pathway.

Evolution of networks and sequences in eukaryotic cell cycle control.

The molecular networks regulating the G1-S transition in budding yeast and mammals are strikingly similar in network structure. However, many of the individual proteins performing similar network roles appear to have unrelated amino acid sequences, suggesting either extremely rapid sequence evolution, or true polyphyly of proteins carrying out identical network roles. A yeast/mammal comparison suggests that network topology, and its associated dynamic properties, rather than regulatory proteins themselves may be the most important elements conserved through evolution. However, recent deep phylogenetic studies show that fungal and animal lineages are relatively closely related in the opisthokont branch of eukaryotes. The presence in plants of cell cycle regulators such as Rb, E2F and cyclins A and D, that appear lost in yeast, suggests cell cycle control in the last common ancestor of the eukaryotes was implemented with this set of regulatory proteins. Forward genetics in non-opisthokonts, such as plants or their green algal relatives, will provide direct information on cell cycle control in these organisms, and may elucidate the potentially more complex cell cycle control network of the last common eukaryotic ancestor.

A theory of hypoellipticity and unique ergodicity for semilinear stochastic PDEs

We present a theory of hypoellipticity and unique ergodicity for semilinear parabolic stochastic PDEs with "polynomial" nonlinearities and additive noise, considered as abstract evolution equations in some Hilbert space. It is shown that if Hörmander's bracket condition holds at every point of this Hilbert space, then a lower bound on the Malliavin covariance operatorμt can be obtained. Informally, this bound can be read as "Fix any finite-dimensional projection on a subspace of sufficiently regular functions. Then the eigenfunctions of μt with small eigenvalues have only a very small component in the image of Π." We also show how to use a priori bounds on the solutions to the equation to obtain good control on the dependency of the bounds on the Malliavin matrix on the initial condition. These bounds are sufficient in many cases to obtain the asymptotic strong Feller property introduced in [HM06]. One of the main novel technical tools is an almost sure bound from below on the size of "Wiener polynomials," where the coefficients are possibly non-adapted stochastic processes satisfying a Lips chitz condition. By exploiting the polynomial structure of the equations, this result can be used to replace Norris' lemma, which is unavailable in the present context. We conclude by showing that the two-dimensional stochastic Navier-Stokes equations and a large class of reaction-diffusion equations fit the framework of our theory.

Phylotastic! Making tree-of-life knowledge accessible, reusable and convenient.

BACKGROUND: Scientists rarely reuse expert knowledge of phylogeny, in spite of years of effort to assemble a great "Tree of Life" (ToL). A notable exception involves the use of Phylomatic, which provides tools to generate custom phylogenies from a large, pre-computed, expert phylogeny of plant taxa. This suggests great potential for a more generalized system that, starting with a query consisting of a list of any known species, would rectify non-standard names, identify expert phylogenies containing the implicated taxa, prune away unneeded parts, and supply branch lengths and annotations, resulting in a custom phylogeny suited to the user's needs. Such a system could become a sustainable community resource if implemented as a distributed system of loosely coupled parts that interact through clearly defined interfaces. RESULTS: With the aim of building such a "phylotastic" system, the NESCent Hackathons, Interoperability, Phylogenies (HIP) working group recruited 2 dozen scientist-programmers to a weeklong programming hackathon in June 2012. During the hackathon (and a three-month follow-up period), 5 teams produced designs, implementations, documentation, presentations, and tests including: (1) a generalized scheme for integrating components; (2) proof-of-concept pruners and controllers; (3) a meta-API for taxonomic name resolution services; (4) a system for storing, finding, and retrieving phylogenies using semantic web technologies for data exchange, storage, and querying; (5) an innovative new service, DateLife.org, which synthesizes pre-computed, time-calibrated phylogenies to assign ages to nodes; and (6) demonstration projects. These outcomes are accessible via a public code repository (GitHub.com), a website (http://www.phylotastic.org), and a server image. CONCLUSIONS: Approximately 9 person-months of effort (centered on a software development hackathon) resulted in the design and implementation of proof-of-concept software for 4 core phylotastic components, 3 controllers, and 3 end-user demonstration tools. While these products have substantial limitations, they suggest considerable potential for a distributed system that makes phylogenetic knowledge readily accessible in computable form. Widespread use of phylotastic systems will create an electronic marketplace for sharing phylogenetic knowledge that will spur innovation in other areas of the ToL enterprise, such as annotation of sources and methods and third-party methods of quality assessment.