Numerical algebraic geometry for model selection and its application to the life sciences

Researchers working with mathematical models are often confronted by the related problems of parameter estimation, model validation, and model selection. These are all optimization problems, well-known to be challenging due to non-linearity, non-convexity and multiple local optima. Furthermore, the challenges are compounded when only partial data is available. Here, we consider polynomial models (e.g., mass-action chemical reaction networks at steady state) and describe a framework for their analysis based on optimization using numerical algebraic geometry. Specifically, we use probability-one polynomial homotopy continuation methods to compute all critical points of the objective function, then filter to recover the global optima. Our approach exploits the geometric structures relating models and data, and we demonstrate its utility on examples from cell signaling, synthetic biology, and epidemiology.Comment: References added, additional clarification

J-factors of short DNA molecules

The propensity of short DNA sequences to convert to the circular form is studied by a mesoscopic Hamiltonian method which incorporates both the bending of the molecule axis and the intrinsic twist of the DNA strands. The base pair fluctuations with respect to the helix diameter are treated as path trajectories in the imaginary time path integral formalism. The partition function for the sub-ensemble of closed molecules is computed by imposing chain ends boundary conditions both on the radial fluctuations and on the angular degrees of freedom. The cyclization probability, the J-factor, proves to be highly sensitive to the stacking potential, mostly to its nonlinear parameters. We find that the J-factor generally decreases by reducing the sequence length ( N ) and, more significantly, below N = 100 base pairs. However, even for very small molecules, the J-factors remain sizeable in line with recent experimental indications. Large bending angles between adjacent base pairs and anharmonic stacking appear as the causes of the helix flexibility at short length scales.Comment: The Journal of Chemical Physics - May 2016 ; 9 page

Probing a non-biaxial behavior of infinitely thin hard platelets

We give a criterion to test a non-biaxial behavior of infinitely thin hard platelets of $D_{2h}$ symmetry based upon the components of three order parameter tensors. We investigated the nematic behavior of monodisperse infinitely thin rectangular hard platelet systems by using the criterion. Starting with a square platelet system, and we compared it with rectangular platelet systems of various aspect ratios. For each system, we performed equilibration runs by using isobaric Monte Carlo simulations. Each system did not show a biaxial nematic behavior but a uniaxial nematic one, despite of the shape anisotropy of those platelets. The relationship between effective diameters by simulations and theoretical effective diameters of the above systems was also determined.Comment: Submitted to JPS

Radiation Induced Damage in GaAs Particle Detectors

The motivation for investigating the use of GaAs as a material for detecting particles in experiments for High Energy Physics (HEP) arose from its perceived resistance to radiation damage. This is a vital requirement for detector materials that are to be used in experiments at future accelerators where the radiation environments would exclude all but the most radiation resistant of detector types.Comment: 5 pages. PS file only - original in WORD Also available at http://ppewww.ph.gla.ac.uk/preprints/97/06

Vibronic interactions in the visible and near-infrared spectra of C60− anions

Electron-phonon coupling is an important factor in understanding many properties of the C60 fullerides. However, there has been little success in quantifying the strength of the vibronic coupling in C60 ions, with considerable disagreement between experimental and theoretical results. We will show that neglect of quadratic coupling in previous models for C60- ions results in a significant overestimate of the linear coupling constants. Including quadratic coupling allows a coherent interpretation to be made of earlier experimental and theoretical results which at first sight are incompatible

The Liouville-type theorem for integrable Hamiltonian systems with incomplete flows

For integrable Hamiltonian systems with two degrees of freedom whose Hamiltonian vector fields have incomplete flows, an analogue of the Liouville theorem is established. A canonical Liouville fibration is defined by means of an "exact" 2-parameter family of flat polygons equipped with certain pairing of sides. For the integrable Hamiltonian systems given by the vector field $v=(-\partial f/\partial w, \partial f/\partial z)$ on ${\mathbb C}^2$ where $f=f(z,w)$ is a complex polynomial in 2 variables, geometric properties of Liouville fibrations are described.Comment: 6 page

Geometry and topology of knotted ring polymers in an array of obstacles

We study knotted polymers in equilibrium with an array of obstacles which models confinement in a gel or immersion in a melt. We find a crossover in both the geometrical and the topological behavior of the polymer. When the polymers' radius of gyration, $R_G$, and that of the region containing the knot, $R_{G,k}$, are small compared to the distance b between the obstacles, the knot is weakly localised and $R_G$ scales as in a good solvent with an amplitude that depends on knot type. In an intermediate regime where $R_G > b > R_{G,k}$, the geometry of the polymer becomes branched. When $R_{G,k}$ exceeds b, the knot delocalises and becomes also branched. In this regime, $R_G$ is independent of knot type. We discuss the implications of this behavior for gel electrophoresis experiments on knotted DNA in weak fields.Comment: 4 pages, 6 figure