Phase diagram of the bosonic double-exchange model

The phase diagram of the simplest approximation to double-exchange systems, the bosonic double-exchange model with antiferromagnetic (AFM) superexchange coupling, is fully worked out by means of Monte Carlo simulations, large-N expansions, and variational mean-field calculations. We find a rich phase diagram, with no first-order phase transitions. The most surprising finding is the existence of a segmentlike ordered phase at low temperature for intermediate AFM coupling which cannot be detected in neutron-scattering experiments. This is signaled by a maximum (a cusp) in the specific heat. Below the phase transition, only short-range ordering would be found in neutron scattering. Researchers looking for a quantum critical point in manganites should be wary of this possibility. Finite-size scaling estimates of critical exponents are presented, although large scaling corrections are present in the reachable lattice sizes.

Dynamical generation of a gauge symmetry in the double-exchange model

It is shown that a bosonic formulation of the double-exchange model, one of the classical models for magnetism, generates dynamically a gauge-invariant phase in a finite region of the phase diagram. We use analytical methods, Monte Carlo simulations and finite-size scaling analysis. We study the transition line between that region and the paramagnetic phase. The numerical results show that this transition line belongs to the universality class of the antiferromagnetic RP^(2) model. The fact that one can define a universality class for the antiferromagnetic RP^(2) model, different from the one of the O(N) models, is puzzling and somehow contradicts naive expectations about universality.

The reliability mathematical model of load-sharing parallel systems based on fatigue cumulative damage

Conventional reliability models for parallel systems are not applicable for the analysis of parallel systems with load transfer and sharing. In this short communication, firstly, the dependent failures of parallel systems are analyzed, and the reliability model of load-sharing parallel system is presented based on Miner cumulative damage theory and the full probability formula. Secondly, the parallel system reliability is calculated by Monte Carlo simulation when the component life follows the Weibull distribution. The research result shows that the proposed reliability mathematical model could analyze and evaluate the reliability of parallel systems in the presence of load transfer.

Bayesian Emulation for Sequential Modeling, Inference and Decision Analysis

The advances in three related areas of state-space modeling, sequential Bayesian learning, and decision analysis are addressed, with the statistical challenges of scalability and associated dynamic sparsity. The key theme that ties the three areas is Bayesian model emulation: solving challenging analysis/computational problems using creative model emulators. This idea defines theoretical and applied advances in non-linear, non-Gaussian state-space modeling, dynamic sparsity, decision analysis and statistical computation, across linked contexts of multivariate time series and dynamic networks studies. Examples and applications in financial time series and portfolio analysis, macroeconomics and internet studies from computational advertising demonstrate the utility of the core methodological innovations.

Chapter 1 summarizes the three areas/problems and the key idea of emulating in those areas. Chapter 2 discusses the sequential analysis of latent threshold models with use of emulating models that allows for analytical filtering to enhance the efficiency of posterior sampling. Chapter 3 examines the emulator model in decision analysis, or the synthetic model, that is equivalent to the loss function in the original minimization problem, and shows its performance in the context of sequential portfolio optimization. Chapter 4 describes the method for modeling the steaming data of counts observed on a large network that relies on emulating the whole, dependent network model by independent, conjugate sub-models customized to each set of flow. Chapter 5 reviews those advances and makes the concluding remarks.

A Molecular-scale Programmable Stochastic Process Based On Resonance Energy Transfer Networks: Modeling And Applications

While molecular and cellular processes are often modeled as stochastic processes, such as Brownian motion, chemical reaction networks and gene regulatory networks, there are few attempts to program a molecular-scale process to physically implement stochastic processes. DNA has been used as a substrate for programming molecular interactions, but its applications are restricted to deterministic functions and unfavorable properties such as slow processing, thermal annealing, aqueous solvents and difficult readout limit them to proof-of-concept purposes. To date, whether there exists a molecular process that can be programmed to implement stochastic processes for practical applications remains unknown.

In this dissertation, a fully specified Resonance Energy Transfer (RET) network between chromophores is accurately fabricated via DNA self-assembly, and the exciton dynamics in the RET network physically implement a stochastic process, specifically a continuous-time Markov chain (CTMC), which has a direct mapping to the physical geometry of the chromophore network. Excited by a light source, a RET network generates random samples in the temporal domain in the form of fluorescence photons which can be detected by a photon detector. The intrinsic sampling distribution of a RET network is derived as a phase-type distribution configured by its CTMC model. The conclusion is that the exciton dynamics in a RET network implement a general and important class of stochastic processes that can be directly and accurately programmed and used for practical applications of photonics and optoelectronics. Different approaches to using RET networks exist with vast potential applications. As an entropy source that can directly generate samples from virtually arbitrary distributions, RET networks can benefit applications that rely on generating random samples such as 1) fluorescent taggants and 2) stochastic computing.

By using RET networks between chromophores to implement fluorescent taggants with temporally coded signatures, the taggant design is not constrained by resolvable dyes and has a significantly larger coding capacity than spectrally or lifetime coded fluorescent taggants. Meanwhile, the taggant detection process becomes highly efficient, and the Maximum Likelihood Estimation (MLE) based taggant identification guarantees high accuracy even with only a few hundred detected photons.

Meanwhile, RET-based sampling units (RSU) can be constructed to accelerate probabilistic algorithms for wide applications in machine learning and data analytics. Because probabilistic algorithms often rely on iteratively sampling from parameterized distributions, they can be inefficient in practice on the deterministic hardware traditional computers use, especially for high-dimensional and complex problems. As an efficient universal sampling unit, the proposed RSU can be integrated into a processor / GPU as specialized functional units or organized as a discrete accelerator to bring substantial speedups and power savings.

Alkenone stable carbon isotopes, carbonate stable carbon and oxygen isotopes, planktic foraminifera boron isotopes and atmospheric CO2 calculations and box model results from Ras il-Pellegrin section, Malta

The development of a permanent, stable ice sheet in East Antarctica happened during the middle Miocene, about 14 million years (Myr) ago. The middle Miocene therefore represents one of the distinct phases of rapid change in the transition from the "greenhouse" of the early Eocene to the "icehouse" of the present day. Carbonate carbon isotope records of the period immediately following the main stage of ice sheet development reveal a major perturbation in the carbon system, represented by the positive d13C excursion known as carbon maximum 6 ("M6"), which has traditionally been interpreted as reflecting increased burial of organic matter and atmospheric pCO2 drawdown. More recently, it has been suggested that the d13C excursion records a negative feedback resulting from the reduction of silicate weathering and an increase in atmospheric pCO2. Here we present high-resolution multi-proxy (alkenone carbon and foraminiferal boron isotope) records of atmospheric carbon dioxide and sea surface temperature across CM6. Similar to previously published records spanning this interval, our records document a world of generally low (~300 ppm) atmospheric pCO2 at a time generally accepted to be much warmer than today. Crucially, they also reveal a pCO2 decrease with associated cooling, which demonstrates that the carbon burial hypothesis for CM6 is feasible and could have acted as a positive feedback on global cooling.

The Rotational Evolution and Magnetospheric Emission of the Magnetic Early B-type Stars

How do the magnetic fields of massive stars evolve over time? Are their gyrochronological ages consistent with ages inferred from evolutionary tracks? Why do most stars predicted to host Centrifugal Magnetospheres (CMs) display no H$\alpha$ emission? Does plasma escape from CMs via centrifugal breakout events, or by a steady-state leakage mechanism? This thesis investigates these questions via a population study with a sample of 51 magnetic early B-type stars. The longitudinal magnetic field \bz~was measured from Least Squares Deconvolution profiles extracted from high-resolution spectropolarimetric data. New rotational periods $P_{\rm rot}$ were determined for 15 stars from \bz, leaving only 3 stars for which $P_{\rm rot}$ is unknown. Projected rotational velocities \vsini~were measured from multiple spectral lines. Effective temperatures and surface gravities were measured via ionization balances and line profile fitting of H Balmer lines. Fundamental physical parameters, \bz, \vsini, and $P_{\rm rot}$ were then used to determine radii, masses, ages, dipole oblique rotator model, stellar wind, magnetospheric, and spindown parameters using a Monte Carlo approach that self-consistently calculates all parameters while accounting for all available constraints on stellar properties. Dipole magnetic field strengths $B_{\rm d}$ follow a log-normal distribution similar to that of Ap stars, and decline over time in a fashion consistent with the expected conservation of fossil magnetic flux. $P_{\rm rot}$ increases with fractional main sequence age, mass, and $B_{\rm d}$, as expected from magnetospheric braking. However, comparison of evolutionary track ages to maximum spindown ages $t_{\rm S,max}$ shows that initial rotation fractions may be far below critical for stars with $M_*>10 M_\odot$. Computing $t_{\rm S,max}$ with different mass-loss prescriptions indicates that the mass-loss rates of B-type stars are likely much lower than expected from extrapolation from O-type stars. Stars with H$\alpha$ in emission and absorption occupy distinct regions in the updated rotation-magnetic confinement diagram: H$\alpha$-bright stars are found to be younger, more rapidly rotating, and more strongly magnetized than the general population. Emission strength is sensitive both to the volume of the CM and to the mass-loss rate, favouring leakage over centrifugal breakout.

Human behaviour recognition in data-scarce domains

This paper presents the novel theory for performing multi-agent activity recognition without requiring large training corpora. The reduced need for data means that robust probabilistic recognition can be performed within domains where annotated datasets are traditionally unavailable. Complex human activities are composed from sequences of underlying primitive activities. We do not assume that the exact temporal ordering of primitives is necessary, so can represent complex activity using an unordered bag. Our three-tier architecture comprises low-level video tracking, event analysis and high-level inference. High-level inference is performed using a new, cascading extension of the Rao–Blackwellised Particle Filter. Simulated annealing is used to identify pairs of agents involved in multi-agent activity. We validate our framework using the benchmarked PETS 2006 video surveillance dataset and our own sequences, and achieve a mean recognition F-Score of 0.82. Our approach achieves a mean improvement of 17% over a Hidden Markov Model baseline.

Legionella pneumophila secretes a mitochondrial carrier protein during infection

The Mitochondrial Carrier Family (MCF) is a signature group of integral membrane proteins that transport metabolites across the mitochondrial inner membrane in eukaryotes. MCF proteins are characterized by six transmembrane segments that assemble to form a highly-selective channel for metabolite transport. We discovered a novel MCF member, termed Legionellanucleotide carrier Protein (LncP), encoded in the genome of Legionella pneumophila, the causative agent of Legionnaire's disease. LncP was secreted via the bacterial Dot/Icm type IV secretion system into macrophages and assembled in the mitochondrial inner membrane. In a yeast cellular system, LncP induced a dominant-negative phenotype that was rescued by deleting an endogenous ATP carrier. Substrate transport studies on purified LncP reconstituted in liposomes revealed that it catalyzes unidirectional transport and exchange of ATP transport across membranes, thereby supporting a role for LncP as an ATP transporter. A hidden Markov model revealed further MCF proteins in the intracellular pathogens, Legionella longbeachae and Neorickettsia sennetsu, thereby challenging the notion that MCF proteins exist exclusively in eukaryotic organisms.

IntelliChair: a non-intrusive sitting posture and sitting activity recognition system

Current Ambient Intelligence and Intelligent Environment research focuses on the interpretation of a subject’s behaviour at the activity level by logging the Activity of Daily Living (ADL) such as eating, cooking, etc. In general, the sensors employed (e.g. PIR sensors, contact sensors) provide low resolution information. Meanwhile, the expansion of ubiquitous computing allows researchers to gather additional information from different types of sensor which is possible to improve activity analysis. Based on the previous research about sitting posture detection, this research attempts to further analyses human sitting activity. The aim of this research is to use non-intrusive low cost pressure sensor embedded chair system to recognize a subject’s activity by using their detected postures. There are three steps for this research, the first step is to find a hardware solution for low cost sitting posture detection, second step is to find a suitable strategy of sitting posture detection and the last step is to correlate the time-ordered sitting posture sequences with sitting activity. The author initiated a prototype type of sensing system called IntelliChair for sitting posture detection. Two experiments are proceeded in order to determine the hardware architecture of IntelliChair system. The prototype looks at the sensor selection and integration of various sensor and indicates the best for a low cost, non-intrusive system. Subsequently, this research implements signal process theory to explore the frequency feature of sitting posture, for the purpose of determining a suitable sampling rate for IntelliChair system. For second and third step, ten subjects are recruited for the sitting posture data and sitting activity data collection. The former dataset is collected byasking subjects to perform certain pre-defined sitting postures on IntelliChair and it is used for posture recognition experiment. The latter dataset is collected by asking the subjects to perform their normal sitting activity routine on IntelliChair for four hours, and the dataset is used for activity modelling and recognition experiment. For the posture recognition experiment, two Support Vector Machine (SVM) based classifiers are trained (one for spine postures and the other one for leg postures), and their performance evaluated. Hidden Markov Model is utilized for sitting activity modelling and recognition in order to establish the selected sitting activities from sitting posture sequences.2. After experimenting with possible sensors, Force Sensing Resistor (FSR) is selected as the pressure sensing unit for IntelliChair. Eight FSRs are mounted on the seat and back of a chair to gather haptic (i.e., touch-based) posture information. Furthermore, the research explores the possibility of using alternative non-intrusive sensing technology (i.e. vision based Kinect Sensor from Microsoft) and find out the Kinect sensor is not reliable for sitting posture detection due to the joint drifting problem. A suitable sampling rate for IntelliChair is determined according to the experiment result which is 6 Hz. The posture classification performance shows that the SVM based classifier is robust to “familiar” subject data (accuracy is 99.8% with spine postures and 99.9% with leg postures). When dealing with “unfamiliar” subject data, the accuracy is 80.7% for spine posture classification and 42.3% for leg posture classification. The result of activity recognition achieves 41.27% accuracy among four selected activities (i.e. relax, play game, working with PC and watching video). The result of this thesis shows that different individual body characteristics and sitting habits influence both sitting posture and sitting activity recognition. In this case, it suggests that IntelliChair is suitable for individual usage but a training stage is required.

Analysis of radiation induced DNA damage by LC-MS/MS in isolated and cellular DNA

Abstract: It is well established that ionizing radiation induces a variety of damage in DNA by direct effects that are mediated by one-electron oxidation and indirect effects that are mediated by the reaction of water radiolysis products, e.g., hydroxyl radicals (•OH). In cellular DNA, direct and indirect effects appear to have about an equal effect toward DNA damage. We have shown that ϒ-(gamma) ray irradiation of aqueous solutions of DNA, during which •OH is the major damaging ROS can lead to the formation several lesions. On the other hand, the methylation and oxidative demethylation of cytosine in CpG dinucleotides plays a critical role in the gene regulation. The C5 position of cytosine in CG dinucleotides is frequently methylated by DNA methyl transferees (DNMTs) and constitutes 4-5% of the total cytosine. Here, my PhD research work focuses on the analysis of oxidative base modifications of model compounds of methylated and non methylated oligonucleotides, isolated DNA (calf-thymus DNA) and F98 cultured cell by gamma radiation. In addition, we identified a series of modifications of the 2-deoxyribose moiety of DNA arising from the exposure of isolated and cellular DNA to ionizing radiation. We also studied one electron oxidation of cellular DNA in cultured human HeLa cells initiated by intense nanosecond 266 nm laser pulse irradiation, which produces cross-links between guanine and thymine bases (G*-T*). To achieve these goals, we developed several methods based on mass spectrometry to analyze base modifications in isolated DNA and cellular DNA.

High-temperature superconductivity in the two-dimensional t-J model: Gutzwiller wavefunction solution

A systematic diagrammatic expansion for Gutzwiller wavefunctions (DE-GWFs) proposed very recently is used for the description of the superconducting (SC) ground state in the two-dimensional square-lattice t-J model with the hopping electron amplitudes t (and t') between nearest (and next-nearest) neighbors. For the example of the SC state analysis we provide a detailed comparison of the method's results with those of other approaches. Namely, (i) the truncated DE-GWF method reproduces the variational Monte Carlo (VMC) results and (ii) in the lowest (zeroth) order of the expansion the method can reproduce the analytical results of the standard Gutzwiller approximation (GA), as well as of the recently proposed 'grand-canonical Gutzwiller approximation' (called either GCGA or SGA). We obtain important features of the SC state. First, the SC gap at the Fermi surface resembles a d(x2-y2) wave only for optimally and overdoped systems, being diminished in the antinodal regions for the underdoped case in a qualitative agreement with experiment. Corrections to the gap structure are shown to arise from the longer range of the real-space pairing. Second, the nodal Fermi velocity is almost constant as a function of doping and agrees semi-quantitatively with experimental results. Third, we compare the

Probabilistic approach of water residence time and connectivity using Markov Chains with application to tidal embayments

Markov Chain analysis was recently proposed to assess the time scales and preferential pathways into biological or physical networks by computing residence time, first passage time, rates of transfer between nodes and number of passages in a node. We propose to adapt an algorithm already published for simple systems to physical systems described with a high resolution hydrodynamic model. The method is applied to bays and estuaries on the Eastern Coast of Canada for their interest in shellfish aquaculture. Current velocities have been computed by using a 2 dimensional grid of elements and circulation patterns were summarized by averaging Eulerian flows between adjacent elements. Flows and volumes allow computing probabilities of transition between elements and to assess the average time needed by virtual particles to move from one element to another, the rate of transfer between two elements, and the average residence time of each system. We also combined transfer rates and times to assess the main pathways of virtual particles released in farmed areas and the potential influence of farmed areas on other areas. We suggest that Markov chain is complementary to other sets of ecological indicators proposed to analyse the interactions between farmed areas - e.g. depletion index, carrying capacity assessment. Markov Chain has several advantages with respect to the estimation of connectivity between pair of sites. It makes possible to estimate transfer rates and times at once in a very quick and efficient way, without the need to perform long term simulations of particle or tracer concentration.