Convex Optimization Algorithms and Statistical Bounds for Learning Structured Models

Thesis (Ph.D.)--University of Washington, 2016-08

Truncation identities for the small polaron fusion hierarchy

We study a one-dimensional lattice model of interacting spinless fermions. This model is integrable for both periodic and open boundary conditions; the latter case includes the presence of Grassmann valued non-diagonal boundary fields breaking the bulk U(1) symmetry of the model. Starting from the embedding of this model into a graded Yang-Baxter algebra, an infinite hierarchy of commuting transfer matrices is constructed by means of a fusion procedure. For certain values of the coupling constant related to anisotropies of the underlying vertex model taken at roots of unity, this hierarchy is shown to truncate giving a finite set of functional equations for the spectrum of the transfer matrices. For generic coupling constants, the spectral problem is formulated in terms of a functional (or TQ-)equation which can be solved by Bethe ansatz methods for periodic and diagonal open boundary conditions. Possible approaches for the solution of the model with generic non-diagonal boundary fields are discussed.

Nonlinear breathing modes at a defect

Recent molecular dynamics (MD) simulations of Cubero et al (1999) of a DNA duplex containing the 'rogue' base difluorotoluene (F) in place of a thymine (T) base show that breathing events can occur on the nanosecond timescale, whereas breathing events in a normal DNA duplex take place on the microsecond timescale. The main aim of this paper is to analyse a nonlinear Klein-Gordon lattice model of the DNA duplex including both nonlinear interactions between opposing bases and a defect in the interaction at one lattice site; each of which can cause localisation of energy. Solutions for a breather mode either side of the defect are derived using multiple-scales asymptotics and are pieced together across the defect to form a solution which includes the effects of the nonlinearity and the defect. We consider defects in the inter-chain interactions and in the along chain interactions. In most cases we find in-phase breather modes and/or out-of-phase breather modes, with one case displaying a shifted mode.

Robust infraestructure design in rapid transit rail systems

Incidents and rolling stock breakdowns are commonplace in rapid transit rail systems and may disrupt the system performance imposing deviations from planned operations. A network design model is proposed for reducing the effect of disruptions less likely to occur. Failure probabilities are considered functions of the amount of services and the rolling stock’s routing on the designed network so that they cannot be calculated a priori but result from the design process itself. A two recourse stochastic programming model is formulated where the failure probabilities are an implicit function of the number of services and routing of the transit lines.

Manipulating the structural and electronic properties of epitaxial NaNbO3 films via strain and stoichiometry

Due to their intriguing dielectric, pyroelectric, elasto-electric, or opto-electric properties, oxide ferroelectrics are vital candidates for the fabrication of most electronics. However, these extraordinary properties exist mainly in the temperature regime around the ferroelectric phase transition, which is usually several hundreds of K away from room temperature. Therefore, the manipulation of oxide ferroelectrics, especially moving the ferroelectric transition towards room temperature, is of great interest for application and also basic research. In this thesis, we demonstrate this using examples of NaNbO3 films. We show that the transition temperature of these films can be modified via plastic strain caused by epitaxial film growth on a structurally mismatched substrate, and this strain can be fixed by controlling the stoichiometry. The structural and electronic properties of Na1+xNbO3+δ thin films are carefully examined by among others XRD (e.g. RSM) and TEM and cryoelectronic measurements. Especially the electronic features are carefully analyzed via specially developed interdigitated electrodes in combination with integrated temperature sensor and heater. The electronic data are interpreted using existing as well as novel theories and models, they are proved to be closely correlated to the structural characteristics. The major results are: -Na1+xNbO3+δ thin films can be grown epitaxially on (110)NdGaO3 with a thickness up to 140 nm (thicker films have not been studied). Plastic relaxation of the compressive strain sets in when the thickness of the film exceeds approximately 10 – 15 nm. Films with excess Na are mainly composed of NaNbO3 with minor contribution of Na3NbO4. The latter phase seems to form nanoprecipitates that are homogeneously distributed in the NaNbO3 film which helps to stabilize the film and reduce the relaxation of the strain. -For the nominally stoichiometric films, the compressive strain leads to a broad and frequency-dispersive phase transition at lower temperature (125 – 147 K). This could be either a new transition or a shift in temperature of a known transition. Considering the broadness and frequency dispersion of the transition, this is actually a transition from the dielectric state at high temperature to a relaxor-type ferroelectric state at low temperature. The latter is based on the formation of polar nano-regions (PNRs). Using the electric field dependence of the freezing temperature, allows a direct estimation of the volume (70 to 270 nm3) and diameter (5.2 to 8 nm, spherical approximation) of the PNRs. The values confirm with literature values which were measured by other technologies. -In case of the off-stoichiometric samples, we observe again the classical ferroelectric behavior. However, the thermally hysteretic phase transition which is observed around 620 – 660 K for unstrained material is shifted to room temperature due to the compressive strain. Beside to the temperature shift, the temperature dependence of the permittivity is nearly identical for strained and unstrained materials. -The last but not least, in all cases, a significant anisotropy in the electronic and structural properties is observed which arises automatically from the anisotropic strain caused by the orthorhombic structure of the substrate. However, this anisotropy cannot be explained by the classical model which tries to fit an orthorhombic film onto an orthorhombic substrate. A novel “square lattice” model in which the films adapt a “square” shaped lattice in the plane of the film during the epitaxial growth at elevated temperature (~1000 K) nicely explains the experimental results. In this thesis we sketch a way to manipulate the ferroelectricity of NaNbO3 films via strain and stoichiometry. The results indicate that compressive strain which is generated by the epitaxial growth of the film on mismatched substrate is able to reduce the ferroelectric transition temperature or induce a phase transition at low temperature. Moreover, by adding Na in the NaNbO3 film a secondary phase Na3NbO4 is formed which seems to stabilize the main phase NaNbO3 and the strain and, thus, is able to engineer the ferroelectric behavior from the expected classical ferroelectric for perfect stoichiometry to relaxor-type ferroelectric for slightly off-stoichiometry, back to classical ferroelectric for larger off-stoichiometry. Both strain and stoichiometry are proven as perfect methods to optimize the ferroelectric properties of oxide films.

Modelos estocásticos de crescimento individual e desenvolvimento de software de estimação e previsão

Os modelos de crescimento individual são geralmente adaptações de modelos de crescimento de populações. Inicialmente estes modelos eram apenas determinísticos, isto é, não incorporavam as flutuações aleatórias do ambiente. Com o desenvolvimento da teoria do cálculo estocástico podemos adicionar um termo estocástico, que representa a aleatoriedade ambiental que influencia o processo em estudo. Actualmente, o estudo do crescimento individual em ambiente aleatório é cada vez mais importante, não apenas pela vertente financeira, mas também devido às suas aplicações nas áreas da saúde e da pecuária, entre outras. Problemas como o ajustamento de modelos de crescimento individual, estimação de parâmetros e previsão de tamanhos futuros são tratados neste trabalho. São apresentadas novas aplicações do modelo estocástico monomolecular generalizado e um novo software de aplicação deste e de outros modelos. ABSTRACT: Individual growth models are usually adaptations of growth population models. Initially these models were only deterministic, that is, they did not incorporate the random fluctuations of the environment. With the development of the theory of stochastic calculus, we can add a stochastic term that represents the random environmental influences in the process under study. Currently, the study of individual growth in a random environment is increasingly important, not only by the financial scope but also because of its applications in health care and livestock production, among others. Problems such as adjustment of an individual growth model, estimation of parameters and prediction of future sizes are treated in this work. New applications of the generalized stochastic monomolecular model and a new software applied to this and other models are presented.

Vehicle routing and location routing problems with minimum latency: algorithms and models

Latency can be defined as the sum of the arrival times at the customers. Minimum latency problems are specially relevant in applications related to humanitarian logistics. This thesis presents algorithms for solving a family of vehicle routing problems with minimum latency. First the latency location routing problem (LLRP) is considered. It consists of determining the subset of depots to be opened, and the routes that a set of homogeneous capacitated vehicles must perform in order to visit a set of customers such that the sum of the demands of the customers assigned to each vehicle does not exceed the capacity of the vehicle. For solving this problem three metaheuristic algorithms combining simulated annealing and variable neighborhood descent, and an iterated local search (ILS) algorithm, are proposed. Furthermore, the multi-depot cumulative capacitated vehicle routing problem (MDCCVRP) and the multi-depot k-traveling repairman problem (MDk-TRP) are solved with the proposed ILS algorithm. The MDCCVRP is a special case of the LLRP in which all the depots can be opened, and the MDk-TRP is a special case of the MDCCVRP in which the capacity constraints are relaxed. Finally, a LLRP with stochastic travel times is studied. A two-stage stochastic programming model and a variable neighborhood search algorithm are proposed for solving the problem. Furthermore a sampling method is developed for tackling instances with an infinite number of scenarios. Extensive computational experiments show that the proposed methods are effective for solving the problems under study.

Numerical study on viscoelastic behavior of polypropylene fiber reinforced concrete subjected to temperature variations using LDPM theory

This thesis is focused on the viscoelastic behavior of macro-synthetic fiber-reinforced concrete (MSFRC) with polypropylene studied numerically when subjected to temperature variations (-30 oC to +60 oC). LDPM (lattice discrete particle model), a meso-scale model for heterogeneous composites, is used. To reproduce the MSFRC structural behavior, an extended version of LDPM that includes fiber effects through fiber-concrete interface micromechanics, called LDPM-F, is applied. Model calibration is performed based on three-point bending, cube, and cylinder test for plain concrete and MSFRC. This is followed by a comprehensive literature study on the variation of mechanical properties with temperature for individual fibers and plain concrete. This literature study and past experimental test results constitute inputs for final numerical simulations. The numerical response of MSFRC three-point bending test is replicated and compared with the previously conducted experimental test results; finally, the conclusions were drawn. LDPM numerical model is successfully calibrated using experimental responses on plain concrete. Fiber-concrete interface micro-mechanical parameters are subsequently fixed and LDPM-F models are calibrated based on MSFRC three-point bending test at room temperature. Number of fibers contributing crack bridging mechanism is computed and found to be in good agreement with experimental counts. Temperature variations model for individual constituents of MSFRC, fibers and plain concrete, are implemented in LDPM-F. The model is validated for MSFRC three-point bending stress-CMOD (crack mouth opening) response reproduced at -30 oC, -15 oC, 0 oC, +20 oC, +40 oC and +60 oC. It is found that the model can well describe the temperature variation behavior of MSFRC. At positive temperatures, simulated responses are in good agreement. Slight disagreement in negative regimes suggests an in-depth study on fiber-matrix interface bond behavior with varying temperatures.

Critical behavior of the susceptible-infected-recovered model on a square lattice

By means of numerical simulations and epidemic analysis, the transition point of the stochastic asynchronous susceptible-infected-recovered model on a square lattice is found to be c(0)=0.176 500 5(10), where c is the probability a chosen infected site spontaneously recovers rather than tries to infect one neighbor. This point corresponds to an infection/recovery rate of lambda(c)=(1-c(0))/c(0)=4.665 71(3) and a net transmissibility of (1-c(0))/(1+3c(0))=0.538 410(2), which falls between the rigorous bounds of the site and bond thresholds. The critical behavior of the model is consistent with the two-dimensional percolation universality class, but local growth probabilities differ from those of dynamic percolation cluster growth, as is demonstrated explicitly.

A new scale-invariant ratio and finite-size scaling for the stochastic susceptible-infected-recovered model

The critical behavior of the stochastic susceptible-infected-recovered model on a square lattice is obtained by numerical simulations and finite-size scaling. The order parameter as well as the distribution in the number of recovered individuals is determined as a function of the infection rate for several values of the system size. The analysis around criticality is obtained by exploring the close relationship between the present model and standard percolation theory. The quantity UP, equal to the ratio U between the second moment and the squared first moment of the size distribution multiplied by the order parameter P, is shown to have, for a square system, a universal value 1.0167(1) that is the same for site and bond percolation, confirming further that the SIR model is also in the percolation class.

Markovian versus non-Markovian stochastic quantization of a complex-action model

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

A spatial stochastic model for rumor transmission

We consider an interacting particle system representing the spread of a rumor by agents on the d-dimensional integer lattice. Each agent may be in any of the three states belonging to the set {0,1,2}. Here 0 stands for ignorants, 1 for spreaders and 2 for stiflers. A spreader tells the rumor to any of its (nearest) ignorant neighbors at rate lambda. At rate alpha a spreader becomes a stifler due to the action of other (nearest neighbor) spreaders. Finally, spreaders and stiflers forget the rumor at rate one. We study sufficient conditions under which the rumor either becomes extinct or survives with positive probability.

Phase diagram of a model for a binary mixture of nematic molecules on a Bethe lattice

We investigate the phase diagram of a discrete version of the Maier-Saupe model with the inclusion of additional degrees of freedom to mimic a distribution of rodlike and disklike molecules. Solutions of this problem on a Bethe lattice come from the analysis of the fixed points of a set of nonlinear recursion relations. Besides the fixed points associated with isotropic and uniaxial nematic structures, there is also a fixed point associated with a biaxial nematic structure. Due to the existence of large overlaps of the stability regions, we resorted to a scheme to calculate the free energy of these structures deep in the interior of a large Cayley tree. Both thermodynamic and dynamic-stability analyses rule out the presence of a biaxial phase, in qualitative agreement with previous mean-field results.

Dynamic transitions in a three dimensional associating lattice gas model

The aggregation of interacting Brownian particles in sheared concentrated suspensions is an important issue in colloid and soft matter science per se. Also, it serves as a model to understand biochemical reactions occurring in vivo where both crowding and shear play an important role. We present an effective medium approach within the Smoluchowski equation with shear which allows one to calculate the encounter kinetics through a potential barrier under shear at arbitrary colloid concentrations. Experiments on a model colloidal system in simple shear flow support the validity of the model in the concentration range considered. By generalizing Kramers' rate theory to the presence of shear and collective hydrodynamics, our model explains the significant increase in the shear-induced reaction-limited aggregation kinetics upon increasing the colloid concentration.

Dynamic transitions in a two dimensional associating lattice gas model

Using Monte Carlo simulations we investigate some new aspects of the phase diagram and the behavior of the diffusion coefficient in an associating lattice gas (ALG) model on different regions of the phase diagram. The ALG model combines a two dimensional lattice gas where particles interact through a soft core potential and orientational degrees of freedom. The competition between soft core potential and directional attractive forces results in a high density liquid phase, a low density liquid phase, and a gas phase. Besides anomalies in the behavior of the density with the temperature at constant pressure and of the diffusion coefficient with density at constant temperature are also found. The two liquid phases are separated by a coexistence line that ends in a bicritical point. The low density liquid phase is separated from the gas phase by a coexistence line that ends in tricritical point. The bicritical and tricritical points are linked by a critical lambda-line. The high density liquid phase and the fluid phases are separated by a second critical tau-line. We then investigate how the diffusion coefficient behaves on different regions of the chemical potential-temperature phase diagram. We find that diffusivity undergoes two types of dynamic transitions: a fragile-to-strong transition when the critical lambda-line is crossed by decreasing the temperature at a constant chemical potential; and a strong-to-strong transition when the critical tau-line is crossed by decreasing the temperature at a constant chemical potential.