The complexity of molecular dynamics simulations necessitates dimension reduction and coarsegraining techniques to enable tractable computation. The generalized Langevin equation (GLE) describes coarse-grained dynamics in reduced dimensions. In spite of playing a crucial role in nonequilibrium dynamics, the memory kernel of the GLE is often ignored because it is difficult to characterize and expensive to solve. To address these issues, we construct a data-driven rational approximation to the GLE. Building upon previous work leveraging the GLE to simulate simple systems, we extend these results to more complex molecules, whose many degrees of freedom and complicated dynamics require approximation methods. We demonstrate the effectiveness of our approximation by testing it against exact methods and comparing observables such as autocorrelation and transition rates.

In this paper, we investigate the long-term behaviour of solutions of the periodic Sigmoid Beverton-Holt equationwhere the a n and δ n are p-periodic positive sequences. Under certain conditions, there are shown to exist an asymptotically stable p-periodic state and a p-periodic Allee state with the property that populations smaller than the Allee state are driven to extinction while populations greater than the Allee state approach the stable state, thus accounting for the long-term behaviour of all initial states. This appears to be the first study of the equation with variable δ. The results are discussed with possible interpretations in Population Dynamics with emphasis on fish populations and smooth cordgrass.

ABSTRACTSelective laser sintering methods are workhorses for additively manufacturing polymer-based components. The ease of rapid prototyping also means it is easy to produce illicit components. It is necessary to have a data-calibrated in-situ physical model of the build process in order to predict expected and defective microstructure characteristics that inform component provenance. Toward this end, sintering models are calibrated and characteristics such as component defects are explored. This is accomplished by assimilating multiple data streams, imaging analysis, and computational model predictions in an adaptive Bayesian parameter estimation algorithm. From these data sources, along with a phase-field model, bulk porosity distributions are inferred. Model parameters are constrained to physically-relevant search directions by sensitivity analysis, and then matched to predictions using adaptive sampling. Using this feedback loop, data-constrained estimates of sintering model parameters along with uncertainty bounds are obtained.

Molecular dynamics (MD) simulations are used in biochemistry, physics, and other fields to study the motions, thermodynamic properties, and the interactions between molecules. Computational limitations and the complexity of these problems, however, create the need for approximations to the standard MD methods and for uncertainty quantification and reliability assessment of those approximations. In this paper, we exploit the intrinsic two-scale nature of MD to construct a class of large-scale dynamics approximations. The reliability of these methods are evaluated here by measuring the differences between full, classical MD simulations and those based on these large-scale approximations. Molecular dynamics evolutions are non-linear and chaotic, so the complete details of molecular evolutions cannot be accurately predicted even using full, classical MD simulations. This paper provides numerical results that demonstrate the existence of computationally efficient large-scale MD approximations which accurately model certain large-scale properties of the molecules: the energy, the linear and angular momenta, and other macroscopic features of molecular motions.

No abstract