Optimization of parallel implementation of UNRES package for coarse-grained simulations to treat large proteins.

We report major algorithmic improvements of the UNRES package for physics-based coarse-grained simulations of proteins. These include (i) introduction of interaction lists to optimize computations, (ii) transforming the inertia matrix to a pentadiagonal form to reduce computing and memory requirements, (iii) removing explicit angles and dihedral angles from energy expressions and recoding the most time-consuming energy/force terms to minimize the number of operations and to improve numerical stability, (iv) using OpenMP to parallelize those sections of the code for which distributed-memory parallelization involves unfavorable computing/communication time ratio, and (v) careful memory management to minimize simultaneous access of distant memory sections. The new code enables us to run molecular dynamics simulations of protein systems with size exceeding 100,000 amino-acid residues, reaching over 1 ns/day (1 µs/day in all-atom timescale) with 24 cores for proteins of this size. Parallel performance of the code and comparison of its performance with that of AMBER, GROMACS and MARTINI 3 is presented.

Newton's method for nonlinear stochastic wave equations driven by one-dimensional Brownian motion.

We consider nonlinear stochastic wave equations driven by one-dimensional white noise with respect to time. The existence of solutions is proved by means of Picard iterations. Next we apply Newton's method. Moreover, a second-order convergence in a probabilistic sense is demonstrated.

Asymptotic Expansion Method with Respect to Small Parameter for Ternary Diffusion Models.

Ternary diffusion models lead to strongly coupled systems of PDEs. We choose the smallest diffusion coefficient as a small parameter in a power series expansion whose components fulfill relatively simple equations. Although this series is divergent, one can use its finite sums to derive feasible numerical approximations, e.g. finite difference methods (FDMs).

Structured populations with diffusion and Feller conditions.

We prove a weak maximum principle for structured population models with dynamic boundary conditions. We establish existence and positivity of solutions of these models and investigate the asymptotic behaviour of solutions. In particular, we analyse so called size profile.