1
|
Chen A, Ren Q, Zhou T, Burrage P, Tian T, Burrage K. Balanced implicit Patankar-Euler methods for positive solutions of stochastic differential equations of biological regulatory systems. J Chem Phys 2024; 160:064117. [PMID: 38353308 DOI: 10.1063/5.0187202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2023] [Accepted: 01/23/2024] [Indexed: 02/16/2024] Open
Abstract
Stochastic differential equations (SDEs) are a powerful tool to model fluctuations and uncertainty in complex systems. Although numerical methods have been designed to simulate SDEs effectively, it is still problematic when numerical solutions may be negative, but application problems require positive simulations. To address this issue, we propose balanced implicit Patankar-Euler methods to ensure positive simulations of SDEs. Instead of considering the addition of balanced terms to explicit methods in existing balanced methods, we attempt the deletion of possible negative terms from the explicit methods to maintain positivity of numerical simulations. The designed balanced terms include negative-valued drift terms and potential negative diffusion terms. The proposed method successfully addresses the issue of divisions with very small denominators in our recently designed stochastic Patankar method. Stability analysis shows that the balanced implicit Patankar-Euler method has much better stability properties than our recently designed composite Patankar-Euler method. Four SDE systems are used to examine the effectiveness, accuracy, and convergence properties of balanced implicit Patankar-Euler methods. Numerical results suggest that the proposed balanced implicit Patankar-Euler method is an effective and efficient approach to ensure positive simulations when any appropriate stepsize is used in simulating SDEs of biological regulatory systems.
Collapse
Affiliation(s)
- Aimin Chen
- School of Mathematics and Statistics, Henan University, Kaifeng 475001, China
| | - Quanwei Ren
- College of Science, Henan University of Technology, Zhengzhou 450001, China
| | - Tianshou Zhou
- School of Mathematics and Statistics, Sun Yat-sen University, Guangzhong 510275, China
| | - Pamela Burrage
- School of Mathematical Sciences, Queensland University of Technology, Brisbane 4001, Australia
| | - Tianhai Tian
- School of Mathematics, Monash University, Clayton 3800, Australia
| | - Kevin Burrage
- School of Mathematical Sciences, Queensland University of Technology, Brisbane 4001, Australia
- Department of Computer Science, University of Oxford, Oxford OX1 3QD, United Kingdom
| |
Collapse
|
2
|
Stuhlmüller NCX, Farrokhzad F, Kuświk P, Stobiecki F, Urbaniak M, Akhundzada S, Ehresmann A, Fischer TM, de Las Heras D. Simultaneous and independent topological control of identical microparticles in non-periodic energy landscapes. Nat Commun 2023; 14:7517. [PMID: 37980403 PMCID: PMC10657436 DOI: 10.1038/s41467-023-43390-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2023] [Accepted: 11/07/2023] [Indexed: 11/20/2023] Open
Abstract
Topological protection ensures stability of information and particle transport against perturbations. We explore experimentally and computationally the topologically protected transport of magnetic colloids above spatially inhomogeneous magnetic patterns, revealing that transport complexity can be encoded in both the driving loop and the pattern. Complex patterns support intricate transport modes when the microparticles are subjected to simple time-periodic loops of a uniform magnetic field. We design a pattern featuring a topological defect that functions as an attractor or a repeller of microparticles, as well as a pattern that directs microparticles along a prescribed complex trajectory. Using simple patterns and complex loops, we simultaneously and independently control the motion of several identical microparticles differing only in their positions above the pattern. Combining complex patterns and complex loops we transport microparticles from unknown locations to predefined positions and then force them to follow arbitrarily complex trajectories concurrently. Our findings pave the way for new avenues in transport control and dynamic self-assembly in colloidal science.
Collapse
Affiliation(s)
- Nico C X Stuhlmüller
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440, Bayreuth, Germany
| | - Farzaneh Farrokhzad
- Experimatalphysik X, Physikalisches Institut, Universität Bayreuth, D-95440, Bayreuth, Germany
| | - Piotr Kuświk
- Institute of Molecular Physics, Polish Academy of Sciences, 60-179, Poznań, Poland
| | - Feliks Stobiecki
- Institute of Molecular Physics, Polish Academy of Sciences, 60-179, Poznań, Poland
| | - Maciej Urbaniak
- Institute of Molecular Physics, Polish Academy of Sciences, 60-179, Poznań, Poland
| | - Sapida Akhundzada
- Institute of Physics and Center for Interdisciplinary Nanostructure Science and Technology (CINSaT), University of Kassel, D-34132, Kassel, Germany
| | - Arno Ehresmann
- Institute of Physics and Center for Interdisciplinary Nanostructure Science and Technology (CINSaT), University of Kassel, D-34132, Kassel, Germany
| | - Thomas M Fischer
- Experimatalphysik X, Physikalisches Institut, Universität Bayreuth, D-95440, Bayreuth, Germany
| | - Daniel de Las Heras
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440, Bayreuth, Germany.
| |
Collapse
|
3
|
Sammüller F, Hermann S, de Las Heras D, Schmidt M. Noether-Constrained Correlations in Equilibrium Liquids. PHYSICAL REVIEW LETTERS 2023; 130:268203. [PMID: 37450808 DOI: 10.1103/physrevlett.130.268203] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/26/2023] [Accepted: 06/12/2023] [Indexed: 07/18/2023]
Abstract
Liquid structure carries deep imprints of an inherent thermal invariance against a spatial transformation of the underlying classical many-body Hamiltonian. At first order in the transformation field Noether's theorem yields the local force balance. Three distinct two-body correlation functions emerge at second order, namely the standard two-body density, the localized force-force correlation function, and the localized force gradient. An exact Noether sum rule interrelates these correlators. Simulations of Lennard-Jones, Yukawa, soft-sphere dipolar, Stockmayer, Gay-Berne and Weeks-Chandler-Andersen liquids, of monatomic water and of a colloidal gel former demonstrate the fundamental role in the characterization of spatial structure.
Collapse
Affiliation(s)
- Florian Sammüller
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Sophie Hermann
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Daniel de Las Heras
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Matthias Schmidt
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| |
Collapse
|
4
|
de Las Heras D, Zimmermann T, Sammüller F, Hermann S, Schmidt M. Perspective: How to overcome dynamical density functional theory. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2023; 35:271501. [PMID: 37023762 DOI: 10.1088/1361-648x/accb33] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/28/2023] [Accepted: 04/06/2023] [Indexed: 06/19/2023]
Abstract
We argue in favour of developing a comprehensive dynamical theory for rationalizing, predicting, designing, and machine learning nonequilibrium phenomena that occur in soft matter. To give guidance for navigating the theoretical and practical challenges that lie ahead, we discuss and exemplify the limitations of dynamical density functional theory (DDFT). Instead of the implied adiabatic sequence of equilibrium states that this approach provides as a makeshift for the true time evolution, we posit that the pending theoretical tasks lie in developing a systematic understanding of the dynamical functional relationships that govern the genuine nonequilibrium physics. While static density functional theory gives a comprehensive account of the equilibrium properties of many-body systems, we argue that power functional theory is the only present contender to shed similar insights into nonequilibrium dynamics, including the recognition and implementation of exact sum rules that result from the Noether theorem. As a demonstration of the power functional point of view, we consider an idealized steady sedimentation flow of the three-dimensional Lennard-Jones fluid and machine-learn the kinematic map from the mean motion to the internal force field. The trained model is capable of both predicting and designing the steady state dynamics universally for various target density modulations. This demonstrates the significant potential of using such techniques in nonequilibrium many-body physics and overcomes both the conceptual constraints of DDFT as well as the limited availability of its analytical functional approximations.
Collapse
Affiliation(s)
- Daniel de Las Heras
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Toni Zimmermann
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Florian Sammüller
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Sophie Hermann
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Matthias Schmidt
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| |
Collapse
|
5
|
Bryer AJ, Rey JS, Perilla JR. Performance efficient macromolecular mechanics via sub-nanometer shape based coarse graining. Nat Commun 2023; 14:2014. [PMID: 37037809 PMCID: PMC10086035 DOI: 10.1038/s41467-023-37801-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2022] [Accepted: 03/30/2023] [Indexed: 04/12/2023] Open
Abstract
Dimensionality reduction via coarse grain modeling is a valuable tool in biomolecular research. For large assemblies, ultra coarse models are often knowledge-based, relying on a priori information to parameterize models thus hindering general predictive capability. Here, we present substantial advances to the shape based coarse graining (SBCG) method, which we refer to as SBCG2. SBCG2 utilizes a revitalized formulation of the topology representing network which makes high-granularity modeling possible, preserving atomistic details that maintain assembly characteristics. Further, we present a method of granularity selection based on charge density Fourier Shell Correlation and have additionally developed a refinement method to optimize, adjust and validate high-granularity models. We demonstrate our approach with the conical HIV-1 capsid and heteromultimeric cofilin-2 bound actin filaments. Our approach is available in the Visual Molecular Dynamics (VMD) software suite, and employs a CHARMM-compatible Hamiltonian that enables high-performance simulation in the GPU-resident NAMD3 molecular dynamics engine.
Collapse
Affiliation(s)
- Alexander J Bryer
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE, 19716, USA
| | - Juan S Rey
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE, 19716, USA
| | - Juan R Perilla
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE, 19716, USA.
| |
Collapse
|
6
|
Renner J, Schmidt M, de Las Heras D. Reduced-variance orientational distribution functions from torque sampling. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2023; 35:235901. [PMID: 36974000 DOI: 10.1088/1361-648x/acc522] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/15/2022] [Accepted: 03/08/2023] [Indexed: 06/18/2023]
Abstract
We introduce a method to sample the orientational distribution function in computer simulations. The method is based on the exact torque balance equation for classical many-body systems of interacting anisotropic particles in equilibrium. Instead of the traditional counting of events, we reconstruct the orientational distribution function via an orientational integral of the torque acting on the particles. We test the torque sampling method in two- and three-dimensions, using both Langevin dynamics and overdamped Brownian dynamics, and with two interparticle interaction potentials. In all cases the torque sampling method produces profiles of the orientational distribution function with better accuracy than those obtained with the traditional counting method. The accuracy of the torque sampling method is independent of the bin size, and hence it is possible to resolve the orientational distribution function with arbitrarily small angular resolutions.
Collapse
Affiliation(s)
- Johannes Renner
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany
| | - Matthias Schmidt
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany
| | - Daniel de Las Heras
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany
| |
Collapse
|
7
|
Sammüller F, de Las Heras D, Schmidt M. Inhomogeneous steady shear dynamics of a three-body colloidal gel former. J Chem Phys 2023; 158:054908. [PMID: 36754804 DOI: 10.1063/5.0130655] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
We investigate the stationary flow of a colloidal gel under an inhomogeneous external shear force using adaptive Brownian dynamics simulations. The interparticle forces are derived from the Stillinger-Weber potential, where the three-body term is tuned to enable network formation and gelation in equilibrium. When subjected to the shear force field, the system develops remarkable modulations in the one-body density profile. Depending on the shear magnitude, particles accumulate either in quiescent regions or in the vicinity of maximum net flow, and we deduce this strong non-equilibrium response to be characteristic of the gel state. Studying the components of the internal force parallel and perpendicular to the flow direction reveals that the emerging flow and structure of the stationary state are driven by significant viscous and structural superadiabatic forces. Thereby, the magnitude and nature of the observed non-equilibrium phenomena differ from the corresponding behavior of simple fluids. We demonstrate that a simple power functional theory reproduces accurately the viscous force profile, giving a rationale of the complex dynamical behavior of the system.
Collapse
Affiliation(s)
- Florian Sammüller
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Daniel de Las Heras
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| | - Matthias Schmidt
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95447 Bayreuth, Germany
| |
Collapse
|
8
|
Muñiz‐Chicharro A, Votapka LW, Amaro RE, Wade RC. Brownian dynamics simulations of biomolecular diffusional association processes. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1649] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Abraham Muñiz‐Chicharro
- Molecular and Cellular Modeling Group Heidelberg Institute for Theoretical Studies (HITS) Heidelberg Germany
- Faculty of Biosciences and Heidelberg Graduate School of Mathematical and Computational Methods for the Sciences (HGS MathComp) Heidelberg University Heidelberg Germany
| | | | | | - Rebecca C. Wade
- Molecular and Cellular Modeling Group Heidelberg Institute for Theoretical Studies (HITS) Heidelberg Germany
- Center for Molecular Biology (ZMBH), DKFZ‐ZMBH Alliance, and Interdisciplinary Center for Scientific Computing (IWR) Heidelberg University Heidelberg Germany
| |
Collapse
|
9
|
Stuhlmüller NCX, Fischer TM, de Las Heras D. Colloidal transport in twisted lattices of optical tweezers. Phys Rev E 2022; 106:034601. [PMID: 36266822 DOI: 10.1103/physreve.106.034601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2022] [Accepted: 08/09/2022] [Indexed: 06/16/2023]
Abstract
We simulate the transport of colloidal particles driven by a static and homogeneous drift force, and subject to the optical potential created by two lattices of optical tweezers. The lattices of optical tweezers are parallel to each other, shifted, and rotated by a twist angle. Due to a negative interference between the potential of the two lattices, flat channels appear in the total optical potential. At specific twist angles, known as magic angles, the flat channels percolate the entire system and the colloidal particles can then be transported using a weak external drift force. We characterize the transport in both square and hexagonal lattices of twisted optical tweezers.
Collapse
Affiliation(s)
- Nico C X Stuhlmüller
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany
| | - Thomas M Fischer
- Experimatalphysik X, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany
| | - Daniel de Las Heras
- Theoretische Physik II, Physikalisches Institut, Universität Bayreuth, D-95440 Bayreuth, Germany
| |
Collapse
|