53,759 research outputs found
The Overlooked Potential of Generalized Linear Models in Astronomy - I: Binomial Regression
Revealing hidden patterns in astronomical data is often the path to
fundamental scientific breakthroughs; meanwhile the complexity of scientific
inquiry increases as more subtle relationships are sought. Contemporary data
analysis problems often elude the capabilities of classical statistical
techniques, suggesting the use of cutting edge statistical methods. In this
light, astronomers have overlooked a whole family of statistical techniques for
exploratory data analysis and robust regression, the so-called Generalized
Linear Models (GLMs). In this paper -- the first in a series aimed at
illustrating the power of these methods in astronomical applications -- we
elucidate the potential of a particular class of GLMs for handling
binary/binomial data, the so-called logit and probit regression techniques,
from both a maximum likelihood and a Bayesian perspective. As a case in point,
we present the use of these GLMs to explore the conditions of star formation
activity and metal enrichment in primordial minihaloes from cosmological
hydro-simulations including detailed chemistry, gas physics, and stellar
feedback. We predict that for a dark mini-halo with metallicity , an increase of in the gas
molecular fraction, increases the probability of star formation occurrence by a
factor of 75%. Finally, we highlight the use of receiver operating
characteristic curves as a diagnostic for binary classifiers, and ultimately we
use these to demonstrate the competitive predictive performance of GLMs against
the popular technique of artificial neural networks.Comment: 20 pages, 10 figures, 3 tables, accepted for publication in Astronomy
and Computin
Steering in computational science: mesoscale modelling and simulation
This paper outlines the benefits of computational steering for high
performance computing applications. Lattice-Boltzmann mesoscale fluid
simulations of binary and ternary amphiphilic fluids in two and three
dimensions are used to illustrate the substantial improvements which
computational steering offers in terms of resource efficiency and time to
discover new physics. We discuss details of our current steering
implementations and describe their future outlook with the advent of
computational grids.Comment: 40 pages, 11 figures. Accepted for publication in Contemporary
Physic
Multi-Architecture Monte-Carlo (MC) Simulation of Soft Coarse-Grained Polymeric Materials: SOft coarse grained Monte-carlo Acceleration (SOMA)
Multi-component polymer systems are important for the development of new
materials because of their ability to phase-separate or self-assemble into
nano-structures. The Single-Chain-in-Mean-Field (SCMF) algorithm in conjunction
with a soft, coarse-grained polymer model is an established technique to
investigate these soft-matter systems. Here we present an im- plementation of
this method: SOft coarse grained Monte-carlo Accelera- tion (SOMA). It is
suitable to simulate large system sizes with up to billions of particles, yet
versatile enough to study properties of different kinds of molecular
architectures and interactions. We achieve efficiency of the simulations
commissioning accelerators like GPUs on both workstations as well as
supercomputers. The implementa- tion remains flexible and maintainable because
of the implementation of the scientific programming language enhanced by
OpenACC pragmas for the accelerators. We present implementation details and
features of the program package, investigate the scalability of our
implementation SOMA, and discuss two applications, which cover system sizes
that are difficult to reach with other, common particle-based simulation
methods
Computational Physics on Graphics Processing Units
The use of graphics processing units for scientific computations is an
emerging strategy that can significantly speed up various different algorithms.
In this review, we discuss advances made in the field of computational physics,
focusing on classical molecular dynamics, and on quantum simulations for
electronic structure calculations using the density functional theory, wave
function techniques, and quantum field theory.Comment: Proceedings of the 11th International Conference, PARA 2012,
Helsinki, Finland, June 10-13, 201
Computational Investigations on Polymerase Actions in Gene Transcription and Replication Combining Physical Modeling and Atomistic Simulations
Polymerases are protein enzymes that move along nucleic acid chains and
catalyze template-based polymerization reactions during gene transcription and
replication. The polymerases also substantially improve transcription or
replication fidelity through the non-equilibrium enzymatic cycles. We briefly
review computational efforts that have been made toward understanding
mechano-chemical coupling and fidelity control mechanisms of the polymerase
elongation. The polymerases are regarded as molecular information motors during
the elongation process. It requires a full spectrum of computational approaches
from multiple time and length scales to understand the full polymerase
functional cycle. We keep away from quantum mechanics based approaches to the
polymerase catalysis due to abundant former surveys, while address only
statistical physics modeling approach and all-atom molecular dynamics
simulation approach. We organize this review around our own modeling and
simulation practices on a single-subunit T7 RNA polymerase, and summarize
commensurate studies on structurally similar DNA polymerases. For multi-subunit
RNA polymerases that have been intensively studied in recent years, we leave
detailed discussions on the simulation achievements to other computational
chemical surveys, while only introduce very recently published representative
studies, including our own preliminary work on structure-based modeling on
yeast RNA polymerase II. In the end, we quickly go through kinetic modeling on
elongation pauses and backtracking activities. We emphasize the fluctuation and
control mechanisms of the polymerase actions, highlight the non-equilibrium
physical nature of the system, and try to bring some perspectives toward
understanding replication and transcription regulation from single molecular
details to a genome-wide scale
- …