-
Biohybrid membrane formation by directed insertion of Aquaporin into a solid-state nanopore
Authors:
Francois Sicard,
A. Ozgur Yazaydin
Abstract:
Technical challenges in molecule sensing and chemical detection have created an increasing demand for transformative materials with high sensitivity and specificity. Biohybrid nanopores have attracted growing interest as they can ideally combine the durability of solid-state nanopores with the precise structure of biological nanopores. Particular care must be taken to control how biological nanopo…
▽ More
Technical challenges in molecule sensing and chemical detection have created an increasing demand for transformative materials with high sensitivity and specificity. Biohybrid nanopores have attracted growing interest as they can ideally combine the durability of solid-state nanopores with the precise structure of biological nanopores. Particular care must be taken to control how biological nanopores adapt to their surroundings once in contact with the solid-state nanopore. Two major challenges are to precisely control this adaptability under dynamic conditions and provide predesigned functionalities that can be manipulated for engineering applications. Here, we report on the computational design of a distinctive class of biohybrid active membrane layer, built from the directed insertion of an aquaporin-incorporated lipid shell into a silica nanopore. First, we describe in detail the mechanisms at play in the insertion of the biological membrane into the solid-state nanopore. Then we analyze the structural stability of the system and demonstrate that its water permeability is comparable to the one measured in the biological environment. Finally, we discuss how the technology implemented could be applicable to environmental and biomedical applications, such as water desalination and drug discovery, where targeting and controlled permeation of small molecules must be efficiently addressed.
△ Less
Submitted 4 April, 2022;
originally announced April 2022.
-
Computational design of armored nanodroplets as nanocarriers for encapsulation and release under flow conditions
Authors:
François Sicard,
Jhoan Toro-Mendoza
Abstract:
Nanocarriers are nanosized materials commonly used for targeted-oriented delivery of active compounds, including antimicrobials and small-molecular drugs. They equally represent fundamental and engineering challenges since sophisticated nanocarriers must show adequate structure, stability, and function in complex ambients. Here, we report on the computational design of a distinctive class of nanoc…
▽ More
Nanocarriers are nanosized materials commonly used for targeted-oriented delivery of active compounds, including antimicrobials and small-molecular drugs. They equally represent fundamental and engineering challenges since sophisticated nanocarriers must show adequate structure, stability, and function in complex ambients. Here, we report on the computational design of a distinctive class of nanocarriers, built from buckled armored nanodroplets, able to selectively encapsulate or release a probe load under specific flow conditions. Mesoscopic simulations offer detailed insight into the interplay between the characteristics of laden surface coverage and evolution of the droplet morphology. First, we describe in detail the formation of \textit{pocket-like} structures in Pickering emulsion nanodroplets and their stability under external flow. Then we use that knowledge to test the capacity of these emulsion-based pockets to yield flow-assisted encapsulation or expulsion of a probe load. Finally, the rheological properties of our model carrier are put into perspective with those of delivery systems employed in pharmaceutical and cosmetic technology.
△ Less
Submitted 18 January, 2021;
originally announced January 2021.
-
Role of Structural Rigidity and Collective Behaviour in the Molecular Design of Gas Hydrates Anti-Agglomerants
Authors:
François Sicard,
Alberto Striolo
Abstract:
Antiagglomerants (AAs) are surface active molecules widely used in the rubber and petroleum industry, among others. In the petroleum industry, it is believed that AAs strongly adsorb to the surface of hydrate particles to prevent the growth of clathrate hydrate within oil pipelines. Small changes in their molecular structures can strongly affect the thermodynamic and kinetic stability of the syste…
▽ More
Antiagglomerants (AAs) are surface active molecules widely used in the rubber and petroleum industry, among others. In the petroleum industry, it is believed that AAs strongly adsorb to the surface of hydrate particles to prevent the growth of clathrate hydrate within oil pipelines. Small changes in their molecular structures can strongly affect the thermodynamic and kinetic stability of the system as a whole. Here we employ molecular dynamics simulations to study the interplay between the modification of the molecular structure, rigidity and collective effects of AAs designed to prevent hydrate agglomeration in the conditions encountered in rocking cell experiments. The AAs are surface-active compounds with a complex hydrophilic head and three hydrophobic tails whose structural rigidity is enhanced with the attachement of a simple aromatic group. We observe that the aromatic group can positively or negatively affect the performance of the AAs, depending on its location along the hydrophobic tail. Our approach is based on first quantifying the molecular mechanisms responsible for the macroscopic performance. Although the mechanisms at play depend on the application, the methodology implemented could be applicable to other high-tech industries, where the agglomeration of small particles must be controlled.
△ Less
Submitted 17 December, 2020;
originally announced December 2020.
-
Assessing Position-Dependent Diffusion from Biased Simulations and Markov State Model Analysis
Authors:
Francois Sicard,
Vladimir Koskin,
Alessia Annibale,
Edina Rosta
Abstract:
A variety of enhanced statistical and numerical methods are now routinely used to extract comprehensible and relevant thermodynamic information from the vast amount of complex, high-dimensional data obtained from intensive molecular simulations. The characterization of kinetic properties, such as diffusion coefficients, of molecular systems with significantly high energy barriers, on the other han…
▽ More
A variety of enhanced statistical and numerical methods are now routinely used to extract comprehensible and relevant thermodynamic information from the vast amount of complex, high-dimensional data obtained from intensive molecular simulations. The characterization of kinetic properties, such as diffusion coefficients, of molecular systems with significantly high energy barriers, on the other hand, has received less attention. Among others, Markov state models, in which the long-time statistical dynamics of a system is approximated by a Markov chain on a discrete partition of configuration space, have seen widespread use in recent years, with the aim of tackling these fundamental issues. Here we propose a general, automatic method to assess multidimensional position-dependent diffusion coefficients within the framework of Markovian stochastic processes and Kramers-Moyal expansion. We apply the formalism to one- and two-dimensional analytic potentials and data from explicit solvent molecular dynamics simulations, including the water-mediated conformations of alanine dipeptide. Importantly, the developed algortihm presents significant improvement compared to standard methods when the transport of solute across three-dimensional heterogeneous porous media is studied, for example, the prediction of membrane permeation of drug molecules.
△ Less
Submitted 13 October, 2020;
originally announced October 2020.
-
Nanoparticles actively fragment armored droplets
Authors:
François Sicard,
Jhoan Toro-Mendoza,
Alberto Striolo
Abstract:
Understanding the complexity of fragmentation processes is essential for regulating intercellular communication in mechanistic biology and developing novel bottom-up approaches in a large range of multiphase flow processes. In this context, self-fragmentation proceeds without any external mechanical energy input allowing one to create efficiently micro- and nanodroplets. Here we examine self-fragm…
▽ More
Understanding the complexity of fragmentation processes is essential for regulating intercellular communication in mechanistic biology and developing novel bottom-up approaches in a large range of multiphase flow processes. In this context, self-fragmentation proceeds without any external mechanical energy input allowing one to create efficiently micro- and nanodroplets. Here we examine self-fragmentation in emulsion nanodroplets stabilized by solid particles with different surface features. Mesoscopic modelling and accelerated dynamics simulations allow us to overcome the limitations of atomistic simulations and offer detailed insight into the interplay between the evolution of the droplet shape and the particle finite-sized effects at the interface. We show that finite-sized nanoparticles play an active role in the necking breakup, behaving like nano-scale razors, and affect strongly the thermodynamic properties of the system. The role played by the particles during self-fragmentation might be of relevance to multifunctional biomaterial design and tuning of signaling pathways in mechanistic biology.
△ Less
Submitted 7 May, 2019;
originally announced May 2019.
-
Control of DNA denaturation bubble nucleation to advance nano-biosensing
Authors:
François Sicard,
Nicolas Destainville,
Philippe Rousseau,
Catherine Tardin,
Manoel Manghi
Abstract:
In the demanding biosensing environment, improving selection efficiency strategies has become an issue of great significance. DNA minicircles containing between 200 and 400 base-pairs, also named microDNA, are representative of the supercoiled DNA loops found in nature. Their short size makes them extremely susceptible to writhe and twist, which is known to play a central role in DNA denaturation.…
▽ More
In the demanding biosensing environment, improving selection efficiency strategies has become an issue of great significance. DNA minicircles containing between 200 and 400 base-pairs, also named microDNA, are representative of the supercoiled DNA loops found in nature. Their short size makes them extremely susceptible to writhe and twist, which is known to play a central role in DNA denaturation. We investigate minicircle lengths and superhelical densities that induce DNA denaturation bubbles of nanometer size and control well-defined long-life. Mesoscopic modeling and accelerated dynamics simulations allow us to study accurately the thermodynamic and dynamical properties associated with the nucleation and closure mechanisms of long-lived denaturation bubbles. Our results pave the way for new types of DNA biosensors with enhanced selectivity for specific DNA binding proteins.
△ Less
Submitted 11 May, 2018;
originally announced May 2018.
-
Emergent Properties of Antiagglomerant Films Control Methane Transport: Implications for Hydrate Management
Authors:
François Sicard,
Tai Bui,
Deepak Monteiro,
Qiang Lan,
Mark Ceglio,
Charlotte Burress,
Alberto Striolo
Abstract:
The relation between collective properties and performance of antiagglomerants (AAs) used in hydrate management is handled using molecular dynamics simulations and enhanced sampling techniques. A thin film of AAs adsorbed at the interface between one flat sII methane hydrate substrate and a fluid hydrocarbon mixture containing methane and n-dodecane is studied. The AA considered is a surface-activ…
▽ More
The relation between collective properties and performance of antiagglomerants (AAs) used in hydrate management is handled using molecular dynamics simulations and enhanced sampling techniques. A thin film of AAs adsorbed at the interface between one flat sII methane hydrate substrate and a fluid hydrocarbon mixture containing methane and n-dodecane is studied. The AA considered is a surface-active compound with a complex hydrophilic head that contains both amide and tertiary ammonium cation groups and hydrophobic tails. At sufficiently high AA density, the interplay between the surfactant layer and the liquid hydrocarbon excludes methane from the interfacial region. In this scenario, we combine metadynamics and umbrella sampling frameworks to study accurately the free-energy landscape and the equilibrium rates associated with the transport of one methane molecule across the AA film. We observe that local configurational changes of the liquid hydrocarbon packed within the AA film are associated with high free-energy barriers for methane transport. The time scales estimated for the transport of methane across the AA film can be, in some cases, comparable to those reported in the literature for the growth of the hydrates, suggesting that one possible mechanism by which AAs delay the formation of hydrate plugs could be providing a barrier to methane transport. Considering the interplay between structural design and collective properties of AAs might be of relevance to improve their performance in flow assurance.
△ Less
Submitted 4 April, 2018;
originally announced April 2018.
-
Computing Transition Rates for Rare Event: When Kramers Theory meets Free Energy Landscape
Authors:
François Sicard
Abstract:
Computing reactive trajectories and free energy (FE) landscapes associated to rare event kinetics is key to understanding the dynamics of complex systems. The analysis of the FE surface on which the underlying dynamics takes place has become central to compute transition rates. In the overdamped limit, most often encountered in biophysics and soft condensed matter, the Kramers' Theory (KT) has pro…
▽ More
Computing reactive trajectories and free energy (FE) landscapes associated to rare event kinetics is key to understanding the dynamics of complex systems. The analysis of the FE surface on which the underlying dynamics takes place has become central to compute transition rates. In the overdamped limit, most often encountered in biophysics and soft condensed matter, the Kramers' Theory (KT) has proved to be quite successful in recovering correct kinetics. However, the additional calculation to obtain rate constants in complex systems where configurational entropy is competing with energy is still challenging conceptually and computationally. Building on KT and the metadynamics framework, the rate is expressed in terms of the height of the FE barrier measured along the minimum FE path and an auxiliary measure of the configurational entropy. We apply the formalism to two different problems where our approach shows good agreement with simulations and experiments and can present significant improvement over the standard KT.
△ Less
Submitted 19 March, 2018; v1 submitted 9 March, 2018;
originally announced March 2018.
-
Buckling in Armored Droplets
Authors:
François Sicard,
Alberto Striolo
Abstract:
The issue of the buckling mechanism in droplets stabilized by solid particles (armored droplets) is tackled at a mesoscopic level using dissipative particle dynamics simulations. We consider spherical water droplet in a decane solvent coated with nanoparticle monolayers of two different types: Janus and homogeneous. The chosen particles yield comparable initial three-phase contact angles, chosen t…
▽ More
The issue of the buckling mechanism in droplets stabilized by solid particles (armored droplets) is tackled at a mesoscopic level using dissipative particle dynamics simulations. We consider spherical water droplet in a decane solvent coated with nanoparticle monolayers of two different types: Janus and homogeneous. The chosen particles yield comparable initial three-phase contact angles, chosen to maximize the adsorption energy at the interface. We study the interplay between the evolution of droplet shape, layering of the particles, and their distribution at the interface when the volume of the droplets is reduced. We show that Janus particles affect strongly the shape of the droplet with the formation of a crater-like depression. This evolution is actively controlled by a close-packed particle monolayer at the curved interface. On the contrary, homogeneous particles follow passively the volume reduction of the droplet, whose shape does not deviate too much from spherical, even when a nanoparticle monolayer/bilayer transition is detected at the interface. We discuss how these buckled armored droplets might be of relevance in various applications including potential drug delivery systems and biomimetic design of functional surfaces.
△ Less
Submitted 8 February, 2017;
originally announced February 2017.
-
Numerical analysis of Pickering emulsion stability: insights from ABMD simulations
Authors:
François Sicard,
Alberto Striolo
Abstract:
The issue of the stability of Pickering emulsions is tackled at a mesoscopic level using dissipative particle dynamics simulations within the Adiabatic Biased Molecular Dynamics framework. We consider the early stage of the coalescence process between two spherical water droplets in decane solvent. The droplets are stabilized by Janus nanoparticles of different shapes (spherical and ellipsoidal) w…
▽ More
The issue of the stability of Pickering emulsions is tackled at a mesoscopic level using dissipative particle dynamics simulations within the Adiabatic Biased Molecular Dynamics framework. We consider the early stage of the coalescence process between two spherical water droplets in decane solvent. The droplets are stabilized by Janus nanoparticles of different shapes (spherical and ellipsoidal) with different three-phase contact angles. Given a sufficiently dense layer of particles on the droplets, we show that the stabilization mechanism strongly depends on the collision speed. This is consistent with a coalescence mechanism governed by the rheology of the interfacial region. When the system is forced to coalesce sufficiently slowly, we investigate at a mesoscopic level how the ability of the nanoparticles to stabilize Pickering emulsions is discriminated by nanoparticle mobility and the associated caging effect. These properties are both related to the interparticle interaction and the hydrodynamic resistance in the liquid film between the approaching interfaces.
△ Less
Submitted 23 March, 2016;
originally announced March 2016.
-
DNA denaturation bubbles: free-energy landscape and nucleation/closure rates
Authors:
François Sicard,
Nicolas Destainville,
Manoel Manghi
Abstract:
The issue of the nucleation and slow closure mechanisms of non superhelical stress-induced denaturation bubbles in DNA is tackled using coarse-grained MetaDynamics and Brownian simulations. A minimal mesoscopic model is used where the double helix is made of two interacting bead-spring rotating strands with a prescribed torsional modulus in the duplex state. We demonstrate that timescales for the…
▽ More
The issue of the nucleation and slow closure mechanisms of non superhelical stress-induced denaturation bubbles in DNA is tackled using coarse-grained MetaDynamics and Brownian simulations. A minimal mesoscopic model is used where the double helix is made of two interacting bead-spring rotating strands with a prescribed torsional modulus in the duplex state. We demonstrate that timescales for the nucleation (resp. closure) of an approximately 10 base-pair bubble, in agreement with experiments, are associated with the crossing of a free-energy barrier of $22~k_{\rm B}T$ (resp. $13~k_{\rm B}T$) at room temperature $T$. MetaDynamics allows us to reconstruct accurately the free-energy landscape, to show that the free-energy barriers come from the difference in torsional energy between the bubble and duplex states, and thus to highlight the limiting step, a collective twisting, that controls the nucleation/closure mechanism, and to access opening time scales on the millisecond range. Contrary to small breathing bubbles, these more than 4~base-pair bubbles are of biological relevance, for example when a preexisting state of denaturation is required by specific DNA-binding proteins.
△ Less
Submitted 12 January, 2015; v1 submitted 15 May, 2014;
originally announced May 2014.
-
Scaling quasi-stationary states in long range systems with dissipation
Authors:
Michael Joyce,
Jules Morand,
François Sicard,
Pascal Viot
Abstract:
Hamiltonian systems with long-range interactions give rise to long lived out of equilibrium macroscopic states, so-called quasi-stationary states. We show here that, in a suitably generalized form, this result remains valid for many such systems in the presence of dissipation. Using an appropriate mean-field kinetic description, we show that models with dissipation due to a viscous damping or due…
▽ More
Hamiltonian systems with long-range interactions give rise to long lived out of equilibrium macroscopic states, so-called quasi-stationary states. We show here that, in a suitably generalized form, this result remains valid for many such systems in the presence of dissipation. Using an appropriate mean-field kinetic description, we show that models with dissipation due to a viscous damping or due to inelastic collisions admit "scaling quasi-stationary states", i.e., states which are quasi-stationary in rescaled variables. A numerical study of one dimensional self-gravitating systems confirms both the relevance of these solutions, and gives indications of their regime of validity in line with theoretical predictions. We underline that the velocity distributions never show any tendency to evolve towards a Maxwell-Boltzmann form.
△ Less
Submitted 15 January, 2014; v1 submitted 13 September, 2013;
originally announced September 2013.
-
Exponents of non-linear clustering in scale-free one dimensional cosmological simulations
Authors:
David Benhaiem,
Michael Joyce,
François Sicard
Abstract:
One dimensional versions of cosmological N-body simulations have been shown to share many qualitative behaviours of the three dimensional problem. They can resolve a large range of time and length scales, and admit exact numerical integration. We use such models to study how non-linear clustering depends on initial conditions and cosmology. More specifically, we consider a family of models which,…
▽ More
One dimensional versions of cosmological N-body simulations have been shown to share many qualitative behaviours of the three dimensional problem. They can resolve a large range of time and length scales, and admit exact numerical integration. We use such models to study how non-linear clustering depends on initial conditions and cosmology. More specifically, we consider a family of models which, like the 3D EdS model, lead for power-law initial conditions to self-similar clustering characterized in the strongly non-linear regime by power-law behaviour of the two point correlation function. We study how the corresponding exponent γdepends on the initial conditions, characterized by the exponent n of the power spectrum of initial fluctuations, and on a single parameter κcontrolling the rate of expansion. The space of initial conditions/cosmology divides very clearly into two parts: (1) a region in which γdepends strongly on both n and κand where it agrees very well with a simple generalisation of the so-called stable clustering hypothesis in three dimensions, and (2) a region in which γis more or less independent of both the spectrum and the expansion of the universe. We explain the observed location of the boundary in (n, κ) space dividing the "stable clustering" region from the "universal" region. We compare and contrast our findings to results in three dimensions, and discuss in particular the light they may throw on the question of "universality" of non-linear clustering in this context.
△ Less
Submitted 6 September, 2013; v1 submitted 28 November, 2012;
originally announced November 2012.
-
Reconstructing the free-energy landscape of Met-enkephalin using dihedral Principal Component Analysis and Well-tempered Metadynamics
Authors:
Francois Sicard,
Patrick Senet
Abstract:
Well-Tempered Metadynamics (WTmetaD) is an efficient method to enhance the reconstruction of the free-energy surface of proteins. WTmetaD guarantees a faster convergence in the long time limit in comparison with the standard metadynamics. It still suffers however from the same limitation, i.e. the non trivial choice of pertinent collective variables (CVs). To circumvent this problem, we couple WTm…
▽ More
Well-Tempered Metadynamics (WTmetaD) is an efficient method to enhance the reconstruction of the free-energy surface of proteins. WTmetaD guarantees a faster convergence in the long time limit in comparison with the standard metadynamics. It still suffers however from the same limitation, i.e. the non trivial choice of pertinent collective variables (CVs). To circumvent this problem, we couple WTmetaD with a set of CVs generated from a dihedral Principal Component Analysis (dPCA) on the Ramachadran dihedral angles describing the backbone structure of the protein. The dPCA provides a generic method to extract relevant CVs built from internal coordinates. We illustrate the robustness of this method in the case of the small and very diffusive Metenkephalin pentapeptide, and highlight a criterion to limit the number of CVs necessary to biased the metadynamics simulation. The free-energy landscape (FEL) of Met-enkephalin built on CVs generated from dPCA is found rugged compared with the FEL built on CVs extracted from PCA of the Cartesian coordinates of the atoms.
△ Less
Submitted 12 November, 2012;
originally announced November 2012.
-
Non-linear gravitational clustering of cold matter in an expanding universe: indications from 1D toy models
Authors:
Michael Joyce,
François Sicard
Abstract:
Studies of a class of infinite one dimensional self-gravitating systems have highlighted that, on the one hand, the spatial clustering which develops may have scale invariant (fractal) properties, and, on the other, that they display "self-similar" properties in their temporal evolution. The relevance of these results to three dimensional cosmological simulations has remained unclear. We show here…
▽ More
Studies of a class of infinite one dimensional self-gravitating systems have highlighted that, on the one hand, the spatial clustering which develops may have scale invariant (fractal) properties, and, on the other, that they display "self-similar" properties in their temporal evolution. The relevance of these results to three dimensional cosmological simulations has remained unclear. We show here that the measured exponents characterizing the scale-invariant non-linear clustering are in excellent agreement with those derived from an appropriately generalized "stable-clustering" hypothesis. Further an analysis in terms of "halos" selected with a friend-of-friend algorithm reveals that such structures are, statistically, virialized across the range of scales corresponding to scale-invariance. Thus the strongly non-linear clustering in these models is accurately described as a virialized fractal structure, very much in line with the "clustering hierarchy" which Peebles originally envisaged qualitatively as associated with stable clustering. If transposed to three dimensions these results would imply, notably, that cold dark matter halos (or even subhalos) are 1) not well modeled as smooth objects, and 2) that the supposed "universality" of their profiles is, like apparent smoothness, an artefact of poor numerical resolution.
△ Less
Submitted 2 March, 2011; v1 submitted 7 December, 2010;
originally announced December 2010.
-
A dynamical classification of the range of pair interactions
Authors:
Andrea Gabrielli,
Michael Joyce,
Bruno Marcos,
Francois Sicard
Abstract:
We formalize a classification of pair interactions based on the convergence properties of the {\it forces} acting on particles as a function of system size. We do so by considering the behavior of the probability distribution function (PDF) P(F) of the force field F in a particle distribution in the limit that the size of the system is taken to infinity at constant particle density, i.e., in the "…
▽ More
We formalize a classification of pair interactions based on the convergence properties of the {\it forces} acting on particles as a function of system size. We do so by considering the behavior of the probability distribution function (PDF) P(F) of the force field F in a particle distribution in the limit that the size of the system is taken to infinity at constant particle density, i.e., in the "usual" thermodynamic limit. For a pair interaction potential V(r) with V(r) \rightarrow \infty) \sim 1/r^a defining a {\it bounded} pair force, we show that P(F) converges continuously to a well-defined and rapidly decreasing PDF if and only if the {\it pair force} is absolutely integrable, i.e., for a > d-1, where d is the spatial dimension. We refer to this case as {\it dynamically short-range}, because the dominant contribution to the force on a typical particle in this limit arises from particles in a finite neighborhood around it. For the {\it dynamically long-range} case, i.e., a \leq d-1, on the other hand, the dominant contribution to the force comes from the mean field due to the bulk, which becomes undefined in this limit. We discuss also how, for a \leq d-1 (and notably, for the case of gravity, a=d-2) P(F) may, in some cases, be defined in a weaker sense. This involves a regularization of the force summation which is generalization of the procedure employed to define gravitational forces in an infinite static homogeneous universe. We explain that the relevant classification in this context is, however, that which divides pair forces with a > d-2 (or a < d-2), for which the PDF of the {\it difference in forces} is defined (or not defined) in the infinite system limit, without any regularization. In the former case dynamics can, as for the (marginal) case of gravity, be defined consistently in an infinite uniform system.
△ Less
Submitted 2 November, 2010; v1 submitted 29 March, 2010;
originally announced March 2010.
-
1-d gravity in infinite point distributions
Authors:
Andrea Gabrielli,
Michael Joyce,
Francois Sicard
Abstract:
The dynamics of infinite, asymptotically uniform, distributions of self-gravitating particles in one spatial dimension provides a simple toy model for the analogous three dimensional problem. We focus here on a limitation of such models as treated so far in the literature: the force, as it has been specified, is well defined in infinite point distributions only if there is a centre of symmetry (…
▽ More
The dynamics of infinite, asymptotically uniform, distributions of self-gravitating particles in one spatial dimension provides a simple toy model for the analogous three dimensional problem. We focus here on a limitation of such models as treated so far in the literature: the force, as it has been specified, is well defined in infinite point distributions only if there is a centre of symmetry (i.e. the definition requires explicitly the breaking of statistical translational invariance). The problem arises because naive background subtraction (due to expansion, or by "Jeans' swindle" for the static case), applied as in three dimensions, leaves an unregulated contribution to the force due to surface mass fluctuations. Following a discussion by Kiessling, we show that the problem may be resolved by defining the force in infinite point distributions as the limit of an exponentially screened pair interaction. We show that this prescription gives a well defined (finite) force acting on particles in a class of perturbed infinite lattices, which are the point processes relevant to cosmological N-body simulations. For identical particles the dynamics of the simplest toy model is equivalent to that of an infinite set of points with inverted harmonic oscillator potentials which bounce elastically when they collide. We discuss previous results in the literature, and present new results for the specific case of this simplest (static) model starting from "shuffled lattice" initial conditions. These show qualitative properties (notably its "self-similarity") of the evolution very similar to those in the analogous simulations in three dimensions, which in turn resemble those in the expanding universe.
△ Less
Submitted 28 September, 2009; v1 submitted 22 December, 2008;
originally announced December 2008.
